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SOIL  EM  PROPERTIES  AND  METAL  DETECTOR  PERFORMANCE: 
THEORY  AND  MEASUREMENT 

Guy  Cross 

Terrascan  Geophysics,  Vancouver,  Canada 


Abstract: 

Metal  detectors  are  unquestionably  the  workhorses  of  humanitarian  demining.  Al¬ 
though  hybrid  dual-sensor  systems  incorporating  a  ground  penetrating  radar  (GPR) 
have  recently  been  introduced  and  are  gaining  acceptance,  a  standard  electromag¬ 
netic  induction  (EMI)  sensor  remains  the  primary  mode  of  detection.  Despite  on¬ 
going  evolution  and  rehnement  of  metal  detector  technology,  however,  the  practical 
performance  of  both  continuous  wave  (FDEM)  and  pulse  induction  (TDEM)  sys¬ 
tems  continues  to  be  restricted  by  so-called  “soil  noise”  effects.  Generally,  “problem” 
or  “difficult”  soils  reduce  signal  to  noise  ratio  and  increase  the  false- detection  rate. 
At  certain  locations,  the  soil  effect  is  so  severe  as  to  render  a  given  metal  detector 
practically  inoperable. 

Although  soil  interference  is  well  established  and  widely  appreciated,  the  precise  ori- 
gin(s)  of  related  noise  has  remained  elusive  as  evidenced  by  the  range  of  common 
descriptions  for  problem  soils,  including  “lateritic  soil”,  “mineralized  soil”,  “conduct¬ 
ing  soil”  and  “magnetic  soil” .  The  absence  of  clarity  and  definitive  understanding  has 
impeded  specification  of  appropriate  methods  and  procedures  for  related  assessment 
of  soil  conditions  and  for  predicting  the  nature  and  extent  of  corresponding  influence 
on  a  given  metal  detector  technology. 

To  address  the  foregoing  situation,  the  study  reported  herein  has  been  part  of  a 
broader  international  effort  to  define  and  quantify  the  effects  of  soil  electromagnetic 
properties  on  the  operation  and  performance  of  metal  detectors  and  allied  (GPR) 
sensors.  In  particular,  the  present  study  is  restricted  to  EMI  metal  detectors  and  our 
aim  has  been  to  characterize  the  nature  and  extent  of  influence  due  to  specific  soil 
electromagnetic  properties  and  related  frequency  dispersion. 

Our  report  begins  with  a  brief  review  of  relevant  electromagnetic  theory  and  sub¬ 
sequently  explores  the  basic  physics  of  EMI  metal  detectors  and  associated  implica¬ 
tions  for  the  relative  significance  of  specific  soil  EM  parameters.  Particular  attention 
is  focused  on  TDEM  sensors  and  related  influence  of  dispersive  soil  electromagnetic 
properties.  The  balance  of  the  report  is  devoted  to  the  theory  and  practice  of  mea¬ 
suring  low-frequency  (100  Hz  -  100  kHz)  soil  electromagnetic  properties.  A  range 
of  practical  instrumentation  and  procedures  are  investigated  with  reference  to  three 
DRDG  reference  soils. 
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1.0  BACKGROUND 


Previous  and  ongoing  conflicts  have  left  approximately  100  million  unexploded  land¬ 
mines  within  the  soils  of  more  than  fifty  countries  worldwide.  Unfortunately,  devel¬ 
oping  nations  that  are  least  prepared  to  cope  with  the  problem  are  disproportionately 
affected.  Vast  tracts  of  land  have  been  effectively  removed  from  agriculture  or  other 
productive  utilization,  and  the  human  cost  is  intolerable.  The  United  Nations  esti¬ 
mates  that  approximately  25,000  civilians  are  killed  or  maimed  by  landmines  each 
year. 

Growing  awareness  of  the  landmine  problem  over  the  past  decade  has  fostered  a  sub¬ 
stantial  international  response,  with  governments  and  non-governmental  organizations 
currently  sponsoring  active  humanitarian  demining  programs  in  virtually  every  mine- 
affected  nation.  However,  demining  is  a  costly  and  time-consuming  process  and  it  is 
currently  estimated  that  for  every  mine  that  is  removed  or  neutralized,  approximately 
20  more  are  laid.  Clearly  there  is  an  urgent  requirement  for  new  and  effective  mine 
detection  technologies  and,  concurrently,  for  the  refinement  of  well-established  and 
proven  methodologies. 

In  particular,  despite  ongoing  development  of  novel  and  alternative  sensor  technolo¬ 
gies,  metal  detectors  remain  the  workhorses  of  humanitarian  demining.  Crucially, 
however,  the  performance  of  both  continuous  wave  and  pulse  induction  metal  detec¬ 
tors  can  be  severely  restricted  by  so-called  “soil-noise”  effects.  Generally,  “problem 
soils”  reduce  detector  sensitivity  and  increase  the  false-detection  rate.  In  certain  lo¬ 
cations,  the  soil-effect  is  so  strong  as  to  render  a  given  detector  practically  inoperable. 

Although  the  adverse  influence  of  certain  soils  are  widely  recognized  and  acknowl¬ 
edged  by  the  demining  community,  there  has  only  recently  been  a  concerted  effort 
to  understand  and  establish  the  responsible  soil  properties  and  underlying  physical 
mechanisms  (Das,  et  ah,  2002,  2003;  ITEP,  2002).  The  current  study  is  a  part  this 
ongoing  effort  to  define  and  quantify  the  influence  of  soil  electromagnetic  properties 
on  the  operation  of  metal  detectors.  Specifically,  our  focus  has  been  to  identify  key 
soil  parameters  and  to  investigate  appropriate  technologies  and  procedures  for  related 
soils  characterization. 

We  begin  with  a  brief  review  of  relevant  electromagnetic  theory  and  subsequently 
explore  the  basic  physics  of  metal  detectors  and  related  implications  for  the  relative 
significance  of  specific  soil  EM  parameters.  The  remainder  of  the  report  is  devoted 
to  the  theory  and  practice  of  measuring  low-frequency  (100  Hz  -  100  kHz)  soil  elec¬ 
tromagnetic  properties.  The  connection  between  measured  soil  electromagnetic  prop¬ 
erties  and  related  physical,  chemical  and  mineralogical  characteristics  of  soils  is  not 
addressed  here  (see  Cross,  2000). 
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2.0  ELECTROMAGNETIC  THEORY  — 

ROLE  OF  SOIL  EM  PROPERTIES 

In  connection  with  a  previons  project  (DRES  CR  2000-091),  Cross  (2000)  reviewed 
the  natnre  and  description  of  soil  electromagnetic  properties  as  relevant  to  the  per¬ 
formance  of  gronnd  penetrating  radar  (GPR)  sensors.  Analysis  emphasized  the 
freqnency-dependent,  complex  natnre  of  electromagnetic  parameters  and  associated 
effective  soil  properties. 

Beginning  with  Maxwell’s  eqnations 


at 

(1) 

VxH  =  J  +  ^, 
at 

(2) 

and  assuming  constitutive  relations 

B  =  //H, 

(3) 

D  =  eE, 

(4) 

J  =  ctE, 

(5) 

relating  electric  E  and  magnetic  H  helds,  associated  electric  displacement  D,  magnetic 
induction  B  and  the  electric  current  density  J,  the  governing  relations 


V2E  =  72e  (6) 

=  (7) 

were  derived  for  harmonically-varying  electromagnetic  helds  within  a  uniform,  linear 
isotropic  soil.  The  associated  propagation  constant 

7  =  [i  277  fn  {a  +  i  27r/e)] 


=  27r/ 


(8) 


encompasses  the  inhuence  of  soil  electromagnetic  properties,  including  the  magnetic 
permeability  fx,  electrical  permittivity  e  and  electrical  conductivity  a.  In  addition, 
it  was  observed  that  these  parameters  are  ordinarily  frequency  dependent  and,  con¬ 
sequently,  complex  valued.  Moreover,  because  observed  dispersion  is  in  each  case 
associated  with  some  form  of  polarization  and  alternating  relaxation,  it  was  shown 
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that  the  frequency  dependence  of  soil  electromagnetic  parameters  can  be  usefully  mod¬ 
elled  by  appropriate  generalizations  of  the  well-known  Debye  spectrum  as  depicted  in 
Figure  1. 

It  was  also  emphasized  that  in  the  practical  process  of  measuring  soil  electromagnetic 
parameters,  and  to  the  extent  that  soil  properties  influence  the  response  of  a  ground 
radar  (or  metal  detector),  it  is  impossible  to  separate  the  combined  influence  of  mul¬ 
tiple  properties  contributing  in-phase.  In  particular,  dehning  complex  parameters 


a  =  o'  ia' 


/  ■  // 
e  =  e  —  le 


/  ■  // 
//  =  //  -i^i  , 


it  was  demonstrated  that  the  propagation  constant  can  be  rewritten  as 

.1/2 


7  =  27r/ 


lie  I  i 


2ir/e 


1 


(9) 

(10) 

(11) 


(12) 


in  terms  of  a  real-valued  effective  conductivity 

a  =  a  2TTfe' 

and  real- valued  effective  permittivity 


~  f  I 

e  —  e 


2./’ 


(13) 


(14) 


each  comprising  the  combined  influence  of  separate  processes  that  are  unresolved 
by  measurement.  The  composite  nature  of  effective  electromagnetic  properties  is 
illustrated  in  Figure  2,  where  measured  high-frequency  electrical  parameters  are  com¬ 
pared  with  effective  properties  predicted  by  Equations  (13)  and  (14).  Theoretically 
predicted  spectra  are  based  on  Debye  models  for  complex  electrical  conductivity  and 
permittivity  as  depicted  in  Figure  1. 

The  propagation  constant  also  incorporates  the  complex  magnetic  permeability 


/  •  // 

fX  =  [X  —  tfX 


(6) 


and  in  connection  with  the  analysis  of  the  ground  radar  measurements,  it  is  commonly 
assumed  that  ii"  =  0,  yielding  an  effective  real- valued  permeability  jl  =  ii' .  While 
there  are  established  and  signihcant  exceptions  to  this  assumption  for  GPR,  it  is 
not  an  unreasonable  simplihcation  at  radar  frequencies.  In  contrast,  however,  it  is 
certainly  not  valid  at  the  lower  operational  frequencies  of  metal  detectors,  particularly 
in  electrically  resistive  soils.  On  the  contrary,  it  is  now  widely  accepted  that  anomalous 
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magnetic  susceptibility  and  related  viscosity  are  in  large  part  responsible  for  observed 
soil-noise  and  associated  degradation  of  metal  detector  performance. 

Finally,  in  the  previous  report  (Cross,  2000),  approximate  forms  of  the  propagation 
constant  for  non-conductive  (cT/27r/e  ^1)  and  conductive  (cj/27r/e  ^1)  media  were 
considered  in  connection  with  representative  soils  property  data.  As  displayed  in 
Figure  3a,  corresponding  governing  relations  from  Equations  (6)  and  (7)  describe 
unattenuated  wave  propagation  for  non-conductive  media  while  for  conductive  media 
the  mode  is  diffusion.  As  suggested  by  Figure  3b,  modelling  of  the  GPR  response  with 
representative  soil  parameters  demonstrated  that  even  at  radar  frequencies  appropri¬ 
ate  for  mine  detection  (~  1.0  GHz),  the  non-conductive  or  dielectric  approximation 
7  ~  i27r/(//e)^/^  is  not  generally  appropriate.  Now,  in  connection  with  metal  de¬ 
tectors,  we  hnd  ourselves  at  the  opposite  end  of  the  spectrum  and  we  consider  the 
low-frequency  or  conductive  approximation  7  ^  (i27r///cT)^/^.  Although  theoretical 
analyses  of  the  electromagnetic  induction  problem  are  nearly  uniform  in  this  so-called 
quasi-static  assumption,  and  while  it  is  certainly  more  viable  than  a  dielectric  as¬ 
sumption  for  GPR,  data  in  Figure  3b  suggest  that  limitations  could  be  signihcant  at 
higher  operational  frequencies  as  well  as  for  time-domain  measurements  However, 
assuming  in  general  that  displacement  currents  are  negligible  at  typical  metal  detector 
frequencies.  Equations  (6)  and  (7)  yield  corresponding  governing  relations 


or,  equivalently. 


E  =  7^  E  =  i2'Kfixa  E 

(16) 

H  =  7^  H  =  i27r///(j  H 

(16) 

=  (17) 


0  an 

=  at- 

(18) 

In  general,  solutions  of  foregoing  diffusion  equations  describe  a  strongly  damped  os¬ 
cillatory  electromagnetic  held  advancing  at  a  diffusion  rate  of  Vd  =  2(7r///iCT)^/^  and 
with  attenuation  rate  a  =  l/dg  =  (tt/z/ct)^/^,  where  dg  is  the  associated  skin  depth. 
This  is,  in  essence,  the  nature  of  the  metal  detector  held  within  the  host  soil.  The 
actual  held  depends  on  the  specihc  geometry,  orientation  and  height  of  the  source  coil, 
time-dependence  of  the  applied  current  waveform  and,  of  course,  soil  electromagnetic 
properties. 


^  Note  Stewart  et  al.  (1994)  observed  that  in  practice,  displacement  currents  become  significant 
as  their  relative  magnitude  exceeds  approximately  5  percent  of  corresponding  conduction  currents 
and  that  this  condition  is  commonly  satisfied  at  frequencies  between  50  and  100  kHz.  Clearly,  the 
transition  frequency  could  be  lower  for  dry  surface  soils. 
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3.0  METAL  DETECTORS  — 

INFLUENCE  OF  SOIL  EM  PROPERTIES 

Metal  detectors  operate  on  the  basis  of  electromagnetic  indnction  (EMI)  as  illnstrated 
in  Fignre  4  In  general,  transmitter  circnitry  impresses  a  time- varying  cnrrent  wave¬ 
form  Iq  npon  the  transmit  coil  to  generate  a  primary  magnetic  held  Hp.  Associated 
magnetic  Hnx  Bp  =  /IHp  passing  throngh  the  receiver  or  search  coil  generates  a  re¬ 
lated  primary  signal  cnrrent  Ip.  The  primary  held  also  penetrates  the  snbsnrface, 
indncing  a  combination  of  magnetization  M.  and  electrical  eddy  cnrrents  Ie  within 
the  host  soil,  proportional  to  bnlk  magnetic  permeability  ni  and  electrical  condnctiv- 
ity  cTi,  respectively.  On  average,  resnlting  secondary  magnetic  helds  associated  with 
time-varying  magnetization  and  eddy  cnrrents  within  the  soil,  combine  to  indnce  a 
secondary  backgronnd  cnrrent  in  the  receive  coil.  Greater  or  lesser  magnetization 
and  eddy  cnrrent  indnction  within  localized  heterogeneities,  having  anomalons  per¬ 
meability /condnctivity,  yield  corresponding  anomalons  secondary  helds  and  a  related 
anomalons  cnrrent  la  in  the  receive  coil.  Conseqnently,  net  cnrrent  in  the  receiver  coil 
is  Is  =  Ip-\-  Ib-\-  Iai  rehecting  the  contribntion  of  primary,  backgronnd  and  anomalons 
magnetic  helds.  As  depicted  in  Fignre  4,  the  anomalons  held  can  be  due  to  the  target, 
where  is  considered  ’’signal”  or  due  to  localized  heterogeneity  of  the  host  soil,  with 
la  considered  ’’noise”.  Of  course,  in  practice  the  anomalous  component  is  generally 
some  combination  of  signal  and  noise. 

There  are  two  principal  classes  of  metal  detectors:  continuous  wave  (CW)  and  pulse 
induction  (PI).  The  essential  difference  is  the  nature  of  the  time- variable  source  cur¬ 
rent  waveform  as  illustrated  in  Figure  5.  The  transmit  coil  of  a  continuous  wave 
or  frequency-domain  (FDEM)  metal  detector  is  driven  by  a  continuous  sinusoidal 
current  waveform.  Through  appropriate  null  adjustment  and  receiver  circuitry,  am¬ 
plitude  and  phase  of  secondary  signal  voltage  Vs  =  vt  +  Va,  due  to  background  and 
anomalous  helds,  are  measured  relative  to  the  primary  signal  voltage  Vp.  Alterna¬ 
tively,  the  response  may  be  quantihed  by  measuring  the  amplitudes  of  secondary 
signal  components  in-phase  and  quadrature  (90  degrees  phase-shifted)  with  respect 
to  the  primary  signal.  In  contrast,  the  pulse  induction  or  time-domain  (TDEM)  metal 
detector  employs  a  transient  source  current  waveform  and  following  termination  of 
current  in  the  transmit  coil  at  time  to,  measures  the  amplitude  of  the  decaying  sec¬ 
ondary  signal  voltage  in  the  receiver  coil.  The  secondary  voltage  waveform  is  sampled 
at  one  or  more  specihed  time-offsets  {tsi,ts2,ts3,  •••)  to  dehne  the  signal  decay  rate. 


2 

For  informative  reviews  of  metal  detector  design  operation  and  application  see  Bruschini  (2002) 
and  Guelle  et  al.  (2004). 
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As  per  standard  treatment  (Grant  and  West,  1965;  Keller  and  Frischknecht,  1966; 
Telford,  et  al.,  1990),  basic  physical  insight  into  metal  detector  response  and  related 
dependence  on  material  electromagnetic  properties  can  be  obtained  by  investigating 
the  response  of  an  eqni valent  Inmped  element  circnit,  inclnding  a  series  combination 
of  resistance  and  indnctance.  The  inflnence  of  material  parameters  and  frequency  are 
investigated  by  quantifying  the  mutual  inductive  coupling  between  the  model  circuit, 
and  the  metal  detector’s  transmitter  and  receiver  coils.  In  particular,  the  model 
comprises  a  shorted  turn,  having  effective  resistance  R  and  effective  inductance  L. 
The  model  is  illustrated  in  Figure  6.  For  now,  we  shall  neglect  any  influence  of  the 
non-conductive,  magnetically  permeable  core. 


3.1  Frequency  Domain  (FDEM)  Response: 

A  time- varying  current 

/  =  /oe*2^^^  (19) 

in  the  transmit  coil  generates  a  primary  time-varying  magnetic  held  Hq  that  couples 
with  both  target  and  receiver  coils.  Associated  time- varying  magnetic  Hux  through 
the  receiver  induces  an  electromotive  force  or  voltage 

vo2  =  —Mo2  —  =  —Mo2—  =  — Mo2i27r/  Jo  (20) 

at  at 

in  the  receiver  coil  where  M02  denotes  the  related  mutual  inductance  ^  between  the 

o 

The  mutual  inductance  between  two  closed  loops  is  a  scaler  constant  defined  by 


Mi2  = 


3*12 

“T’ 


where 


$12  —  /  B2  •  da2  —  IJ-2  I  Hq  •  da.2 
J  A.2  d  A.2 

is  the  magnetic  flux  generated  by  time- varying  current  I  in  loop  1  and  threading  loop  2.  The  vector 
magnetic  flux  density  B2  is,  of  course,  related  to  the  magnetic  field  Hi,  produced  by  time  varying 
current  I  in  loop  1  by  the  constitutive  relation  B2  —  /t2Hi,  where  p.2  denotes  the  magnetic 
permeability  of  the  (isotropic)  medium  surrounding  loop  2  and  dSL2  is  a  local  unit  vector  normal  to 
the  cross-sectional  area  A2  of  loop  2.  The  reciprocity  theorem  yields  M21  =  M\2- 
Note  that  Equation  (2)  is  a  form  of  Faraday’s  induction  law 


V  = 


Eo  •  dsn 


'C2 


dL 

dt 


Bo  •  dsio 


'  a.2 


(i$i2 
dt  ’ 


that  is,  the  electromotive  force  V  around  loop  2,  given  by  the  curl  or  line  integral  of  the  related 
electric  field  vector  E2  around  the  circuit  C2  of  loop  2,  is  equal  to  the  negative  time  rate  of  change 
of  magnetic  flux  B2  threading  the  loop.  In  accordance  with  Lenz’s  law,  the  negative  sign  implies 
that  the  induced  voltage  yields  a  current  in  loop  2  that  opposes  the  source  current  in  loop  1. 
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transmit  coil  (0)  and  the  model  circnit  (1). 

Simnltaneonsly,  primary  magnetic  flnx  throngh  the  target  coil  indnces  a  voltage 

^^01  =  -Moi  ^  =  -Moi  i27Tf  lo  (21) 

where  Mqi  denotes  the  corresponding  mntnal  inductance  between  transmitter  and 
target  coils.  An  associated  current  Ii  is  developed  in  accordance  with  a  generalized 
Ohm’s  law 


=  =  ^ ^  =  -Moi  ,  ,  Jo 

Z  R  +  i27TffXiL  R  +  i2'KfL 


(22) 


where  Z  =  R  +  i27r/L  denotes  the  effective  complex  impedance  of  the  model  circuit. 
Note  that  R  =  Ie  is  the  model  eddy  current.  The  eddy  current,  in  turn,  produces  a 
secondary  magnetic  held  at  the  receiver  coil,  inducing  a  related  voltage 

vi2  =  =  -M01M12  Iq  (23) 

dt  R  +  i27r/L 

where  M12  is  the  mutual  inductance  between  the  model  circuit  and  the  receiver  coil. 

Finally,  assume  that  by  appropriate  electronics  and  nulling,  the  metal  detector  mea¬ 
sures  the  voltage  ratio 


^  _  -M01M12  (27r/)"/oe^^"^^ 

V02  —Mo2i27rf{R  +  i27rfL)lQe'^‘^'^R 

M01M12  i{27TfL/R){l-i27TfL/R) 

M02L  1  +  {27TfL/Ry 

_  M01M12  {2'nfl/Rf  +  i{27TfL/R) 

M02L  1  +  {27TfL/Ry 


which  is  commonly  expressed  as 


V12 

V02 


M01M12  /'a‘^  +  ia\ 

M02L  V 1 + oO 


(24) 


(25) 


where  K^y  =  M01M12/M02T  is  the  so-called  coupling  coefficient  and 


Fa{a) 


ia\  a‘^  ia 

1  +  J  1  +  1  + 


(26) 


is  the  associated  response  function  in  terms  of  response  parameter  a  =  2'Kf  L/R. 


Note  that  while  the  coupling  coefficient  is  a  purely  geometrical  relation  dependent 
on  transmitter  receiver  configuration  and  model  geometry,  the  complex  inductive 
response  function  is  obviously  dependent  upon  frequency  and  the  electromagnetic 
properties  of  the  model  circuit  or  equivalent  soil  as  characterized  by  the  response 
parameter  (effective  induction  number)  a.  Although  the  negative  polarity  indicated 
in  Equations  (24)  and  (25)  is  more  commonly  incorporated  in  the  coupling  coefficient, 
it  is  instructive  and  perhaps  more  appropriate  to  retain  it  as  an  integral  part  of  the 
response  function. 

Real  and  imaginary  parts  of  the  resulting  induction  response  function  are  depicted  in 
Figure  7a,  together  with  corresponding  amplitude  and  phase  spectra 

(27) 

Note  that  for  a  — )■  0,  the  so-called  resistive  limit,  the  negative-valued  quadrature 
component  dominates  as  the  overall  response  diminishes  toward  zero  for  either  a 
static  source  held  (/  =  0)  or  a  resistive  (i?  =  cx)  or  cr  =  0)  ground.  Equivalently, 
as  the  response  amplitude  trends  toward  zero,  the  secondary  signal  voltage  (^12)  is 
phase-shifted  or  lagged  by  90  degrees  relative  to  the  primary  signal  voltage  (no2)-  In 
contrast,  as  frequency  increases  and/or  the  anomalous  ground  becomes  more  conduc¬ 
tive,  both  in-phase  (real)  and  quadrature  (imaginary)  components  increase,  initially 
the  quadrature  more  rapidly  and,  subsequently,  the  in-phase  until  the  two  contribu¬ 
tions  are  equivalent  (Re{Fcr}  =Im{Fcr}  =  —0.5)  for  a  =  1.  Notice  that  this  condition 
is  associated  with  a  relative  phase-shift  of  45  (-135)  degrees  between  primary  and 
secondary  coil  voltages,  or  what  is  effectively  the  transition  between  resistive  and 
inductive  response  regimes. 

With  further  increase  in  frequency  or  conductivity  a  — )■  cx),  the  response  function 
becomes  purely  real-valued,  reaching  an  absolute  maximum  {\Ffj\  =  1.0)  referred  to 
as  the  inductive  limit.  The  secondary  signal  voltage  (^12)  becomes  saturated  and  180 
degrees  out  of  phase  (Re{Fcr}  =  —1.0)  with  the  primary  (^02)-  Practically  speaking, 
in  the  inductive  limit,  eddy  currents  are  restricted  to  the  surface  of  the  conductive 
medium  (soil)  and  the  associated  secondary  held  effectively  cancels  the  primary  held 
within  the  conductor.  In  contrast,  for  a  resistive  ground  and/or  at  lower  frequencies, 
the  eddy  currents  penetrate  more  deeply. 

Now,  beyond  the  standard  treatment  of  an  anomalously  conductive  ground,  we  ex¬ 
tend  the  model  to  consider  the  added  influence  of  anomalous  magnetic  permeability. 
Again,  neglecting  any  interaction  between  the  shorted  turn  and  permeable  core,  the 
core  is  magnetized  parallel  to  the  primary  magnetic  held,  yielding  a  net  magnetic 
moment  that  reinforces  the  primary  held  at  the  receiver  coil  In  particular,  assum¬ 
ing  real- valued  (lossless)  permeability,  magnetization  of  the  permeable  core  produces 

^  Note  that  within  a  magnetically  polarizable  medium  (e.g.  a  magnetically  susceptible  soil),  the 
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a  secondary  magnetic  field  that  is  pnrely  in-phase  with  the  primary  held  and  in¬ 
dependent  of  freqnency.  Thns,  adding  a  corresponding  real-valned  magnetization 
component  —  1)  to  the  indnctive  response  in  Eqnation  (9),  we  obtain 

the  net  non-interacting  response  function 


Fa+^JL  —  Ffj  -\-  F^  — 


F  1) 


a 


1  +  a'^ 


+ 


t  a 


1  -|-  ’ 


(28) 


where  is  a  geometry-dependent  factor  that  appropriately  scales  the  inductive 
coupling  coefficient  Kcr  and  Hr  =  m/mo  denotes  the  relative  magnetic  permeability 
Corresponding  response  parameters  are  displayed  in  Figure  7b,  assuming  an  arbitrary 
relative  permeability  such  that  =  0.333Re  {F'cr(cx))}. 

In  contrast  with  the  purely  inductive  case  in  Figure  7a,  the  net  low-frequency  or 
resistive  {a  — )■  0)  response  |Fcr+^|  is  now  dominated  by  the  magnetization  response, 
with  zero  phase-shift  between  (^12)  and  (no2)-  Again,  with  increasing  frequency  or 
conductivity,  the  absolute  response  amplitude  increases,  albeit  more  gradually,  and 
there  is  a  polarity  reversal  of  the  in-phase  component  consistent  with  transition  from 
magnetization  to  induction  dominated  response.  Equivalently,  the  phase-shift  be¬ 
tween  secondary  and  primary  signal  voltages  transitions  from  zero  to  -180  degrees, 
approaching  the  inductive  limit  for  a  — )■  cx). 

vector  magnetization  IVI  is  defined  as  the  difference  between  the  observed  magnetic  flux  density 
B  =  //H  and  the  corresponding  free-space  magnetic  flux  density  Bq  —  MqH,  normalized  by 
free-space  permeability  Hq.  Thus  we  have 

^  ^  B  -  /tpH  ^  HrlJ-pH  -  MqH 
Mo  Mo 

=  {Hr  -  1)H 

=  XvH 


where  flj.  =  m/ Mo  is  the  relative  permeability  and  Xv  —  (Mr  —  1)  is  the  related  volume  magnetic 
susceptibility.  Now  for  a  volume  limited  distribution  of  magnetized  material,  the  net  vector  magnetic 
moment  A4  is  just  the  integrated  magnetization 


M=  /  Mdv. 

Jv 

Conversely,  the  magnetization  is  simply  the  magnetic  moment  per  unit  volume.  Of  course,  at  suf¬ 
ficient  distance  from  the  magnetized  source  material,  the  resulting  field  is  indistinguishable  from 
that  of  a  infinitesimal  magnetic  dipole  with  equivalent  moment  and  direction,  or  for  that  matter  an 
equivalent  current  loop  like  that  of  our  model. 

Note  that  a  constant  magnetization  response  r  ^  implies  that  related  secondary  signal  voltage 
generated  in  the  receiver  coil  is  directly  proportional  to  frequency  as  for  the  primary  signal. 
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The  non-interacting  model,  however,  predicts  this  indnction  limited  response  remains 
diminished  by  the  offsetting  inffnence  of  freqnency-independent  magnetization.  In 
fact,  however,  on  acconnting  for  interaction  between  core  and  coil,  the  indnctive  re¬ 
sponse  is  also  enhanced  by  magnetization  within  the  core  and  vice  versa.  Related 
effects  are  modelled  in  Fignre  8a,  where  polarity  of  the  response  fnnction  is  reversed 
to  be  consistent  with  common  convention.  Note  that  while  transition  to  the  indnctive 
limit  is  snbstantially  more  gradnal,  the  limiting  high-freqnency  response  {a  — )■  cx)) 
for  an  anomalous  conductive-magnetic  ground  remains  purely  inductive,  consistent 
with  a  -180  degree  phase-shift.  In  effect,  the  secondary  magnetic  held  associated 
with  saturation  eddy  currents  virtually  cancels  the  primary  held  within  the  medium 
(permeable  core)  and,  consequently,  nullihes  related  inhuence  of  magnetization.  In 
practice,  magnetization  associated  with  increased  permeability  enhances  both  pri¬ 
mary  and  induced  helds  with  equivalent  net  effect.  Moreover,  because  effective  model 
inductance  L  increases  with  core  permeability,  the  entire  response  is  correspondingly 
displaced  as  a  function  of  induction  number  a.  In  particular,  transition  from  resistive 
to  inductive  response  regimes  occurs  for  an  induction  number  a  >  1,  increasing  with 
increasing  permeability  (susceptibility). 

Comparison  with  Figure  8b  indicates  that  the  modelled  response  is  consistent  with 
that  expected  for  a  horizontal  coplanar  coil  system  over  a  homogeneous  conductive- 
permeable  soil  (Huang  and  Won,  2000)  In  particular,  response  curves  are  displayed 
in  Figure  8b  for  an  effective  induction  number  a^f.  Depicted  response  characteristics 
assume  a  coil  separation  s  =1.66  m  and  are  similar  for  more  closely  spaced  or  coaxial 
coils.  In  general,  as  coil  separation  decreases,  peak  quadrature  response  and  related 
in-phase  inhection  occur  for  a  lower  effective  induction  number.  Corresponding  mag¬ 
nitude  and  phase  spectra  are  depicted  in  Figure  8.5. 

In  effect,  the  response  parameter  a  (modelled  as  a  =  27t f  L/R)  is  equivalent  to  the 
attenuation  rate  as  dehned  in  Chapter  2.0.  Although  we  have  referred  to  a  rather 
loosely  as  the  “induction  number” ,  the  actual  induction  number  is  more  generally 

r]  =  as=  -^  =  (29) 

where,  dg  =  l/a  =  (7r///(j)“^/^  is  the  corresponding  skin  depth  as  introduced  in 
Section  2.0  and  s  represents  an  appropriate  dimensional  parameter.  For  the  case  of  a 
bistatic  dual-coil  system  considered  in  Figure  8/8.5,  s  denotes  the  separation  between 
coplanar  coil  axes.  Assuming  s  =  1.0  to  facilitate  scaling.  Figure  9  provides  some 
indication  of  the  connection  between  soil  electromagnetic  parameters  and  induction 
number  as  related  to  associated  FDEM  response  characteristics  modelled  in  Figure 

®  It  is  emphasized  that  the  range  of  modelled  permeabilities  in  Figures  7  and  8/8.5  greatly  exceeds 
that  of  natural  soils  and  other  geological  materials.  Model  permeability  is  exaggerated  for  purpose 
of  illustration. 
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8.  In  particular,  Figure  9  charts  electrical  conductivity  cr  as  a  function  of  magnetic 
susceptibility  x  =  (h/^o  —  1)  for  77  =  a  =1.0,  0.1  and  0.01  and  with  operating 
frequency  (/  <  100  kHz)  as  a  parameter. 

First,  it  is  generally  observed  that  the  induction  number  is  comparatively  insensitive 
to  magnetic  susceptibility  within  the  nominal  range  of  electromagnetic  parameters 
for  natural  soils.  In  fact,  this  is  to  be  expected  given  the  equal  weight  of  magnetic 
permeability  and  electrical  conductivity  and  the  relatively  minimal  deviation  of  per¬ 
meability  //  =  //o(l  +  x)  on  fhe  scale  of  natural  soil  susceptibility  variation.  However, 
this  is  not  to  say  that  instrument  response  is  insensitive  to  magnetic  susceptibility 
fluctuation.  In  fact,  as  suggested  by  Figure  8  and  for  soils  possessing  a  signihcant 
magnetic  susceptibility,  related  magnetization  yields  the  bulk  of  measured  signal  for 
resistive  soils  and  at  lower  frequencies,  where  a  <<  1.  Although  the  induction  num¬ 
ber  depends  primarily  on  the  frequency-conductivity  fa  product,  it  is  clear  that  for 
typical  operating  frequencies  the  vast  majority  of  soils  satisfy  a  <  1.  In  fact,  the 
chart  for  a  =  1  indicates  that  the  entire  range  of  nominal  soil  electromagnetic  prop¬ 
erties  satishes  the  condition  a  <  1.0  for  frequencies  less  than  /=100  kHz.  Similarly, 
while  a =0.1  and  0.01  are  clearly  more  limiting  conditions,  even  soils  with  signihcant 
electrical  conductivity  approaching  1.0  mS/m  satisfy  the  condition  a  <  0.01  for  an 
operating  frequency  of  /=30  kHz. 

Of  course,  the  actual  induction  number  77,  related  response  characteristics  and  as¬ 
sociated  soil  inhuence  for  a  given  FDEM  device  are  dependent  on  the  specihc  coil 
conhguration  and  operating  frequency.  Consequently,  beyond  foregoing  general  de¬ 
velopment,  it  is  useful  to  consider  specihc  theory-predicted  results  for  a  horizontal 
coaxial  coil  arrangement  as  representative  of  common  metal  detectors. 

In  general,  for  source  current  /  impressed  on  a  horizontal,  single-turn  circular  coil  of 
radius  o  positioned  height  ^  =  —h  above  a  homogeneous,  linear  and  isotropic  soil, 
the  resulting  electromagnetic  held  comprises  the  following  orthogonal  electric  and 
magnetic  held  components  as  developed  in  Section  4  (Equations  4.86-4.88)  of  Ward 
and  Hohmann  (1987) 


=  -m  filial 


^-uo{z+h) 


—  Ji  (Ao)  Ji  ( Ap)  d\ , 

Uq 


(30) 


and 


Hp  = 


al 


^-uo{z+h) 


+  TteC 


uo{z-h) 


A  Ji(Ao)  Ji(Ap)  dX 


^-uo{z+h)  ^^^^uo{z-h) 


— Ji(Ao)  Jo{Xp)  dX  . 

Uo 


(31) 

(32) 


Here,  E^f),  Hp  and  Hz  denote  the  azimuthal  electric  held,  the  radial  magnetic  held 
and  the  vertical  magnetic  held,  respectively.  Jo()  and  Ji()  denote  Bessel  functions 
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of  first  kind  and  integer  orders  0  and  1,  respectively,  and  A  is  the  Hankel  transform 
variable.  The  consistent  portion  of  the  integrand  in  brackets  comprises  an  initial  term 
associated  with  the  primary  held  and  a  second  term  dne  to  the  soil-related  held,  where 


HiUq  -  HoUi 
HiUq  -I-  HoUi 

denotes  the  ehective  rehection  coefhcient,  —  7I,  and 

1/2 

[-127: fukCJkY^'^ 


Ik  =  ‘27:  f 


t-lk  ^k 


O'k 


27r/efc 


-  1 


(33) 


(34) 


is  the  qnasi-static  (see  Section  2.0)  propagation  constant  for  air  {k  - 
{k  =  1),  respectively.^  For  air,  we  have  (Tq  =  0,  //q  =  47r  x  10“^ 
uq  =  X,  yielding 

//I  A  -  fiQUi 
HiX  +  Houi 


=  0)  and  for  soil 
H/m  and,  thns 

(35) 


For  the  specihc  case  of  a  coincident  horizontal  (coaxial)  coil  system,  Eqnation  (32) 
yields  the  following  relation  for  the  ratio  of  secondary  Hg  to  primary  Hp  magnetic 
helds  threading  the  receiver  coil  (Won  and  Hnang,  2004).® 

CXD 

rxB  Ji(Ao)dA.  (37) 

Associated  in-phase  and  qnadratnre  response  characteristics  are  displayed  as  fnnctions 
of  indnction  nnmber  r]  =  (7r///cT)^/^o  in  Fignre  10a.  Note  that  in  this  case,  the 
dimensional  parameter  s  =  o  is  the  coil  radins. 

Eqnivalent  results  are  described  by  Bruschini  (2004)  for  the  corresponding  secondary 
signal  voltage 


v{f)  =  -i27:fiiQ7:a^I 


tte  e  [Ji(Ao)]^  dA  , 


(38) 


that  follows  from  Equation  30  or  32  with  .2  =  —h  and  p  =  o.^  Related  (absolute)  re¬ 
sponse  curves  are  displayed  in  the  right-hand  panel  of  Figure  10  for  both  non-magnetic 


Note  the  corresponding  induction  number  p  is  the  product  of  the  real  part  of  the  quasi-static 
propagation  constant  Re\p)'\  =  OL  =  \ / dg  =  (tt/z/Ct)^/^  and  the  dimensional  parameter  S. 

O 

The  corresponding  field  ratio  for  horizontal  coplanar  coils  (Figure  8b)  is 
Hg/Hp  =  s®  rxE  A^  Jo(Ao)  dX  (37.5). 

^  See  Section  3.3.1  for  details.  In  general,  integral  expressions  (37)  and  (38)  are  evaluated  nu¬ 
merically.  See  for  example  Anderson  (1979,  1982);  Johansen  and  Sorensen  (1979). 
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and  magnetic  soils  with  finite  electrical  conductivity.  Note  that  the  corresponding  re¬ 
sponse  parameter  a  =  27Tfnaa^  is  related  to  the  induction  number  r]  by  the  relation 
a  =  2rf. 

In  general,  complex  response  characteristics  in  Figure  10  are  consistent  with  those 
predicted  by  the  heuristic  model  (Figures  7-8)  and  illustrate  the  very  significant  in¬ 
fluence  of  soil  magnetic  susceptibility.  More  specifically,  results  indicate  that  the 
resistive  limit  is  extended  to  higher  induction  numbers  with  increasing  magnetic  sus¬ 
ceptibility  and  that,  within  the  resistive  limit  the  in-phase  component  of  the 
response  is  predominant  and  proportional  to  magnetic  susceptibility. 

In  particular,  for  single-turn  coincident  coils  located  height  h  above  a  non-conducting 
and  non-dispersive  magnetic  soil.  Equation  38  yields  (Das,  2004;  2006)  the  corre¬ 
sponding  secondary  signal  voltage 


Vi2  =  —i2TTflJ,QTTa‘^I 


X 


L2  +  X 


m{h)  ~  —iTTfuoTTO^IXi^ih), 


(39) 


where  x  =  (//r  — 1)  denotes  the  volume  magnetic  susceptibility  and  m{h)  is  a  stand-off 
dependent  coupling  factor  analogous  to  in  foregoing  development  (see  details  in 
Section  3.3.3).  For  natural  soils  (x  <<  2),  we  obtain  the  indicated  approximation 
with  Vi2  proportional  to  x  and  in  phase  with  the  corresponding  primary  signal  voltage 
vo2  =  -f27r///o7ro^/  [Ji(Ao)]^  dA. 

As  regards  practical  constraints  on  the  resistive  limit,  the  lower  left  diagram  in  Figure 
10  (Won  and  Huang,  2004)  displays  the  accuracy  of  a  resistive-limit  assumption  as  a 
function  of  soil  magnetic  susceptibility  x  and  frequency-conductivity  product  fa.  It 
is  clear  that  even  for  modest  magnetic  susceptibility  x  =  1-0  x  10“^  SI,  related  error 
is  only  marginally  in  excess  of  30%  for  fa=1000  (for  example,  ct=100  mS/m  for  /=10 
kHz). 

Response  characteristics  are  also  substantially  dependent  on  sensor  height  or  stand-off 
above  grade  as  illustrated  in  Figure  11  from  Bruschini  (2004).  In  particular,  it  is  ob¬ 
served  that  the  resistive  limit  for  given  magnetic  permeability  is  effectively  shifted  to 
a  significantly  lower  response  parameter  value  as  normalized  sensor  height  Hn  =  h/a 
increases.  In  addition,  it  is  noted  that  (within  the  resistive  limit)  the  absolute  value 
of  the  in-phase  response  decreases  by  more  than  two  orders  of  magnitude  on  increas¬ 
ing  the  sensor  roughly  from  surface  {Hn  =  0.01  as  per  Figure  10)  to  height  above 
grade  equal  to  coil  radius  {Hn  =  1-0).  In  comparison,  corresponding  reduction  in 
quadrature-phase  response  is  relatively  limited.  Relative  influence  reverses,  however, 
as  the  response  parameter  approaches  the  inductive  limit. 

As  a  common  rule  of  thumb,  the  resistive  limit  for  a  non-magnetic  ground  is  ?7  <  0.02  (e.g. 
Frischknecht,  1987). 
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A  final  consideration  for  freqnency  domain  devices  is  the  inflnence  of  dispersive  soil 
electromagnetic  properties.  In  theory,  associated  effect  conld  be  investigated  by  ap¬ 
propriate  generalization  of  the  foregoing  henristic  model  (i.e.  substitnting  a  com¬ 
plex  coil  impedance  R  —  i/27TfC  for  R  and  allowing  a  complex  core  permeability 
Hr  =  Hr  ~  However,  to  hrst  order,  related  qnadratnre  inflnence  is  predictable 

in  a  qnalitative  sense.  Compared  with  an  equivalent  non-dispersive  soil,  electrical 
polarization  produces  a  relative  increase  in  the  in-phase  response  (for  common  fa) 
while  magnetic  viscosity  yields  an  effective  increase  in  the  quadrature  response. 

For  example,  generalizing  the  previous  result  in  Equation  (39)  for  a  non-conductive 
soil  with  frequency  dependent  complex  magnetic  susceptibility  xif)  =  x'  Equa¬ 
tion  (38)  yields  a  complex- valued  secondary  voltage  Vi2{t)  =  Vi2  with  in-phase 

and  quadrature  signal  components 

^i2(/)  = -^27r///o7ro^d  ^  ^  ^,2  rn{h)  ^ -tt  f  hqt^o^ I  x"  m{h)  (40) 

and 

r-y'2  ,  2y'  + 

^12(/)  = -^27r///o7ro^/  ^  2)2  +  ^ -mf  HoT^a^I  x' m{h),  (41) 

respectively  (see  Das  2004;  2006).  Note,  as  indicated,  the  real  part  of  the  secondary 
signal  voltage  v'12  is  approximately  proportional  to  the  quadrature  susceptibility  x" 
and  lags  the  primary  signal  voltage  vq2  by  90  degrees. 

In  general,  quadrature  parameters  are  small  compared  with  corresponding  real  parts 
and,  consequently,  soil  electromagnetic  dispersion  is  anticipated  to  have  relatively  lim¬ 
ited  influence  on  frequency  domain  (FDEM)  devices  at  a  given  operating  frequency. 
For  multi- frequency  sensors  and/or  in  connection  with  compensation  adjustments, 
however,  influence  can  be  signihcant.  More  signihcant,  however,  is  the  pivotal  in¬ 
fluence  of  soil  electromagnetic  dispersion  and  related  viscosity  on  the  response  and 
performance  of  time  domain  metal  detectors.  Related  response  characteristics  and 
implications  are  discussed  in  the  following  Section  and  in  greater  detail  in  Section 
3.3. 

3.2  Time  Domain  (TDEM)  Response: 

As  demonstrated  by  Grant  and  West  (1965),  the  same  heuristic  model  employed  to 
investigate  the  FDEM  response  can  be  extended  to  provide  useful  insight  on  the  nature 
of  the  time-domain  response  and  related  influence  of  soil  electromagnetic  properties. 
Neglecting  the  charging  portion  of  the  applied  transmit  current  waveform,  we  assume 
a  step  termination  of  constant  current  and  take  the  termination  as  zero  time  reference. 
Thus,  we  write  the  drive  current  waveform  as 

/  =  Jo  [1 -'u(t)],  (42) 
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where 


r  0;  t  <0 

u{t)  =  <  t/At]  0>t<At  (43) 

I  1;  t>At 

defines  a  step  function  u{t)  with  hnite  rise-time  At.  Again,  we  begin  by  ignoring  the 
influence  of  the  permeable  core.  As  per  Equation  (21),  the  induced  voltage  in  the 
model  circuit  is  given  by 


Voi  —  —Mqi  —  —  —Mqi  ~  —  Mqi  Iq 


(44) 


Now,  as  per  Kirchoff’s  circuit  conventions,  we  must  have  that  the  induced  emf  is 
balanced  by  related  potential  drops  across  effective  resistance  R  and  inductance  L, 
yielding 


Mqi  Iq 


du{t) 

dt 


(45) 


where  R  =  Ie  denotes  the  model  eddy  current  induced  in  the  shorted  turn  (compare 
with  the  equivalent  frequency  domain  expression  in  Equation  (22).  Of  course,  prior  to 
constant  current  termination  dl/dt  =  Uoi  =  /i  =  0.  Subsequently,  the  source  current 
is  terminated  and  falls  to  zero  as  IqR/ At).  Thus,  during  the  finite  termination  interval 
At,  we  have  du{t)/dt  =  1/At  and  the  previous  equation  yields 


Mqi  Iq 

At 


Rh  +  l"^ 

dt 


or,  equivalently  the  differential  equation 


(wf + 


Mqi  Iq 

RAt  ’ 


(46) 


(47) 


yielding  solution 


Mqi  Iq 

RAt 


1  _  e-(-R/b)i 


(48) 


under  condition  /i  =  0  for  t  =  0.  However,  given  t  «  L / R  ior  t  <  At  the  approx¬ 
imation  ^  1  +  X  leads  to  the  following  linear  relation  for  initial  “eddy  current” 
build-up  in  the  model  circuit 

_  Mqi  Iq  t 

Subsequently,  for  t  >  At,  when  there  is  no  longer  current  flowing  in  the  transmit  coil, 
we  are  left  with  the  homogeneous  (no  source  term)  equivalent  of  Equation  (47) 


(£/fl)f +  /:=0. 


(50) 
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with  purely  exponential  solution 


L 


(81) 


where  initial  amplitude  is  determined  by  the  condition  Ii  =  MqiIq/L,  for  t  =  At 
from  Equation  (49).  Thus,  we  hnd  that  upon  terminating  a  constant  source  current 
in  the  transmit  coil,  there  is  initially  a  linear  build-up  of  the  “eddy  current”  within 
the  model  circuit  that  parallels  the  linear  decay  of  the  source  current.  Subsequently, 
this  eddy  current  decays  exponentially  with  time  constant  r  =  L/R  depending  upon 
the  equivalent  lumped  circuit  parameters  of  the  model.  Finally,  of  course,  the  induced 
time-varying  current  within  the  model  circuit  generates  a  secondary  magnetic  held 
that,  in  turn,  produces  a  secondary  signal  voltage  (from  Equations  49  and  51) 


dll  M01M12 

V12  —  —Mi2—r-  — - ~ - Iq 

dt  L 


S(t)  - 


(52) 


where 


d{t) 


0;  t  <  0 

1/At;  0>t<At 
0;  t>  At 


(53) 


represents  an  impulse  having  amplitude  1/At  and  duration  At.  The  corresponding 
primary  signal  voltage  is 


V02  —  —M02—I0  [1  —  u{t)]  —  M02I0  5{t). 


(54) 


The  predicted  time-dependent  response  is  illustrated  in  Figure  12a.  It  is  useful  to 
consider  the  connection  between  these  results  and  corresponding  frequency  domain 
response  spectra  displayed  in  Figure  7a.  In  particular,  it  is  noted  that  the  time-domain 
decay  constant  r  =  L/R  is  closely  related  to  the  transition  between  resistive  and 
inductive  response  regimes.  Specihcally,  we  have  r  =  l/27r/  for  a  =  1.  Comparing 
the  secondary  time-domain  signal  voltage  Vi2{t)  in  Equation  (52)  with  the  equivalent 
frequency  domain  signal  Vi2{f)  given  by  Equation  (23),  it  is  emphasized  that  the 
two  signals  do  not  constitute  a  time-frequency  pair.  Rather,  Vi2{f)  is  the  Fourier 
transform  of  the  impulse  response  of  the  soil  model,  while  Vi2{t)  is  the  corresponding 
step  response  (neglecting  hnite  rise  time  At).  Consequently,  ignoring  polarity,  Vi2{f) 
is  equivalently  the  Fourier  transform  of  the  time  derivative  of  Vi2{t).  That  is 

d 

vi2{f)  =  J  —[vi2{t)]e~^‘^^^^dt,  (55) 

or,  via  the  derivative  property  of  the  Fourier  transform, 

/cx:) 

vi2{t)e~''‘^^^^dt 

-00  5  (56) 

=  i27TfR{vi2{t)} 


17 


where  J^{g{t)}  denotes  the  forward  Fourier  transform  of  an  arbitrary  time-dependent 
functional  argument  g{t).  Inversely,  it  follows  that  the  time-domain  step  response  is 
obtained  via  the  inverse  transform 


where  T  ^{g{f)}  denotes  the  inverse  Fourier  transform  of  an  arbitrary  frequency- 
dependent  functional  argument  g{f). 

Now,  as  for  previous  analysis  of  the  frequency  domain  response,  let’s  consider  the 
qualitative  influence  of  magnetization  within  the  permeable  core.  Again,  we  ignore 
interaction  between  the  core  and  the  shorted  turn.  Prior  to  terminating  the  source 
current  (t  <  0),  there  is  constant  magnetization  M  =  of  the  core  proportional 

to  the  primary  held  Hq.  Moreover,  it  follows  that  related  magnetic  hux  threading 
the  receiver  coil  is  also  in  static  state  and,  consequently,  produces  no  related  signal 
voltage  in  the  receiver  coil.  However,  on  termination  of  source  current  at  t  =  0, 
this  Hux  decreases  as  inductive  Hux  associated  with  increasing  “eddy  current”  in  the 
model  coil  increases  for  0  <  t  <  At.  In  general,  these  Huxes  are  180  degrees  out  of 
phase,  so  that  the  time  rate  of  change  of  magnetization  Hux  supports  that  due  to 
increasing  inductive  Hux.  The  result,  as  suggested  in  Figure  12b,  is  that  the  initial 
impulsive  signal  voltage  induced  in  the  receiver  coil  during  Hnite  termination  time  At 
is  enhanced  by  declining  magnetization  in  the  permeable  core.  Thus,  for  0  >  t  <  At, 
the  total  signal  voltage  is 


Vl2 


MoiM: 
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Iq  +  —  1) 


(58) 


where  denotes  an  appropriate  coupling  coefHcient  and  5{t)  is  given  by  Equation 
(53).  Assuming  real- valued  (nondispersive)  permeability,  however,  magnetization  is 
zero  for  t  >  At  with  no  further  inHuence  on  the  time  domain  response  that  is  subse¬ 
quently  due  only  to  decaying  “eddy-current”  within  the  model  circuit. 

Although  the  foregoing  conclusion  is  consistent  with  the  non-interacting  frequency- 
domain  response  depicted  in  Figure  7b,  it  is  evident  that  on  allowing  interaction 
between  coil  and  core,  magnetization  within  the  core  will  continue  to  exert  a  second 
order  inHuence  as  long  as  decaying  eddy  current  yields  a  Hnite  magnetic  Held  within 
the  core.  Consequetly,  as  per  prior  frequency- domain  analysis,  related  inHuence  is 
largely  restricted  to  the  induction-dominated  response  regime  for  a  >  1.0.  In  view  of 
relations  described  by  Equations  (55-57),  the  relatively  gradual  decay  of  eddy  currents 
and  limited  deviation  from  Hr  =  1.0,  however,  foregoing  analysis  suggests  that  second- 
order  eHects  within  a  conductive-permeable  soil  are  relatively  minor.  In  practice,  for 
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common  soils,  the  result  is  equivalent  to  a  marginal  increase  in  effective  conductivity 
(see  Section  3.3.5). 

Finally,  it  is  emphasized  that  foregoing  analysis  of  the  TDEM  response  has  been 
restricted  to  constant,  real-valued  soil  electromagnetic  parameters.  However,  as  dis¬ 
cussed  in  the  previous  Section,  soil  properties  are  generally  frequency-dependent,  dis¬ 
persive  quantities.  In  particular,  for  soils  with  signihcant  magnetic  viscosity  (quadra¬ 
ture  permeability),  magnetization  does  not  fall  to  zero  in  termination  time  At  as 
indicated  in  Figure  12b  but,  rather,  decays  more  gradually  as  for  the  eddy  current. 
Consequently,  soils  having  significant  quadrature  permeability  can  have  a  substantial 
and,  in  fact,  dominant  influence  on  the  decay  rate  of  the  time-domain  response. 

Recent  research  and  investigations  (Guelle,  2002;  Pasion,  et  ah,  2002,  Billings,  et  ah, 
2003a,  2003b;  Borry,  et  al,  2003;  Das  2004,  2006)  have  conhrmed  the  significance  of 
soil  magnetic  viscosity  and  identihed  related  influence  as  the  primary  factor  restrict¬ 
ing  the  performance  of  time-domain  PI  metal  detectors.  For  example.  Das  (2006) 
demonstrated  that  even  for  an  extreme  soil  electrical  conductivity  of  ex  =  5.0  S/m, 
a  viscous  magnetic  soil  with  only  moderate  static  susceptibility  x  =  0.0005  SI  yields 
a  substantially  greater  response  at  times  exceeding  10  /rs  as  illustrated  in  Figure  13. 
Thus,  in  light  of  growing  consensus  on  the  singular  significance  of  magnetic  viscosity, 
the  following  section  is  devoted  specifically  to  characterizing  the  nature  of  the  viscous 
magnetic  response  and,  concurrently,  to  investigating  the  analogous  and  potentially 
significant  influence  of  electrical  conductivity  dispersion  on  TDEM  sensors. 


3.3  Dispersive  Soil  EM  Properties  and  TDEM  Response: 


3.3.1  General  Formulation: 

We  begin  by  returning  to  Equations  (30-32)  for  the  electromagnetic  held  due  to  source 
current  /  impressed  on  a  horizontal  single-turn  circular  coil  of  radius  o  at  height 
^  =  —h  above  a  homogeneous,  linear  and  isotropic  soil.  The  frequency-dependent 
voltage  induced  in  a  coaxial,  coplanar  single-turn  receiver  coil  of  radius  b  follows  as 
the  integral  of  the  azimuthal  electric  held  (Equation  30)  around  the  circular  coil 

/p2,7r 

E-d£  =  bJ  E^d(t)  (59) 

or,  equivalently  via  Faraday’s  law,  as  the  time  rate  of  change  of  integrated  magnetic 
flux  density  (Equation  32)  through  the  coil 


v{f)  =  -  J  ^-ds  =  -i27r///o 


2n 


HzP  def)  dp. 


(60) 


0  ^0 
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The  result  is 


^(/) 


—i27TfnoTT  abl 


^-uo{z+h)  ^^^^uo{z-h) 


Ji(Ao)  Ji(A6)  dX 


(61) 


and,  on  setting  h  =  a,  z  =  h  =  0,  we  have  for  coincident,  surface  deployed  coils 


where,  again. 


Tte  = 


pOO 

/  [1 -|- tte]  [Ji(Ao)]^dA, 

(62) 

Jo 

piX-  poui 

PiX  +  PoUi 

(38) 

denotes  the  effective  reflection  coefficient  at  the  boundary  between  air  and  soil. 


Finally,  the  corresponding  time-domain  step  response  is  obtained  by  dividing  the 
previous  equation  by  i27r/  (equivalent  to  integrating  the  associated  impulse  response) 
and  evaluating  the  Fourier  transform  with  respect  to  frequency 


v{t)  =  —HoTTO^I 


[1 -I- tte]  d/ [Ji(Ao)]^  dA 


(63) 


or,  equivalently,  by  Laplace  transformation  with  substitution  s  =  i2TTf. 

In  what  follows,  we  review  a  range  of  solutions  for  specified  soil  property  models, 
and  provide  a  comparative  assessment  of  related  influence  on  the  TDEM  response 
for  a  small-scale,  horizontal  coincident- coil  configuration  at  the  surface  of  a  uniform 
half-space. 


3.3.2  Non-Magnetic,  Non-Dispersive  Conductive  Soil: 

For  reference,  we  begin  with  a  uniform  non-magnetic,  non-dispersive  conductive  soil. 
Assuming  ai  =  a^c  and  pi  =  po,  we  have 

A  —  ui  2A  ,  , 

1  -I-  rxE  =  1-1-  - - =  - -  (64) 

and,  consequently,  from  Equation  (63) 

pOO  poo  r\  \ 

v{t)  =  —poTTO^I  /  /  - -  d/ [Ji(Ao)]^  dA  (65) 

Jo  J-oo  X  +  Ui 


with  ui  =  [A^  -|-  127: f poCTdc]-  The  problem  was  initially  treated  by  Lee  and  Lewis 
(1973)  and,  subsequently,  by  Raiche  and  Spies  (1981).  On  evaluating  the  Laplace 
transform  of  1/(A  -|-  Ui)  and  redefining  the  integration  variable  A'  =  oA,  the  resulting 
solution  is 


v{t)  =  Vi{t) 


—2pQ^aI 

t 


S{t), 


(66) 
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where 


S{t) 


a/^  A' erfc(-\/T  A')  X' [Ji{X')]^  dX' 


(67) 


is  an  integral  function  of  normalized  time  r  =  t/afXoa^,  with  erfc()  denoting  the 
complementary  error  fnnction.  With  snbstitntion  r  =  tA'^,  S{t)  is  readily  evaluated 
numerically  using  a  Gauss-Laguerre  quadrature  rule.  The  result  is  depicted  in  Figure 
14.  It  is  observed  that  for  r  <  0.005,  S{t)  becomes  constant  leaving  vi{t)  (x  and 
independent  of  conductivity.  In  particular,  the  early-time  asymptote  is  v{t)  =  noal /2t 
(hm^^o5'(T)  =  1/407). 

More  signihcantly,  for  r  >  10,  S{t)  (X  t  yielding  v{t)  x  t  The  computed 
response  is  depicted  in  Figure  15  for  three  values  of  adc-  Absolute  signal  voltage 
is  normalized  for  arbitrary  source  current  I  (and  by  displaying  the  absolute  value 
of  the  response,  we  effectively  chart  the  response  for  the  typical  case  of  a  step-off 
source  current).  As  predicted,  transition  from  to  dependence  occurs  at 

later  time  with  increasing  soil  conductivity.  Note  for  adc  =  0.1  mS/m  asymptotic 
early  time  dependence  is  particulalry  evident  with  transition  to  dependence 
at  approximately  t  =  10“^*^  —  10“^  ns. 


Employing  a  power  series  expansion  for  the  Bessel  function  in  Equation  (67)  and 
integrating  term  by  term,  Raiche  and  Spies  (1981)  obtained  the  following  late-time 
approximation  for  S{t) 


1  ^  (-1)^  (2fc+2)!  mV 

{k -\- 1)\  (k -\- 2)\  {2k -\- h)  0'^/ 


(68) 


Retaining  only  the  initial  term  for  k  =  0,  one  obtains  the  common  late-time  (r  >  10) 
approximation 


The  foregoing  approximation  V2{t)  is  compared  with  the  previous  result  vi{t)  in  Fig¬ 
ure  16  for  t  >  10“^°  s.  Clearly,  for  the  majority  of  metal  detectors  with  measurement 
windows  in  the  range  of  ten  to  several  hundred  microseconds,  this  late-time  approxi¬ 
mation  is  perfectly  adequate. 

3.3.3  Non- Conductive,  Non-Dispersive  Magnetic  Soil: 

For  a  non-conductive  soil  ai  =  0,  we  have  71  =  0,  =  A  and,  consequently, 

1  ,  1  ,  MoA(//r  -  1)  1  ,  X 

1  -I-  rTE  =  1-1 - — - —  =  1-1 - - 

PoX{pr  +  1)  X  +  2 
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where,  //r  =  Mi/Mo  =  Mdc/^o  represents  the  relative  magnetic  permeability  and  y  = 
Hr  —  lis  the  volnme-specihc  magnetic  snsceptibility.  In  particnlar,  for  a  non-dispersive 
magnetic  soil,  x  =  Mdc/ho  —  1  =  Xdc  and  Eqnation  (63)  yields 

v{t)  =  -noTTo^I  - S{t)  [  [Ji(Ao)]^  dA,  (71) 

V  Xdc  +  Jo 

where  the  Dirac  delta  fnnction  d{t)  is  dehned  as  follows 


d(t)  =0;  t  7^  0 


1. 


As  anticipated,  the  response  for  a  non-condnctive,  non-dispersive  magnetic  soil  is 
limited  to  an  impnlse  at  the  instant  of  sonrce  cnrrent  termination  (t=0).  There  is  no 
snstained  transient  response.  In  practice,  becanse  detection  circnitry  is  designed  to 
avoid  the  primary  switching  transient,  no  signal  is  registered. 


Das  (2004,  2006)  treated  the  more  general  case  of  coaxial,  coplanar  coils  at  height  h 
above  the  air-soil  interface.  In  general,  the  stand-off  dependent  conpling  coefficient  is 


M{h)=  I  Ji(Ao)  Ji(A6)e“2^^dA 


2 


(72) 


where  =  4o6/[(o  -|-  6)^  =  4h^]  and  K  and  E  represent  complete  elliptic  integrals 
of  the  first  and  second  kind  respectively.  For  the  specific  case  of  snrface-deployed 
{h  =  0)  coincident  (6  =  o)  coils,  we  have 

pOO 

M{0)=  /  [Ji{Xa)]^  dX  =  oo 

Jo 

Again,  for  instantaneons  sonrce  cnrrent  termination,  the  impnlsive  response  is  infinite. 
In  practice,  of  conrse,  the  response  is  finite  for  finite  termination  time.  The  main 
finding,  however,  is  that  a  non-condnctive  and  non-dispersive  magnetic  soil  yields  no 
snstained  transient  response.  The  same  cannot  be  said  for  a  dispersive  or  viscons 
magnetic  soil. 


(73) 


3.3.4  Non-Conductive,  Dispersive  Magnetic  Soil: 

As  Das  (2004)  has  demonstrated,  a  nsefnl  approximation  for  the  response  over  a  non- 
condnctive,  dispersive  magnetic  soil  is  obtained  by  generalizing  the  previons  resnlt  for 
a  frequency-dependent  magnetic  susceptibility.  In  particular,  a  well  established  model 
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for  susceptibility  dispersion  of  a  viscous  magnetic  soil  (Richter,  1937;  Chikazumi, 
1965;  Lee,  1984)  assumes  a  log-uniform  distribution  of  grain-volume  related  relaxation 
constants.  The  corresponding  magnetic  susceptibility  is  given  by 


Xif)  Xdc 


1  f  1  +  i27r/T2  \ ' 

ln(T2/Ti)  ^  Vl +  ^27r/TiyJ  ’ 


(74) 


where  ti  and  T2  denote,  respectively,  lower  and  upper  band-limits  on  the  time-constant 
distribution.  For  ti  =  —  At,  T2  =  t^  +  At  and  in  the  limit  At  — )■  0,  Equation  (74) 

reduces  to  the  standard  Debye  dispersion  relation  (Debye,  1929) 


xif) 


Xdc 

1  +  i2ir/T^ 


(75) 


and  the  accompanying  transient  decay  of  magnetization  is  purely  exponential.  In 
general,  for  an  arbitrary  time-constant  distribution  7)^(t),  with  T^i^)  dr  =  1,  we 
have 


xif)  Xdc 


TAr) 

1  +  i27r/T 


dr, 


(76) 


which  yields  Equation  (74)  for  the  specihc  1/T-scaled  logarithmic  time-constant  dis¬ 
tribution 


r  0,  T  <  Ti; 

T'xi^)  =  '7)t(lnT)/T  =  <  l/Tln(T2/Ti),  Ti  <  T  <  T2;  ,  (77) 

I  0,  T  >  T2. 


The  nature  of  related  frequency  dependence  is  displayed  in  Figure  17  for  log(T2/Ti) 
ranging  between  0.0  ixFD  ~  90%)  and  100.0  ixFD  ~  2%).^^  In  general,  together  with 
available  empirical  evidence  (e.g.  Dearing,  et  ah,  1996),  the  foregoing  model  suggests 
that  a  relatively  broad  time-constant  bandwidth  log(T2/Ti)  >  10  ixFD  <  15  —  20%) 
is  the  practical  reality.  Indeed,  it  is  well  recognized  (Nagata,  1961;  Dunlop  and 
Ozdemir,  1997)  that  relatively  minute  variation  in  the  grain  volume  of  magnetic 
particles  (in  proximity  of  the  stable  single-domain/superparamagnetic  transition)  is 
associated  with  a  comparatively  enormous  swing  in  related  viscous  time-constants. 


Now,  noting  that  for  the  vast  majority  of  soils  x  <<  2,  we  have  from  Equation  (70) 


1  +  T‘te  1  + 


xif) 

2 


(78) 


A  commonly  referenced  measure  of  frequency  dependence  of  magnetic  susceptibility  is  defined 
3-8  XFD  —  \iXLF  ~  Xhf) /Xlf\  X  100,  where  XLF  and  XFLF  denote  low-frequency  and  high- 
frequency  susceptibility  measured  at  frequencies  and  fjj,  spanning  a  decade.  See  Figure  17.5 
for  connection  between  XFD  and  common  dispersion  models,  including  log-uniform  (Richter,  1937  - 
Equation  74),  Cole-Cole  (Cole  and  Cole,  1941)  and  log-normal  (Wagner,  1913;  Yager,  1936)  models. 
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and,  therefore, 


v{t)  ^  —no7:a‘^IM{0) 


from  Equation  (63).  Finally,  substituting  s 
transforms  yields 


i  e"'/'  df 

(79) 

=  i27Tf  and  applying  standard  Laplace 


1  /  1  +  t27r/T2  \ 

ln(T2/Ti)  ^  V 1  +  i27r/Ti  / 


n(t)  —hqtto^I M.{0) 

and  assuming  ti  «  t  «  T2, 


m 


Xdc 


21n(T2/Ti)  t 


1  _  p-t/T-2 


(80) 


n(t)  —fj,Q7:a^I M{0) 


Xdc 


21n(r2/Ti)  (' 


(81) 


Note  that  in  addition  to  an  instantaneous  primary  impulse  associated  with  source 
current  termination,  magnetic  viscosity  leads  to  a  sustained  transient  response  having 
dependence.  Thus,  in  comparison  with  late-time  dependence  due  to  eddy 

currents  in  a  non- magnetic  and  non-dispersive  conductive  soil  (Equation  69),  signal 
decay  is  far  more  gradual  for  a  viscous  magnetic  soil,  resulting  in  a  substantially 
enhanced  and  potentially  anomalous  response. 

Although  Ad(0)  =  oo  for  coincident  coils,  a  hnite  approximation  is  obtained  for  b  ^  a. 
In  particular.  Das  (2006)  compared  the  late-time  response  predicted  by  Equation  (81) 
with  that  for  a  purely  conductive  soil  as  per  Equation  (69).  As  previously  noted,  it 
was  observed  that  even  for  an  extreme  soil  conductivity,  a  viscous  magnetic  soil  with 
only  moderate  static  susceptibility  yields  a  substantially  greater  response  over  the 
typical  measurement  time  range.  We  shall  return  to  the  case  of  a  viscous  magnetic 
soil  but,  hrst,  we  consider  the  more  general  case  of  a  magnetic  and  conductive  soil. 


3.3.5  Non-Dispersive  Conductive-Magnetic  Soil: 

Lee  (1984)  treated  the  common  case  of  a  soil  having  both  hnite  electrical  conductivity 
and  magnetic  susceptibility.  Beginning  with  Equation  (30)  and  considering  a  single¬ 
turn  coil  of  hnite  cross-sectional  radius  g,  a  corresponding  held  relation  was  obtained 
for  step  source  current  with 


l  +  rTE=  - .  (82) 

PiX  +  PoUl 

The  corresponding  time-domain  response  is 


v{t) 


—  2^0770^1 


2p\\ 

PiX-\-  PqUi 


A2TT  ft 

- df  [Ji{\a)f  dAd\ 

T7Q^ 


(83) 
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where  A  represents  the  cross-sectional  area  of  a  coincident  receiver  loop  (condnctor). 

Assnming  a  non-dispersive  soil  with  /j,i  =  1  +  Xdc,  the  integral  with  respect  to  /  is 
evalnated  by  appropriate  contonr  integration  and  the  integral  over  the  cross-section 
of  the  coil  is  approximated  by  Trp^.  Finally,  the  integral  with  respect  to  A  is  evalnated 
by  expanding  the  Bessel  fnnctions  as  a  power  series  and  integrating  term  by  term. 
For  small  Xdc  Lee  obtained  the  following  late-time  (r  =  t/anoo^  »  1)  response  by 
Taylor  series  approximation  abont  Xdc  =  0  (see  also  Ignetik,  et  ah,  1985) 

lipy/^al  (  a\iQO?  \  _  19/roXdeV^Qd  /  a\iQO? 

20t  \  t  J  280t  V  t 


V{t)  =  V3{t)  ^ 


Comparing  this  result  with  Equation  (69),  it  is  observed  that  the  hrst  term  is  identical 
to  the  late-time  response  for  a  non-magnetic  soil.  Moreover,  the  second  term  associ¬ 
ated  with  frequency-independent  susceptibility  Xdc  is  characterized  by  an  equivalent 
time  dependence.  Thus,  as  Lee  suggested,  the  response  can  be  rewritten  as 


//oVFa/  /  _  y/ir  I  5/2 

20t  \  t  J  ~  20 


(85) 


where 


cra  =  o- 


19  Xdc 
14 


2/3 


(86) 


denotes  a  correspondingly  enhanced  apparent  conductivity.  In  effect,  the  influence 
of  non-dispersive  susceptibility  on  the  late-time  response  is  indistinguishable  from 
a  marginal  increase  in  electrical  conductivity.  Induced  magnetization  enhances  the 
decaying  inductive  held  by  a  factor  =  //i//ro  =  1  +  Xdc-  As  indicated  by  Equation 
(86),  however,  the  effect  is  minor  for  typical  soil  susceptibilities  and  the  resulting  late¬ 
time  response  does  not  deviate  appreciably  from  that  predicted  by  Equation  (69). 
Related  inhuence  is  illustrated  in  Figure  18  and  in  Figure  19  for  a  representative 
range  of  Xdc-  As  demonstrated  in  previous  sections,  the  principal  signature  of  non- 
dispersive  magnetic  susceptibility  is  an  impulsive  signal  accompanying  source  current 
termination  at  t  =  0. 


3.3.6  Conductive,  Dispersive  Magnetic  Soil: 

More  generally,  Lee  (1984)  considered  a  frequency-dependent  magnetic  permeabil¬ 
ity.  In  particular,  assuming  the  susceptibility  dispersion  model  in  Equation  (74)  and 
Figure  17,  we  have  the  related  soil  magnetic  permeability 


Ml 


Mo 


+  Xdc 


1 

ln(T2/Ti) 


In 


/ 1  -Fi27r/T2\l\ 

\1+  i27r/Ti  y  J  y  ’ 


(87) 


in  Equation  (82).  More  specihcally,  Lee  (1984)  expands  Equation  (82)  via  Taylor 
series  about  {pi  —  //o)/mo  and  evaluates  the  resulting  held  equation  via  methods 
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similar  to  those  employed  for  a  non-viscous  permeability.  The  resulting  late-time 
(r  =  t/anoo^  »  1)  response  is 


v{t) 


fioV^al  / (Ta/roo^  \ 

2ot  V  t  y 
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19  Mo  XdcV^al 

( cJMoo^ 

\  V’(5/2) 
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V  i  ) 

_ln(T2/Ti) 

ln(T2/Ti)_ 

2//0  XdcT^al  /  a noo?  \  _  //p  Xdc  al  ln(2  a/g) 
3tln(T2/Ti)  V  t  /  2t  ln(T2/Ti) 


(88) 


where  the  digamma  function  'ip{5/2)  ^  0.703157.  Note  that  the  initial  term  is  con¬ 
sistent  with  Equations  (85-86),  comprising  the  late-time  response  for  a  soil  having 
apparent  conductivity  aa,  including  the  influence  of  non-viscous  magnetic  suscepti¬ 
bility  (for  T2/T1  — )■  00,  Equation  (88)  reduces  to  Equation  (85))  Remaining  terms 
reflect  the  influence  of  magnetic  viscosity  and  related  frequency-dependent  permeabil¬ 
ity  as  described  by  Equation  (87).  In  general,  for  practical  coil  geometries,  the  coil 
radius  o  is  large  compared  with  its  cross-sectional  radius  g  and  the  late-time  response 
is  ultimately  dominated  by  the  hnal  term  with  dependence,  compared  with 
dependence  for  nondispersive  soils.  Thus,  as  previously  discussed,  magnetic  viscosity 
yields  a  substantially  enhanced  and  potentially  anomalous  response.  The  value  of 
a/g,  combined  with  the  electrical  conductivity  adc,  magnetic  susceptibility  Xdc  and 
the  viscous  magnetic  time-constants  ti,  T2,  determine  the  precise  transition  from 
to  dependence. 

The  predicted  response  is  displayed  in  Figure  20  for  Xdc  =  0.001,  log(T2/Ti)  =  20 
{Xfd  10%),  o/p=100  and  the  same  three  conductivity  values  in  Figure  18  and 
previous  hgures.  Decreasing  the  ratio  a/g  by  a  factor  of  10  only  marginally  reduces 
the  predicted  response.  In  general,  with  increasing  electrical  conductivity,  predicted 
transition  from  to  dependence  is  effectively  delayed  as  per  Equation  (88). 

However,  assuming  typical  coil  dimensions  for  hand-held  systems  and  practical  limits 
on  soil  electrical  conductivity,  the  viscous  magnetic  response  is  generally  predominant 
within  the  practical  measurement  range  {t  >10  /rs).  We  shall  return  to  this  issue  in 
following  discussion.  First,  we  consider  the  relative  influence  of  specihc  magnetic 
parameters. 

While  magnetic  dispersion  is  an  essential  condition  for  sustained  decay.  Figure 
21  suggests  that  the  predicted  response  is  not  as  sensitive  to  the  degree  of  underlying 

12 

Related  expressions  (31)  and  (35)  in  Lee  (1984)  are  not  entirely  consistent,  with  (35)  (in  the 
limit  Ti  — >■  0)  implying  a  marginally  higher  apparent  conductivity  tj^^n 

(31)  and  current  Equations  (86)  and  (88).  The  effect  is  indistinguishable  on  the  scale  of  current 
diagrams  and  insignificant  in  the  present  context. 
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dispersion  as  it  is  to  the  absolute  level  of  magnetic  susceptibility.  In  particular,  Figure 
21a  displays  the  influence  of  increasing  time-constant  bandwidth  from  log  (T2/T1)  =  10 
(xfd  >  15%)  to  log  (T2/T1)  =  100  {xfd  ~  2%),  a  representative  range  for  the  major¬ 
ity  of  magnetic  soils  that  are  nominally  viscous.  It  is  observed  that  the  corresponding 
order  of  magnitude  reduction  in  related  response  is  relatively  minor  compared  with 
the  effect  of  varying  the  static  (low-frequency)  susceptibility  Xdc  over  a  representative 
range  spanning  several  orders  of  magnitude  for  a  nominal  time-constant  bandwidth 
(Figure  21b).  Thus,  although  magnetic  dispersion  is  a  prerequisite  condition  for  the 
critical  transition  from  to  t~^  decay  rate,  the  onset  and  level  of  viscosity  dom¬ 

inated  response  is  more  strongly  influenced  by  absolute  magnetic  susceptibility.  In 
effect,  absolute  susceptibility  scales  underlying  viscosity  to  yield  net  frequency  de¬ 
pendence  that  can  be  considerable  despite  a  relatively  modest  Xfd-  Indeed,  this 
would  appear  to  conhrm  the  direct  influence  of  both  absolute  susceptibility  and  re¬ 
lated  frequency  dependence  (Ax  =  xlf  —  Xhf)  on  ground  reference  height  (GRH) 
for  a  standard  calibrated  metal  detector  as  reported  by  Guelle,  et  al.  (2006). 

Returning  to  Equation  (81),  it  is  useful  to  compare  the  foregoing  result  with  the 
late-time  approximation 

v{t)  =  v^{t)  ^  -//oTTO^I M(0)  2  ln^4%i)  T 


Note  that  in  place  of  A1(0)  =  cx)  as  predicted  by  Equation  (72)  for  6  =  o,  h  =  0,  Equa¬ 
tion  (88)  implies  A1(0)  =  ln(2o/p)/7ro.  The  resulting  approximation  is  illustrated  in 
Figure  22.  As  expected,  there  is  excellent  agreement  with  predicted  response 
over  the  practical  measurement  range  (10 //s  —  1000 //s).  Again,  although  Equation 
(89)  hinges  on  a  dispersive  magnetic  susceptibility,  the  expression  also  emphasizes 
the  modulating  and  potentially  predominant  influence  of  absolute  susceptibility  Xdc 
(largely  reflecting  composition  and  concentration  of  the  viscous  soil  magnetic  frac¬ 
tion). 

Equating  expressions  (89)  and  (85),  we  obtain  the  following  approximate  relation 


tv  ~  Mo 


Vf  ln(T2/Ti) 

10  Xdc  ln(2o/p) 


(90) 


for  the  transition  time  from  to  dependence  for  a  conductive  and  viscous 

magnetic  soil.  Although  the  resulting  expression  inevitably  overestimates  the  actual 
transition  time,  it  is  an  adequate  approximation  and  provides  additional  insight  on 
related  influence  of  speciflc  soil  electromagnetic  parameters.  In  particular.  Figure  23 
illustrates  the  offsetting  influence  of  time-constant  bandwidth  log(T2/Ti)  and  absolute 
low-frequency  susceptibility  Xdc- 
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Given  well  established  bounds  on  Xdc  and  accepting  log(T2/Ti)  >  10  {xfd  <  15  — 
20%)  as  an  empirical  lower  limit,  it  remains  of  interest  to  consider  a  practical  upper 
limit  on  time-constant  bandwidth  and  related  implications.  In  theory,  so  long  as  T2/T1 
remains  hnite,  the  corresponding  response  (Equation  88)  continues  to  be  viscosity 
dominated  and  display  dependence  for  t  >  ty-  However,  as  practical  matter, 
increasing  log(T2/Ti)  requires  a  proportional  rise  in  Xdc  to  maintain  a  given  ty  or 
related  signal  level  and,  consequently,  effective  magnetic  viscosity  v  Xdc/ log  (T2/T1) 
is  the  prefered  parameter  for  gauging  the  net  influence  of  soil  magnetic  dispersion. 

Equation  (90)  also  indicates  that  transition  time  ty  scales  with  coil  radius  in  approx¬ 
imate  accordance  with  a  standard  quasi-static  scaling  relation  (e.g.  Frischknecht, 
1987) 

ty  ~  ty  ;  (91) 

where  ty  denotes  viscous  transition  time  for  scaled  coil  radius  o'.  Figure  24  displays 
ty  as  a  function  of  coincident  coil  radius  o,  with  low-frequency  conductivity  adc  as  a 
parameter.  Trendlines  are  flt  over  the  range  o  =  0.1  —  1.0  m.  Deviation  (reduction) 
from  standard  squared  dependence  on  o'/o  increases  as  the  ratio  a/g  increases  with 
decreasing  coil  radius.  As  previously  noted,  however,  the  influence  is  subtle. 

The  corresponding  scaling  relation  for  conductivity  (cr  =  cxa  ~  Cdc)  is 

ty  =  ty  f ^  •  (92) 

\^dc  J 


3.3.7  Non-Magnetic,  Dispersive  Conductive  Soil: 

Finally,  Lee  (1981)  and  El-Kaliouby,  et  al.  (1995,  1997)  have  investigated  the  in¬ 
fluence  of  frequency-dependent  electrical  conductivity  on  the  TDEM  response  of  a 
surface- deployed  coincident-coil  system.  Assuming  a  non-magnetic  soil,  we  have  as 
per  Equations  (35)  and  (64) 


1  -|-  tte  =  1  T - -  (93) 

A  +  mi  ^  ^ 

In  practice,  magnetic  viscosity  is  estimated  as  U  Ax/log(/H//L),  with  Ax  =  xlf  - 
XHF  denoting  the  differential  susceptibility  measured  between  measurement  frequencies  /x,  and 
fdj.  In  particular,  for  measurement  frequencies  spanning  a  decade  {ffl/ fu  —  10)i  ^  Ax- 
Owing  to  normalization  by  XLFi  related  parameter  XFD  —  ^X/XLF  ^  100  provides  improved 
discrimination  of  intrinsic  dispersion,  related  to  grain-size  distribution  of  soil  magnetic  content,  but 
fails  to  reflect  the  net  extent  of  effective  magnetic  viscosity  (see  also  Sections  3.4  and  4.5.1). 
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with  ui  =  [A^  +  i27r///oC’'(/)]^^^-  Substituting  Equation  (93)  in  Equation  (63)  and 
ignoring  the  impulsive  primary  response  associated  with  termination  of  the  source 
current  at  t  =  0,  yields  the  step  response  due  to  soil 

v{t)  =  fioTTO^I  /  /  - - -  [Ji{Xa)]^  df,  (94) 

Jo  J-oo  A  +  Ml 


Lee  (1981)  evaluates  the  integral  with  respect  to  A,  obtaining  a  power  series  solution 
that  is  subsequently  integrated  with  respect  to  complex  frequency  /  =  -Re*^/27r  by 
way  of  appropriate  contour  integration. 

Assuming  a  Cole-Cole  model  for  electrical  conductivity  dispersion  (Cole  and  Cole, 
1941;  Pelton,  et  ah,  1978) 


^  —  ^dc 
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1  -I-  [1  -  m]  [t27r/Tcr]‘ 
the  resulting  asymptotic  integral  solution  is 
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(96) 


S(fi)  =  E 
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(97) 

R  is  the  contour  integration  variable  and  angle  ip  =  tt/ A  dehnes  the  branch  cut 
geometry.  In  addition,  we  have 


F{R)  = 


fJ'O^dcR 


1-1-25-1-7^ 
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an  =  sin(7r/2  +  ip)  +  ip  +  (2n  -|-  3)  (7r/2  +  ip  +  </>)/2, 
(3n  =  sin(7r/2  F  ip)  F  ip  F  (2n  -|-  2)  (7r/2  Fip  F  </>)/2, 


with 

7  =  {TaR)p 

5  =  7  cos  [(7r/2  F  ipi)c] , 

?7  =  7sin  [(7r/2  Fip)c], 

(1  —  m)  77 
1  -I-  (1  —  m)  5 

With  the  change  of  integration  variable  r  =  Rt  simp,  the  response  is  readily  evalu¬ 
ated  numerically  employing  a  Gauss-Laguerre  quadrature  formulation  (El-Kaliouby, 
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et  al.,  1995).  Results  presented  by  Lee  (1981)  for  specific  Cole-Cole  parameters,  were 
confirmed  independently  by  Raiche  (1983)  by  way  of  an  alternative  computational 
method.  As  a  check  on  our  numerical  evaluation  of  Equation  (96)  the  same  results 
were  reproduced  prior  to  investigating  response  characteristics  on  the  scale  of  hand¬ 
held  sensors.  In  general,  the  response  predicted  by  Equation  (96)  is  characterized 
by  a  late-time  polarity  reversal  and  diminished  decay  rate.  Specihc  characteristics 
depend  on  the  nature  and  degree  of  electrical  polarization  and  associated  dispersion. 


Frequency  dependence  of  the  electrical  conductivity  as  specihed  by  Equation  (95)  de¬ 
pends  on  the  static  value  adc,  the  corresponding  high-frequency  value  aoo  (or  charge- 
ability  m  =  ((Too  —  o'dc)/o'oo),  the  Cole-Cole  distribution  parameter  c  and  reference 
time-constant  as  illustrated  in  Figure  25. 


In  particular,  for  c 
relation 


1.0,  Equation  (95)  yields  the  related  Debye-like  dispersion 


cr{f)  =  CTc 


^  dc 


1  +  {1  —  m)  i2TTfTa 


(98) 


and  the  accompanying  transient  decay  of  induced  electrical  polarization  is  purely 
exponential.  For  an  arbitrary  time  constant  distribution  7^(t),  with  Tair)  dr  = 
1,  we  obtain  the  general  relation. 


^(/)  ^oo  (^oo  ^  dc)  I 

Jo 


TAr) 


1  -|-  (1  —  m)  i27TfT 
A  l/r-scaled,  logarithmic  time-constant  distribution  having  form 

7^(lnr)  1  sin[(l  — c)7r] 


dr. 


(99) 


r.(r)  = 


27rT  cosh  [cln(T/Tcr)]  —  cos  [(1  —  c)7r 


(100) 


yields  the  Cole-Cole  dispersion  relation  of  Equation  (95).  Figure  25  displays  predicted 
normalized  frequency  dependence  for  adc  =  -01  mS/m,  m  =  0.3,  =  0.0001  s  and 

for  the  Cole-Cole  distribution  parameter  c  ranging  between  0.01  and  1.0. 

The  late-time  response  predicted  by  Equation  (96)  for  a  non-magnetic,  electrically 
polarizable  soil  is  depicted  in  Figure  26  (together  with  previous  models)  for  three 
values  of  static  conductivity.  Note  that  predicted  polarity  reversal  occurs  between 
0.1  and  1.0  and  that  corresponding  later-time  decay  rates  are  substantially  re¬ 
duced  compared  with  the  dependence  observed  for  a  non-dispersive  soil.  As  Lee 
(1981)  noted  for  t/anoo^  >  10,  the  corresponding  late-time  response  is  well  approxi¬ 
mated  by  retaining  only  the  initial  (n  =  0)  term  of  the  asymptotic  series  of  Equation 
(97).  For  parameters  in  Figure  26,  the  deviation  between  single-term  and  multi-term 
approximations  is  minimal  for  t  >  1.0  /ds. 

For  the  specihc  case  adc  =  0.01  S/m,  Figure  27  provides  some  indication  of  the 
inhuence  of  individual  Cole-Cole  parameters  on  the  predicted  response.  For  all  but 
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the  most  rapidly  acting  polarization  processes  (r^r  >  10“®  s),  polarity  reversal  and 
related  decay  rate  rednction  occur  later  as  the  associated  time-constant  r^r  increases 
(Figure  27a).  Increasing  chargeability  m  has,  in  general,  the  opposite  effect  (Figure 
27b),  leading  to  earlier  polarity  reversal.  Dependence  on  the  Cole-Cole  distribution 
parameter  c  is  more  complicated.  For  specihc  parameter  values  identified  in  Figure 
27c,  polarity  reversal  occurs  at  minimum  time  for  c  ~  0.2  and  increases  for  both  lesser 
and  greater  values.  The  general  nature  and  relative  influence  of  Cole-Cole  parameters 
on  polarity  reversal  time  tn  is  illustrated  in  Figure  28  for  representative  parameter 
ranges  and  reference  values  as  per  Figure  27. 

The  corresponding  rate  of  increase  in  related  decay  rates,  following  peak  negative¬ 
valued  response  is  also  signihcant.  Related  variation  with  Cole- Cole  parameter  values 
is  displayed  in  Figure  29.  In  general,  it  is  observed  that  polarization  parameters  have 
considerable  and  complicated  influence  over  the  transition  from  polarity  reversal  to 
late-time  dependence.  In  particular,  for  specific  parameter  values  considered. 

Figure  29  suggests  that  sustained  response  (low  decay  rate)  immediately  following 
polarity  reversal  is  associated  with  increasing  Ta,  moderate  m  and  increasing  c.  It  is 
also  observed,  however,  that  these  characteristics  do  not  generally  lead  to  an  early 
phase  reversal. 

Finally,  it  is  important  to  appreciate  that  response  characteristics  are  also  dependent 
on  coil  radius  and  conductivity,  with  polarity  reversal  occurring  later  as  the  product 
of  coil  radius  and  conductivity  increases.  Figure  30  displays  polarity  reversal  time 
tR  as  a  function  of  coincident  coil  radius  o,  with  low-frequency  conductivity  adc  as  a 
parameter.  Lower  curves  are  for  c  =  0.5  and  upper  curves  for  c  =  0  with  remaining 
parameters  as  per  the  reference  model  in  Figure  27.  Results  imply  a  generalized 
scaling  relation 

4  =  tR  ,  (101) 

where,  denotes  the  reversal  time  for  scaled  coil  radius  o'.  Note  that  for  an  in- 
hnitely  broad  time  constant  distribution  (c  =  0),  Equation  (101)  reduces  to  the 
standard  quasi-static  relation  =  tR{a'/a)‘^  with  squared  dimensional  dependence 
(e.g.  Frischknecht,  1987)  as  Equation  (96)  reverts  to  Equations  (66-69)  for  a  non- 
polarizable  soil  with  low-frequency  conductivity  a  =  2a dd^  —  RT')  (Lee,  1981).  The 
corresponding  scaling  relation  for  conductivity  follows  as 

\^dc  / 

The  principal  hnding,  however,  is  that  induced  electrical  polarization  produces  a  sig¬ 
nal  polarity  reversal  and  related  decay  rate  reduction  within  the  typical  measurement 
range  for  metal  detectors  and  related  small-scale  TDEM  sensors.  Moreover,  it  is 
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evident  on  the  basis  of  theory-predicted  results  in  foregoing  Figures  21,  26  and  27 
that  specihc  soil  electromagnetic  parameters  can  potentially  combine  to  yield  an  ef¬ 
fective  response  that  is  dominated  by  electrical  polarization  despite  limited  magnetic 
dispersion. 

3.4  Summary: 

Metal  detectors,  both  frequency  (FDEM)  and  time-domain  (TDEM)  varieties  are 
signihcantly  influenced  by  soil  electromagnetic  properties  and  predominantly  by  soil 
magnetic  characteristics.  In  general,  while  FDEM  systems  are  principally  affected  by 
the  absolute  levels  of  soil  magnetic  susceptibility  and  electrical  conductivity,  TDEM 
systems  are  more  strongly  influenced  by  related  dispersion  and  associated  viscosity. 

In  particular,  where  soils  incorporate  a  substantial  magnetic  component  having  vis¬ 
cous  susceptibility,  the  late-time  TDEM  response  is  generally  enhanced  with  decay 
rate,  compared  with  for  a  non-magnetic  or  non-viscous  magnetic  soil. 

With  few  exceptions,  the  model  of  soil  magnetic  dispersion  described  by  Equations 
(74)  and  (87),  and  leading  to  Equations  (88)  and  (89)  response),  appears  to  be 
generally  supported  by  available  data.  In  particular,  recent  multi-frequency  measure¬ 
ments  of  soil  magnetic  susceptibility  (West  and  Bailey,  2005;  Preetz  and  Igel,  2005; 
and  see  Section  4.5.1)  are  largely  consistent  with  the  following  approximate  relations 
(for  I/27rT2  <<  /  <<  I/27rTi,  ti  <<  t  «  T2) 

Xdc  _  dx'if)  _  2  x"{f)  _ _ I  dM{t)  _  ^ 

log(T2/Ti)  d\og  f  TT  loge  Hodlogt 

(Mullins  and  Tite,  1973;  Dabas,  et  ah,  1992),  where  the  right-most  expression  dehnes 
magnetic  viscosity  u  as  the  rate  of  change  of  time-dependent  magnetization  M{t) 
normalized  by  the  primary  source  held  HQ{t). 

General  validity  of  foregoing  relations  and  well-established  empirical  limits  on  xfd 
rehect  natural  grain-size  variation  and  imply  a  correspondingly  broad  time-constant 
distribution  (log(T2/Ti)  >  10).  It  is  also  observed  that  the  magnitude  of  soil  mag¬ 
netic  viscosity  depends  on  both  the  extent  of  underlying  dispersion  (grain-size/time- 
constant  distribution)  and  the  absolute  susceptibility  (composition  and  concentration) 
of  related  soil  magnetic  content  (i.e.  u  Xdc/ log  (T2/T1)). 

Signihcantly,  it  is  vastly  more  common  to  encounter  a  highly  magnetic  soil  with  low 
viscosity,  than  a  highly  viscous  soil  with  low  susceptibility.  The  explanation  is  that 
substantial  viscosity  {xfd  >  2%)  is  largely  attributable  to  a  signihcant  hne-grained 
superparamagnetic  fraction  that  generally  carries  a  correspondingly  higher  intrinsic 
magnetic  susceptibility  (Maher,  1988;  Forster,  et  ah,  1994;  Dearing,  et  ah,  1996).  In 
contrast,  a  soil  incorporating  a  substantial  concentration  of  stable  single-domain  or 
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multi- domain  magnetic  material  can  possess  a  considerable  magnetic  susceptibility 
with  limited  or  negligible  viscosity. 

In  theory,  only  hnite  magnetic  dispersion  is  required  to  produce  the  sustained  and 
characteristic  response.  However,  where  magnetic  dispersion  is  limited  {xfd  < 
2%),  related  signihcance  is  increasingly  dependent  on  the  level  of  associated  suscep¬ 
tibility  (concentration  of  related  viscous  magnetic  constant)  to  support  or  effectively 
amplify  the  related  signature.  Where  viscous  magnetic  content  is  limited  and  re¬ 
lated  susceptibility  is  low,  there  is  increasing  potential  for  the  related  signature  to 
be  superceded  by  the  background  conductive  response  (i.e.  tv,  given  by  Equa¬ 

tion  (90)  exceeds  the  effective  measurement  time).  Clearly,  in  the  limiting  case  of 
a  non-dispersive  magnetic  susceptibility,  no  amount  of  associated  magnetic  content 
yields  dependence  and  whatever  level  of  susceptibility  only  produces  an  effective 
enhancement  of  electrical  conductivity  and  related  response  as  per  Equations  (85-86). 

Finally,  foregoing  hndings  demonstrate  that  appreciable  dispersion  of  soil  electri¬ 
cal  conductivity  can  also  produce  anomalous  and  sustained  response  characteristics 
within  the  typical  measurement  range  of  hand-held  TDEM  systems.  In  particular, 
signihcant  electrical  chargeability  yields  a  characteristic  signal  polarity  reversal  and 
time-dependent  decay  rate  reduction.  Numerical  modelling  on  the  basis  of  Equations 
(96-97)  predicts  dependence  (Figure  29).  Most  signihcantly,  results 

suggest  that  for  soils  with  sufficient  chargeability,  related  induced  polarization  could 
potentially  dominate  the  TDEM  response  for  soils  having  non-dispersive  or  limited 
magnetic  susceptibility. 

In  general  additional  work  is  required  to  characterize  the  practical  extent  and  signif¬ 
icance  of  low-frequency  electrical  dispersion  in  soils.  However,  limited  available  data 
(e.g.  Ogilvy  and  Kuzmina,  1972;  Mehran  and  Arulanandan,  1977;  Iliceto,  et  ah,  1982; 
Olhoeft,  1985,  1987)  suggest  that  soil  electrical  chargeability  is  rarely  sufficient  to  ex¬ 
ceed  the  influence  of  even  moderate  magnetic  viscosity  (see  further  related  discussion 
in  Section  5.0). 

In  general,  the  remainder  of  the  current  report  is  devoted  to  a  brief  review  of  rele¬ 
vant  theory  and  methods  for  measurement  and  characterization  of  low-frequency  soil 
electromagnetic  properties. 
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4.0  MEASUREMENT  OF  SOIL  EM  PROPERTIES: 


As  per  foregoing  discussion  of  metal  detectors  and  the  basic  principles  of  induction 
electromagnetic  measurement  systems,  it  often  useful  and  instructive  to  model  elec¬ 
tromagnetic  characteristics  of  the  soil  in  terms  of  an  equivalent  electric  circuit.  As 
we  move  from  held  sensors  to  the  laboratory,  a  circuit  theory  perspective  is  less  an 
analogy  and  more  a  reality.  In  general,  laboratory  measurement  of  low-frequency 
soil  electromagnetic  properties  employs  a  range  of  sample  holders  or  hxtures  that 
are  effectively  large-scale  equivalents  of  fundamental  circuit  components,  including 
resistors,  capacitors  and  inductors.  The  sample  soil  simply  replaces  the  equivalent 
resistive,  dielectric  or  magnetic  material  within  the  component’s  interior.  The  hxture 
is  subsequently  driven  by  a  controlled  and  calibrated  source  signal  (voltage/current) 
with  simultaneous  detection  of  one  or  more  response  signals.  Finally,  appropriate 
signal  conditioning  and  processing  is  applied  to  source  and  response  signals  to  derive 
required  signal  attributes  and  to  compensate  for  hxture  and  measurement  system 
response  characteristics  prior  to  display  and/or  storage  of  measured  data. 

Again,  measurements  can  be  acquired  in  either  time  or  frequency  domain  and  a  con¬ 
siderable  range  of  laboratory  instrumentation  has  been  reported.  Due  to  the  relative 
instrumental  complexities  of  time-domain  measurements,  laboratory  characterization 
of  soil  electromagnetic  parameters  is  largely  carried  out  by  applying  a  continuous  si¬ 
nusoidal  source  signal  over  a  range  of  discrete  frequencies.  Typically,  the  relative  gain 
and/or  phase  of  a  resulting  response  signal  is  quantihed  with  reference  to  the  source 
signal  to  yield  a  vector  measurement  the  sample’s  complex  impedance  or  other  related 
parameters.  In  some  cases,  sample  related  deviation  of  hxture’s  resonant  frequency 
or  related  phase  characteristics  are  detected  to  yield  equivalent  effective  impedance 
parameters. 

Although  consideration  is  also  given  to  appropriate  technologies  for  held  measurement 
of  soil  electromagnetic  properties,  the  current  study  has  largely  focused  on  laboratory 
characterization.  In  addition  to  assessing  the  performance  of  selected  special-purpose 
instrumentation,  our  primary  approach  has  been  to  investigate  the  use  of  standard 
off-the-shelf  impedance  analyzers,  multi- frequency  LCR  meters  and  similar  widely 
available  technologies  for  measurement  of  both  electrical  and  magnetic  soil  properties. 
In  following  sections,  relevant  measurement  theory  is  developed  in  connection  with  the 
simplest  of  sample  hxture  conhgurations  to  demonstrate  the  practical  characterization 
of  soil  electromagnetic  properties  by  measurement  of  related  lumped  circuit  impedance 
parameters.  We  begin  by  dehning  the  frequency-dependent  complex  impedance. 
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4.1  Complex  Impedance- Admittance: 

In  general,  the  complex  impedance  of  an  arbitrary  two-terminal  network  is  defined  as 


Z{f)  =  R  +  iX, 


where  R{f)  denotes  the  eqnivalent  resistance  and 


X{f)  =  27r/L 


1 


(104) 


(105) 


represents  the  eqnivalent  reactance,  comprising  inductive  and  capacitive  contributions 
associated  with  equivalent  inductance  L  and  capacitance  C,  respectively.  Alterna¬ 
tively,  as  indicated  in  Figure  31,  the  complex  impedance  may  be  written  in  polar 
form 

Z{f)  =  \Z{f)\e^^^f\  (106) 

with  associated  frequency-dependent  magnitude 

|Z(/)|  =  (.R2  ^X2)1/2  (107) 


and  corresponding  phase  angle 


»(/)  =  tan  ' 


(108) 


relative  to  the  resistive  axis  in  the  complex  impedance  plane.  Note  also  that  while 
capacitive  reactance  leads  the  resistance  by  90  degrees,  inductive  reactance  lags  by 
90  degrees. 

The  ratio  of  resistive  to  reactive  components  of  the  complex  impedance 

R 

D  =  tan  5  =  —  (109) 

X 

is  referred  to  as  the  dissipation  factor  or  loss  tangent  tan  S,  where  S  is  the  complement 
of  the  phase  angle  9.  The  reciprocal  quantity  Q  =  1/D  is  termed  the  quality  factor. 
Finally,  in  connection  with  the  combination  of  impedances  in  parallel  (Figure  31),  it 
is  advantageous  to  make  use  of  the  inverse  impedance  or  the  admittance 


y{f) 


1 


1 

R  +  iX 


R-iX 


G  +  iB, 


where  G{f)  =  1/R  denotes  the  equivalent  conductance  and 


B{f)  =  27:  fC 


1 


(110) 


(111) 


35 


represents  the  equivalent  susceptance,  comprising  capacitive  and  inductive  contribu¬ 
tions  associated  with  parallel  equivalent  capacitance  C  and  inductance  L,  respectively. 
Note  that,  in  general,  measured  impedance  parameters  are  frequency  dependent  and 
typically  impure.  Consequently,  as  for  basic  material  parameters,  we  refer  to  these 
measured  quantities  as  equivalent  or  effective  parameters.  Figure  32  illustrates  basic 
phase  relations  for  pure  and  combined  circuit  elements. 

4.2  Electrical  Parameters  -  Theory: 

In  practice,  electrical  and  magnetic  properties  of  a  given  soil  sample  are  measured 
independently  using  appropriate  hxture  conhgurations  to  selectively  enhance  the  re¬ 
spective  response.  In  particular,  while  soil  magnetic  properties  are  investigated  by 
applying  a  primary  magnetic  held,  soil  electrical  parameters  are  commonly  measured 
in  response  to  a  primary  electric  held. 

Expanding  on  the  treatment  by  von  Hippel  (1965),  we  consider  a  time- varying  source 
potential 

V  =  Voe^^^^\  (112) 

applied  to  the  terminals  of  a  parallel-plate  sample  hxture  as  pictured  in  Figure  33a. 
Prior  to  placing  soil  in  the  hxture,  the  applied  potential  produces  a  charging  current 

where  Q  =  CqV  denotes  the  charge  developed  on  a  given  plate  electrode  and  Cq  =  7ceo 
is  the  geometric  capacitance  of  the  air-hlled  sample  hxture,  with  eo  representing  the 
electrical  permittivity  of  free-space  and  7c  =  A/ 1  a  geometrical  constant  given  by  the 
ratio  of  cross-sectional  area  (A)  to  separation  {1)  between  electrodes.  It  follows  from 
Equation  (113)  that  the  charging  current  leads  the  applied  potential  by  7r/2  radians 
or  90  degrees  (Figure  32).  Finally,  neglecting  fringing  ehects  at  edges,  the  associated 
electric  held  within  the  hxture  volume  is 

E  =  -VE  =  V/l  =  (114) 

^0 

where  q  =  QjA  denotes  the  time  dependent  surface  charge  density  on  a  given  plate 
electrode. 

Now,  for  a  sample  hlled  soil-box,  the  ehective  capacitance,  C  =  (e/eo)C'o  =  eye, 
is  enhanced  by  the  relative  electrical  permittivity  =  e/eo,  where  e  represents  the 
intrinsic  soil  permittivity.^^  Thus,  from  Equation  (113),  we  obtain  the  following 

Note  the  corresponding  dielectric  constant  is  the  special  case  defined  as  K  =  e^(0)/eo,  where 
e^(0)  denotes  the  static  real  electrical  permittivity  and  Cq  =  8.85  X  10  F/m  denotes  the 
electrical  permittivity  of  free  space. 
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expression  for  what  is  effectively  the  displacement  cnrrent  associated  with  bonnd- 
charge  polarization  processes  within  the  sample  soil 


/e  =  t27r/e7cV  =  (115) 

where 

^  dE 

is  the  associated  displacement  cnrrent  density.  In  addition,  there  is  a  free-charge 
condnction  cnrrent  predicted  by  Ohm’s  law 


I.= 


V 

R 


V 


-V  =  7,aV, 
P 


(117) 


where  R  =  pi/ A  denotes  the  effective  electrical  resistance,  p  is  the  electrical  resis¬ 
tivity  and  a  =  1/p  is  the  corresponding  electrical  condnctivity.  In  contrast  with 
the  displacement  cnrrent,  the  condnction  cnrrent  and  associated  condnction  cnrrent 
density 

=  ^  =  aE,  (118) 

within  the  sample  volnme,  are  in-phase  with  the  sonrce  potential  (Fignre  32). 
Finally,  combining  Equations  (115)  and  (117)  we  obtain  the  net  circuit  relation 


I  =  la  E  Ie  =  E  i‘27Tfe)  V 


(119) 


and  from  Equations  (116)  and  (118)  the  corresponding  held  relation 


J  =  Ja  E  Je  =  (cr  +  t27r/e)  E. 


(120) 


Now,  as  previously  described,  it  is  generally  observed  that  soil  electrical  properties 
display  significant  frequency  dependence  and,  consequently,  a  and  e  must  be  complex¬ 
valued  parameters 

a  =  a'  +  ia"  e  =  e'  -  ie" .  (121  -  122) 

Therefore,  we  have 


a  -\-  i27Tfe  =  {a  +  ia")  +  {e  —  ie") 

=  {(7'  E  277 fe")  +  t27r/ 

and,  as  a  result.  Equations  (119)  and  (120)  become 

7  =  7c  (d  +  i27Tfe)  V,  (123) 
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where 


J  =  (cj  +  i27Tfe)  E, 


(124) 


a  =  a  +  277 fe"  e  =  e  -\ - -  (125  —  126) 

27TJ 

denote  real-valued  composite  effective  parameters  (as  per  Section  2.0  Equations  13- 
14).  In  particular,  while  complex  conductivity  and  complex  permittivity  are  associ¬ 
ated  with  free  and  bound  charge  processes,  respectively,  real-valued  effective  conduc¬ 
tivity  and  effective  permittivity  are  similarly  associated  with  loss  (power  dissipation) 
and  polarization  (energy  storage)  processes.  For  example,  it  is  observed  that  viscous 
rotational  displacement  current  associated  with  polarization  of  dipolar  molecules  con¬ 
tributes  in-phase  with  the  Ohmic  loss  current  associated  with  mobile  ion  transport. 

Similarly,  the  effective  polarization  of  the  mobile  charge  distribution  associated  with 
interfacial  polarization  processes  is  in-phase  and  indistinguishable  from  the  dielectric 
polarization  of  bound  atomic  and  molecular  charge.  It  is  only  possible  to  separate  and 
quantify  related  contributions  to  composite  effective  parameters  based  on  the  relative 
phase  of  associated  currents. 

In  practice,  it  is  also  important  to  appreciate  that  in  addition  to  reflecting  material 
parameters  representative  of  the  sample  soil,  related  in-phase  and  quadrature  signal 
components  also  include  parasitic  contributions  associated  with  hxture  hardware  and 
cabling  (Haruta,  2000).  The  influence  of  higher  order  electromagnetic  fields  is  as¬ 
sumed  to  be  negligible.  For  example,  magnetic  fields  associated  with  time  varying 
current  density  and  related  induced  eddy  currents  are  assumed  to  have  no  appreciable 
influence  even  for  strongly  magnetic  and/or  conductive  soils. 

Recognizing  Equation  (123)  as  a  combination  of  parallel  currents,  we  measure  in- 
phase  and  quadrature  components  of  the  net  current  /  with  gain-phase  reference  to 
source  potential  V  to  yield  the  corresponding  frequency-dependent  complex  admit¬ 
tance  Y{f)  =  I /V  over  a  range  of  discrete  source  frequencies  between  100  Hz  and 
100  kHz. 

For  the  present  case  of  a  sample-fllled  soil-box.  Equation  (123)  yields 

Y{f)  =  ^  =  G  +  iB  =  j^{a  +  i27Tfe),  (127) 

and,  therefore,  assuming  that  the  effective  susceptance  of  the  parallel-plate  fixture  is 
purely  capacitive,  B  =  i27TfC,  we  have^® 


a  =  a  +  27Tfe'  =  —  =  — ^ 
7c  'JcR 


1 

p 


(128) 


15 


Note  that  the  equivalent  parallel  capacitance  Cp  is  related  to  the  equivalent  series  capacitance 
Cg  by  Cp  =  (/s(l  ff  D^)  where  D  =  R/X  =  G/B  denotes  the  dissipation  factor. 
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(129) 


~  t  I 

€  —  €  -\- 


27r/ 


C 

7c' 


where  p  =  l/a  defines  the  corresponding  effective  electrical  resistivity. 

Again,  it  is  emphasized  that  measnred  complex  impedance  spectra  can  also  reflect 
snperimposed  response  signatnres  associated  with  fixtnre  hardware  and  cabling.  Elec¬ 
trode  polarization  phenomena  are  also  potentially  significant  sonrces  of  error.  Meth¬ 
ods  to  reduce  and  compensate  for  these  parasitic  responses,  including  four-terminal 
vs.  two-terminal  measurement  configurations,  are  discussed  further  in  Section  4.4. 
For  now,  we  turn  our  attention  to  the  measurement  of  soil  magnetic  properties. 


4.3  Magnetic  Parameters  -  Theory: 

In  direct  analogy  with  the  previous  section,  we  consider  the  same  time-varying  source 
potential 

V  =  Voe^^^^\  (130) 

in  this  case,  applied  to  the  terminals  of  a  solenoidal  sample  fixture  as  depicted  in 
Figure  33b.  Here,  in  accordance  with  Faraday’s  induction  law  V  =  Lo{dI/dt),  the 
applied  potential  produces  a  magnetizing  current 

where  Lq  =  7^//o  denotes  the  inductance  of  the  air-cored  solenoid,  7^  is  a  coil  specific 
geometric  constant  and  po  =  47r  x  10“^  H/m  is  the  magnetic  permeability  of  free- 
space.  In  particular,  7^  =  N‘^Aj£,  where  A  denotes  the  cross-sectional  area  of  the 
solenoid,  £  is  it’s  length  and  N  is  the  number  of  turns.  Finally,  note  that  in  contrast 
with  the  case  of  the  parallel-plate  fixture.  Equation  (131)  indicates  that  the  drive 
current  lags  the  applied  potential  by  7r/2  radians  or  90  degrees  (Figure  32). 

Now,  upon  introducing  a  coaxial  sample  core  into  the  fixture,  the  effective  inductance, 
L  =  {p/po)Lo  =  is  enhanced  by  the  relative  magnetic  permeability  pr  =  p/po 
and  assuming  an  associated  frequency-dependent,  complex- valued  soil  magnetic  per¬ 
meability 

p  =  p'  —  ip” ,  (132) 

we  obtain  the  general  relation 


1 

p'  —  ip" 

for  the  net  current,  within  the  sample-cored  solenoid.  By  way  of  appropriate  instru¬ 
mentation,  the  in-phase  and  quadrature  components  of  the  net  current  /  are  detected 
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and  measured  with  gain-phase  reference  to  the  source  potential  V  to  yield  the  corre¬ 
sponding  frequency-dependent  complex  impedance 

Z(f)  =  j=R  +  iX  =  27r/7^  (//"  +  i//')  (134) 

over  a  range  of  discrete  source  frequencies.  Assuming  the  reactive  impedance  of 
the  solenoid  hxture  to  be  purely  inductive,  X  =  27r/L,  soil  magnetic  parameters 
are  related  to  measured  equivalent  impedance  parameters  according  to  the  following 
expressions 

//'  =  —  = //o(x  +  1)  (135) 

7a. 

m"  =  ^  f  =  Mo  x",  (136) 

27r/7^ 

where  x'  and  x"  denote  corresponding  in-phase  and  quadrature  components  of  the  as¬ 
sociated  complex,  frequency-dependent  magnetic  susceptibility  x(f)  =  x' if)- 
Consequently,  complex  magnetic  susceptibility  components  follow  directly  from  the 
difference  between  related  impedance  parameters  measured  before  and  after  insertion 
of  the  sample  core  as 

,  _  AL  _AL 
^  7^//o  Lq 

„  _  AR  _  AR 

^  27r/7^//o  27r/Lo 

where  AL  =  L  —  Lq,  AR  =  R  —  Rq  and  Rq  denotes  the  resistance  of  the  air-cored 
solenoid. 

Note  that  a  magnetic  loss  current  associated  with  the  quadrature  permeability  (sus¬ 
ceptibility)  and  related  viscosity  is  phase  shifted  by  180  degrees  and  effectively  reduces 
the  apparent  resistive  loss  current  Ic  =  V/ Rc  the  solenoid’s  windings.  In  addition, 
there  is  a  potentially  signihcant  back  current  associated  with  the  development  of  in¬ 
ductive  eddy  currents  within  an  electrically  conductive  sample  core.  Finally,  although 
higher  order  electromagnetic  helds  are  generated  within  the  sample  soil  in  connection 
with  time- variation  of  related  secondary  helds,  related  inhuence  on  the  measured 
response  is  presumed  to  be  negligible. 

Finally,  it  is  noted  that  an  alternative  approach  is  to  design  a  solenoid  having  speci- 
hed  open-coil  resonant  frequency  /q  =  l/27r(LoC'o)^/^  (where  Cq  represents  the  coil’s 
equivalent  distributed  capacitance)  and  to  similarly  derive  the  permeability  (suscepti¬ 
bility)  by  measuring  the  deviation  A/*  =  /*  —  /q  ,  where  f*  =  l/27r(LC'o)^/^  denotes 
the  resulting  resonant  frequency  on  insertion  of  the  sample  core.  In  practice,  the  LC 
oscillator  circuit  is  rather  more  complicated,  however  improved  accuracy  and  precision 
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make  resonant  frequency  methods  the  standard  approach  for  dedicated  susceptibil¬ 
ity  meters  and  we  shall  return  to  these  methods  in  connection  with  description  of 
commercially  available  susceptibility  instruments.  First,  we  turn  our  attention  to  the 
practical  measurement  of  soil  electrical  conductivity. 

4.4  Electrical  Parameters  -  Measurement: 

As  discussed  in  previous  sections  (see  Section  2.0  and  including  footnote  1),  it  is  as¬ 
sumed  that  the  metal  detector  response  is  adequately  described  under  a  quasistatic 
assumption  that  displacement  currents  associated  with  electrical  permittivity  e  are 
negligible  compared  with  conduction  currents  associated  with  conductivity  a.  Conse¬ 
quently,  the  following  discussion  is  limited  to  measurement  of  the  electrical  conduc¬ 
tivity  and  the  corresponding  reciprocal  quantity,  electrical  resistivity  p. 

As  a  point  of  departure,  methods  investigated  for  laboratory  and  held  measurement 
of  electrical  conductivity  may  be  viewed  as  generalizations  and  extensions  of  ASTM 
Standards  G57-06  and  G187-05  (ASTM,  2005,  2006). 

4.4.1  Laboratory  Methods  and  Procedures: 

Laboratory  measurement  of  direct  current  soil  electrical  conductivity  is  carried  out 
routinely  in  connection  with  a  wide  range  of  agricultural  and  engineering  applications. 
Although  specihcations  vary  considerably,  the  sample  hxture  is  generally  referred  to 
as  a  soil  box  and  roughly  approximates  the  parallel  plate  capacitor  model  discussed  in 
Section  4.2.  In  particular,  we  have  adopted  a  standard  format  (Figure  34)  with  geom¬ 
etry  designed  such  that  electrical  resistance  measured  in  four-electrode  conhguration 
yields  soil  electrical  resistivity  directly  in  Ohm-cm  units  The  corresponding  design 
constraint  is  7c  =  Ajly  =  1.0,  where  A  denotes  the  area  of  end-plate  electrodes  and 
ly  =  xi  is  the  separation  between  interior  potential  electrodes.  Resulting  resistivity 
(conductivity)  follows  as  p  =  l/cr  =  jcP  =  P- 

Alternatively,  the  interior  potential  electrodes  can  be  removed,  leading  to  standard 
two-electrode  resistivity  p  =  l/cr  =  7ci?  =  {A/£)R  as  per  Equation  (128),  where 
£  =  xi  +  2x2  is  the  full  length  between  end-plate  electrodes.  Goncurrently,  however, 
it  is  important  to  appreciate  that  the  current  soil  box  design  is  inconsistent  with 
the  assumption  7c  >>  1.0  for  an  ideal  parallel  plate  capacitor  and,  consequently. 
Equation  (127)  for  the  associated  electrical  permittivity  is  invalid.  In  particular,  it 
can  be  demonstrated  that  for  small  7c,  the  basic  relation  C  =  7ce  no  longer  holds 
and,  rather,  the  effective  capacitance  approaches  a  constant  value  as  hxture  length 
increases.  As  noted  above,  however,  our  present  interest  is  limited  to  the  electrical 
conductivity  and  there  is  no  similar  limitation  on  the  relation  R  =  p/yc  =  I/TcC’"- 

Soil  box  fixtures  employed  for  investigations  reported  here  are  commercially  available  from  MC 
Miller  Inc.  (www.mcmiller.com).  Similar  purpose-built  fixtures  are  available  from  other  suppliers 
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The  primary  advantage  of  the  present  soil  box  design  is  the  inclnsion  of  interior  po¬ 
tential  electrodes  to  permit  fonr-electrode  measurements  to  minimize  the  influence 
of  electrode  polarization  effects  at  end-plate  electrodes  as  illustrated  in  Figure  35. 
Notice  that  while  a  simple  two-electrode  Impedance  (or  Admittance)  measurement 
given  by  the  vector  ratio  of  source  voltage  Vg  and  source  current  A ,  neglects  the  influ¬ 
ence  of  the  polarization  held  Ep  within  the  sample,  the  corresponding  four-electrode 
measurement  samples  the  net  effective  held  E  =  Eg  —  Ep  via  related  voltage  drop 
V  across  interior  potential  electrodes.  Finally,  it  is  important  to  emphasize  the  dis¬ 
tinction  between  “four-terminal”  and  “four-electrode”  measurements.  In  particular, 
it  is  noted  that  both  two-electrode  and  four-electrode  measurements,  as  previously 
described,  are  acquired  in  four-terminal  conhguration,  with  signal  current  and  signal 
voltage  detected  separately  in  parallel  circuits. 

Unfortunately,  existing  standards  for  soil  electrical  resistivity  characterization  typ¬ 
ically  call  for  measurement  of  static  or  “low-frequency”  resistivity  pdc  and,  conse¬ 
quently,  commercially  produced  electrical  resistivity  meters  invariably  employ  a  sin¬ 
gle  low-frequency  or  commutated  DC  voltage  source  with  no  capability  for  multi¬ 
frequency  measurements.  For  this  reason,  laboratory  measurements  of  frequency- 
dependent  electrical  resistivity  (conductivity)  have  generally  employed  impedance, 
LCR  or  gain-phase  analyzers  developed  for  electronic  test  and  measurement  applica¬ 
tions. 

In  general,  multi- frequency  impedance  analyzers  and  LCR  meters  are  conhgured  for 
four-terminal  measurement.  Figure  36  displays  a  generic  four-terminal  impedance 
analyzer  as  conhgured  for  both  two-electrode  and  four-electrode  soil  box  impedance 
measurements.  Unfortunately,  experimentation  and  related  investigation  reveals  that, 
in  practice,  the  design  of  common  impedance  analyzers  and  LCR  meters  is  not 
conducive  to  accurate  impedance  measurements  in  the  desired  four-electrode  con¬ 
hguration.  On  the  contrary,  LCR  meters  are  engineered  primarily  for  testing  and 
quality  control  on  two-lead  electronic  components  and,  consequently,  as  commonly 
conhgured^^,  the  four-terminal  measurement  circuit  assumes  only  minimal  contact 
impedance  Zh  =  Z^c  +  '^hp  between  He  (high  current)  and  Hp  (high  potential)  ter¬ 
minals  and,  similarly,  Zl  =  ZicEZip  between  Lc  (low  current)  and  Lp  (low  potential) 
terminals  (Haruta,  2000). 
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The  standard  design  is  an  autobalancing  bridge  circuit  and  is  considerably  more  complicated 
in  detail  than  portrayed  in  Figures  (37-38).  The  source  signal  generated  by  a  precision  oscillator  is 
applied  at  current  terminal  Hq  .  Resulting  vector  potential  drop  across  the  sample  volume  is  sensed 
at  the  Hp  terminal  with  reference  to  a  common  shield  conductor  maintained  at  virtual  ground  by 
a  feedback  controlled  balance  oscillator  with  sensing  at  the  L p  terminal.  Test  current  returning  via 
the  shield  conductor  is  gauged  by  monitoring  current  through  a  calibrated  range  resistance  at  the 
Lc  terminal  (see  Haruta,  2000  for  details). 
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The  corresponding  two-electrode  soil-box  confignration  is  depicted  in  Fignre  37.  In 
contrast,  as  illnstrated  in  Fignre  38,  the  fonr-electrode  soil-box  measnrement  conhg- 
nration  implies  Zu  >>  +  Z^p  and  Zl  »  Zic  Zip^  with  Zh  ^  Zl  on  the  same 

order  of  magnitnde  as  the  sample  impedance  Zs-  Unfortnnately,  the  foregoing  condi¬ 
tion  strongly  restricts  the  fnnction  of  common  fonr-terminal  impedance  analyzers  and 
LCR  meters  for  fonr  electrode  resistivity  (conductivity)  measurement  (particularly  for 
low-moisture  (resistive)  soil  condition). 

In  view  of  the  foregoing  limitation  on  LCR  meters  and  required  portability  for  held 
characterization  of  soil  electrical  characteristics,  we  have  simultaneously  investigated 
an  alternative  approach  employing  a  stand-alone  signal  generator  and  digital  multi¬ 
meters  (DMMs).  We  shall  refer  to  related  measurements  as  GVI  (generator- voltage- 
current)  measurements.  In  effect,  related  four-terminal  measurement  conhgurations, 
as  illustrated  in  Figure  39,  are  identical  to  those  for  the  integrated  impedance  ana¬ 
lyzer  with  the  advantage  that  voltage  and  current  measurement  circuits  are  entirely 
independent.  An  isolation  ampliher  is  added  at  the  voltmeter  input  to  avoid  related 
circuit  interaction,  particularly  at  high  frequencies. 

Relative  performance  characteristics  of  investigated  measurement  conhgurations  are 
illustrated  in  Figure  40  for  untreated  Vancouver  tap  water  and  for  both  75  ml  and 
250  ml  soil-box  hxtures.  Note  that,  in  general,  all  investigated  methods  yield  results 
in  good  agreement  with  water  conductivity  as  monitored  via  a  YSI  Model  34  conduc¬ 
tivity  meter.  The  observed  shift  between  75  ml  and  250  ml  hxtures  is  attributable 
to  a  related  shift  in  tap  water  conductivity.  Given  the  uniform  single-phase  nature 
of  water,  no  substantial  electrophysical  or  electrochemical  polarization  processes  are 
anticipated  within  the  sample  volume  and,  consequently,  it  is  assumed  that  observed 
frequency  dependence  is  attributable  to  a  combination  of  hxture,  cabling  and  instru¬ 
ment  effects.  Finally,  it  is  noted  that  a  conductivity  of  approximately  a  ~  2.5  mS/m 
[p  ~  400  Ohm-m)  is  representative  of  moderately  resistive  soils  for  which  parasitic 
signatures,  including  electrode  polarization  effects,  are  generally  more  pronounced. 

We  hrst  consider  results  for  two-electrode  measurement  conhgurations.  In  particu¬ 
lar,  as  anticipated,  two-electrode  LGR  measurements  conhrm  a  roughly  frequency- 
independent  resistivity  for  water  samples,  with  the  exception  of  a  minor  apparent  re¬ 
duction  at  the  upper  band  limit  of  100  kHz.  In  contrast,  corresponding  two-electrode 
GVI  measurements  display  a  pronounced  decrease  in  apparent  electrical  resistivity  for 
frequencies  approaching  and  exceeding  10  kHz.  Experimentation  conhrms,  however, 
that  observed  effects  are  related  to  input  impedance  characteristics  of  the  voltmeter 
(DVM)  and  are  largely  eliminated  by  isolating  the  DVM  from  the  measurement  cir¬ 
cuit.  Isolation  is  accomplished  by  insertion  of  a  unity-gain  instrumentation  ampliher 
at  the  voltmeter  input  as  indicated  in  Figure  39.  Although  resulting  measurements 
conhrm  improved  response  characteristics,  there  remains  a  minor  apparent  increase 
in  measured  resistivity  at  high  frequencies  attributed  to  underestimation  of  the  test 
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current.  As  required,  calibrated  compensation  adjustments  can  be  applied  to  correct 
for  related  minor  roll-off  of  the  current  meter’s  frequency  response. 

As  regards  four-electrode  results,  it  is  observed  that  similar  high  frequency  effects  are 
encountered  for  GVI  measurements.  However,  again,  isolation  of  the  DVM  from 
impedance  characteristics  of  the  potential  measurement  circuit  (including  sample 
impedance),  largely  eliminates  related  measurement  error.  More  signihcantly,  four- 
electrode  LCR  measurements  illustrate  the  severely  restricted  function  of  these  in¬ 
struments  in  four-electrode  conhguration.  In  particular,  it  is  noted  that  the  HP4192A 
impedance  analyzer  yields  adequate  sample  characterization  only  over  a  limited  fre¬ 
quency  range  between  approximately  100  Hz  and  4  kHz.  Measurement  error  is  unac¬ 
ceptable  for  both  lower  and  higher  frequencies. 

Foregoing  results  for  untreated  Vancouver  tap  water  illustrate  both  the  limitations  of 
standard  LCR  meters  and  the  viability  of  the  considerably  lower  cost  GVI  approach. 
Although  the  latter  method  yields  only  scalar  measurement  of  apparent  electrical 
parameters  (versus  phase- sensitive  vector  impedance/admittance  provided  by  LCR 
meters),  the  advantage  of  four-electrode  measurement  capability  is  signihcant.  It  is 
notable  that  only  minor,  if  any,  electrode  polarization  effects  are  observed  for  tap 
water  samples  and  this  is  likely  attributable  to  a  relatively  high  ion  mobility  in  prox¬ 
imity  to  electrodes.  There  is  some  suggestion  of  a  minor  polarization  effect  in  the 
case  of  the  75  ml  hxture  where  it  is  observed  that  four-electrode  measurements  yield 
a  marginally  lower  resistivity  than  corresponding  two-electrode  measurements.  How¬ 
ever,  the  opposite  deviation  is  observed  in  the  case  of  the  250  ml  hxture,  suggesting 
that  minor  offsets  between  2-electrode  and  4-electrode  measurements  are  of  question¬ 
able  signihcance.  In  contrast,  as  we  shall  see,  electrode  polarization  effects  can  be 
clear  and  considerable  for  a  soil-hlled  hxture  and,  particularly  at  lower  soil  moisture 
levels. 

For  example.  Figure  41  compares  the  previous  result  for  Vancouver  tap  water  with  a 
local  loam  soil  having  similar  conductivity,  but  also  displaying  substantial  frequency 
dispersion.  The  gravimetric  moisture  content  w  =  M^/Mg  is  approximately  15%, 
where  and  Mg  denote  mass  of  water  and  mass  of  the  lab-dry  solid  fraction,  re¬ 
spectively.  In  contrast  with  results  for  Vancouver  tap  water  samples,  measurements 
on  the  loam  soil  display  clear  evidence  of  both  intrinsic  dispersion  and  electrode 
polarization.  In  particular,  two-electrode  resistivity  measurements  clearly  exceed  cor¬ 
responding  four-electrode  measurements,  particularly  at  lower  frequencies,  and  imply 
an  apparently  higher  level  of  overall  frequency  dispersion.  In  general,  these  features 
are  the  typical  signatures  of  electrode  polarization. 

As  discussed  in  Section  3.3.7,  a  common  measure  of  soil  electrical  polarization  and 
related  dispersion  is  the  electrical  chargeability  dehned  by 

^  ^  Pdc  -  Poo  ^ 

Pdc  ^oo 
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where  pdc  and  aac  denote  low-frequency  or  static  electrical  parameters  and  poo  and 
(Too  represent  corresponding  high-frequency  values.  For  example,  for  the  loam  soil  in 
Figure  41,  Equation  (137)  yields  a  chargeability  of  roughly  m  =  0.16  for  4-electrode 
GVI  measurements  compared  with  approximately  m  =  0.20  for  2-electrode  LCR  data. 

Frequency  dependent  resistivity /conductivity  for  three  DRDC  reference  soils  are  dis¬ 
played  in  Figures  42-44  on  the  same  scale  as  Figures  40  and  41.  In  view  of  evident 
limitations,  four-electrode  LCR  measurements  are  not  considered.  Results  for  the 
DRDC  Australia-A  soil  in  Figure  42  illustrate  that  in  general,  as  conductivity  in¬ 
creases,  the  absolute  level  of  electrode  polarization  effects  diminish.  Thus,  given  the 
direct  relation  between  moisture  and  conductivity  for  a  given  soil,  it  follows  that  elec¬ 
trode  polarization  phenomena  are  observed  to  increase  as  moisture  content  decreases 
for  a  particular  soil  (and  potentially  above  some  threshold  moisture  level). 

The  Australia-A  soil  also  displays  relatively  limited  absolute  dispersion.  However, 
due  to  normalization  by  absolute  conductivity /resistivity,  it  is  clear  that  electrical 
chargeability  can  remain  signihcant.  In  particular,  although  only  subtle  dispersion  is 
observable  on  the  present  scale,  chargeabilities  in  the  range  m  =  0.11  —  0.17  are  com¬ 
parable  to  those  for  the  loam  soil  in  Figure  41  with  comparatively  evident  frequency 
dependence. 

The  DRDC  Cambodia-1  soil  in  Figure  43  is  considerably  more  resistive  for  a  given 
moisture  content  than  the  previous  Australia-A  soil.  Moreover,  consistent  with  prior 
observations,  electrode  polarization  effects  are  more  substantial  and  increase  consid¬ 
erably  with  increasing  resistivity  as  moisture  decreases.  In  addition,  comparison  of  re¬ 
sults  for  Australia-A  and  Cambodia-1  soils  for  a  consistent  moisture  level  of  w  =  15% 
illustrates  that  the  extent  of  electrode  polarization  is  influenced  considerably  by  soil 
resistivity  and  related  characteristics  (ionic  content,  texture,  clay  content,  etc.)  as 
well  as  by  moisture  content. 

Results  for  the  Cambodia- 1  soil  again  demonstrate  that  greater  absolute  dispersion 
does  not  translate  directly  to  higher  chargeability.  Although  absolute  dispersion  is 
substantially  greater  than  observed  for  the  Australia-A  soil,  corresponding  charge- 
abilities  for  four-terminal  GVI  measurements  are  considerably  lower. 

Finally,  results  for  the  DRDC  (UNMAC)  Bosnia  soil  (Figure  44)  reveal  the  relatively 
conductive  nature  of  the  soil  and  further  reinforce  previously  observed  relations  be¬ 
tween  electrical  resistivity,  moisture  level,  electrode  polarization,  absolute  dispersion 
and  chargeability.  In  particular,  it  noted  that  for  a  moisture  level  of  approximately 
16%,  there  is  practically  no  visible  indication  of  electrode  polarization  on  the  present 
scale  and  no  signihcant  chargeability.  It  is  also  notable  for  the  Bosnia  soil  and,  as 
a  general  rule,  that  degraded  measurement  accuracy  for  GVI  measurements  at  high- 
frequency  is  less  apparent  with  increasing  soil  conductivity.  In  fact,  as  predictable  on 
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the  basis  of  related  prior  discussion,  the  functionality  of  LCR  and  impedance  analyzers 
in  four-electrode  conhguration  is  also  much  improved  for  conductive  soil  samples. 

Dependent  on  actual  hxture  volume  and  the  relative  scale  of  soil  heterogeneity,  it 
may  be  advisable  (depending  on  bulk  soil  texture)  to  exclude  coarser-grained  frac¬ 
tions  from  measurement.  In  general,  as  reported  in  Cross  (2000),  the  nature  of  soil 
electrical  properties  is  largely  determined  by  hner-grained  soil  fractions  due  to  compar¬ 
atively  enormous  surface  area  and  related  moisture  and  pH-dependent  ionic  exchange 
capacity /mobility.  Consequently,  for  predominantly  hne-textured  soils,  the  influence 
of  coarser-grained  gravel  components  is  ordinarily  minimal.  In  connection  with  test¬ 
ing  of  DRDC  reference  soils,  as  reported  above,  minor  gravel  fractions  were  excluded 
via  dry  sieving  at  7(^8  (6.73  mm)  mesh  (consistent  with  subsequent  magnetic  analysis 
-  Section  4.5.1). 

Because  soil  electrical  conductivity  and  related  frequency  variation  are  strongly  de¬ 
pendent  on  moisture  level,  soil  samples  should  be  well  sealed  and  related  laboratory 
characterization  should  be  carried  out  as  soon  after  sampling  as  possible.  In  addition 
to  assessing  as- received  electrical  characteristics,  the  range  and  variation  of  moisture 
dependent  electrical  characteristics  can  be  established  by  subsequently  saturating  the 
sample  with  distilled  (deionized)  water  and  acquiring  a  series  of  measurements  as  soil 
moisture  is  reduced,  initially  via  air  dying  and  ultimately  via  incremental  oven  treat¬ 
ment.  Note,  where  samples  are  allowed  to  dry  within  the  measurement  hxture,  due 
consideration  should  be  given  to  potential  inhuence  of  differential  moisture  loss  from 
the  open  surface  of  the  contained  sample.  The  advantage,  however,  is  that  soil  fabric 
remains  constant.  Alternatively,  samples  can  be  removed  from  the  hxture  between 
successive  measurements  to  ensure  uniform  drying.  Whatever  the  method,  related 
procedure  should  be  described  and  sample  moisture  level  should  be  determined  and 
recorded  in  connection  with  all  associated  measurements. 

Where  samples  are  allowed  to  air-dry  within  the  sample  hxture,  an  additional  consid¬ 
eration  is  the  potential  inhuence  of  increasing  sample-electrode  contact  impedances 
due  to  sample  contraction  and  related  physical  separation  between  electrodes  and 
sample  soil.  While  related  inhuence  at  current  electrodes  is  relatively  insignihcant 
and  can  be  overcome  by  adding  a  small  amount  of  distilled  water  at  the  electrode¬ 
sample  contact,  corresponding  effects  (particularly  signihcant  capacitive  impedance) 
at  potential  electrodes  are  relatively  difhcult  to  address  and  represent  a  limiting  source 
of  measurement  error  at  low  moisture  levels. 

A  hnal  procedural  consideration  is  the  way  in  which  sample  soil  is  placed  into  the  soil 
box  or  other  hxture  and  related  porosity / permeability  inhuence  on  measured  soil  elec¬ 
trical  properties.  As  regards  measurement  under  as-sampled  moisture  condition,  it  is 
generally  recommended  that  sample  soil  be  gradually  transfered  to  the  hxture  with  ap¬ 
plication  of  limited  shaking,  mild  mechanical  shocks  (i.e.  lightly  impacting  hxture  on 
bench  top)  and  limited  direct  pressure  to  achieve  moderate  and  uniform  compaction. 
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In  particular,  sample  soil  should  not  be  compacted  in  layers  and,  in  general,  stratifi¬ 
cation  should  be  avoided.  Where  as-sampled  conditions  are  not  of  particular  interest 
and,  generally,  in  connection  with  characterization  of  moisture-dependent  electrical 
properties,  sample  soil  should  be  transfered  to  the  hxture  in  lab-saturated  (soil  paste) 
condition  for  optimum  uniformity  of  fabric.  Despite  optimum  procedure,  however, 
it  should  be  appreciated  that  lab-measured  soil  electrical  properties  can  differ  sub¬ 
stantially  from  in  situ  characteristics  and,  wherever  possible,  in  situ  measurement  is 
preferred.  In  contrast,  soil  magnetic  properties  (excluding  remanent  magnetization) 
depend  primarily  on  soil  mineralogy  and  temperature  and  are  largely  unaffected  by 
soil  fabric  and/or  related  moisture  level/ distribution.  Consequently,  as  discussed  in 
Section  4.5,  lab-measured  soil  magnetic  susceptibility  is  adequately  corrected  for  rela¬ 
tive  compaction  via  comparatively  straight  forward  density  normalization  (assuming 
reasonable  constraints  on  in  situ  bulk  density). 

4.4.2  Laboratory  Instrumentation: 

Sample  Fixture: 

There  is  clearly  considerable  latitude  on  specifying  an  appropriate  sample  hxture. 
However,  as  described  in  the  previous  section,  it  is  recommended  that  the  hxture 
be  designed  to  facilitate  four-electrode  measurements  for  purpose  of  avoiding  the 
inhuence  of  electrode  polarization.  In  addition,  it  is  noted  that  potential  electrodes 
should  be  sufhciently  removed  from  current  electrodes  to  avoid  the  region  of  polarized 
charge  accumulation  and  should  be  separated  by  a  sufhcient  interval  to  ensure  a 
sizable  potential  gradient.  As  a  rule  of  thumb,  the  ratio  of  end-plate  electrode  area 
A  to  related  separation  £  should  satisfy  A/£  <1.0  and  should  not  exceed  A/£  =  10. 

Good  quality  connectors  should  be  employed  to  minimize  related  contact  impedances 
and  electrodes  should  be  fabricated  from  brass,  stainless  steel  or  other  corrosion  re¬ 
sistant  materials.  It  is  noted  that  non-polarizable  electrodes  (e.g.  copper/copper 
sulphate)  could  also  be  considered  to  avoid  possible  polarization  effects  at  interior 
potential  electrodes.  Although  it  is  anticipated  that  related  inhuence  is  minimal,  this 
assumption  has  not  been  investigated  or  conhrmed.  Again,  Figure  34  provides  ap¬ 
proximate  specihcation  for  soil-box  hxtures  employed  here.  Additional  guidance  is 
provided  in  ASTM  G57-06  and  G187-05. 

Impedance  Analyzers  and  Multi-Frequency  LCR  Meters: 

There  is  a  wide  range  of  commercially  available  impedance  analyzers  and  multi¬ 
frequency  LGR  meters.  Figure  45  pictures  two  precision  LGR  meters  utilized  for 
the  current  study.  For  purpose  of  measuring  frequency-dependent  soil  electrical  con¬ 
ductivity/resistivity  as  described  in  the  previous  section,  we  recommend  the  following 
general  specihcations. 
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First,  the  frequency  range  should  extend  from  10  Hz  to  100  kHz,  with  discrete  test 
frequencies  of  100  Hz,  1  kHz,  10  kHz  as  a  minimum  requirement.  Specihed  test 
signal  level  should  be  >1.0  Y rms  and  preferably  >5.0  Y rms  into  50  fl.  Typically, 
the  source  voltage  level  is  regulated  and  selectable,  with  corresponding  test  current 
dependent  on  sample  load.  Some  instruments  also  permit  current  regulation  with 
compliance  to  a  maximum  specihed  source  voltage.  In  general,  high  quality,  precision 
instruments  are  capable  of  accurate  impedance  measurements  at  test  current  levels 
as  low  as  1.0  //A.  Recommended  measurement  ranges  are  Conductance  -  G  (|y|,6*): 
<1.0  //S  to  >100.0  mS  (parallel  circuit  mode)  and  Resistance  -  R  {\Z\,  9):  <10.0  fl  to 
>1.0  MQ  (series  circuit  mode).  As  regards  measurement  accuracy,  we  recommend  a 
basic  accuracy  of  1%.  However,  actual  performance  is  a  function  of  load  impedance 
and  instrument-specihc  characteristics  as  functions  of  frequency  should  be  reviewed. 
Display  resolution  should  be  at  least  4.5  digits  and  multi-fold  measurement  averaging 
capability  is  a  signihcant  asset. 

Finally,  as  a  practical  matter,  a  particularly  signihcant  performance  parameter  is  the 
ehective  input  impedance,  which  must  be  sufficiently  high  to  avoid  loading  ehects  in 
connection  with  highly  resistive  (dry)  soils.  As  a  minimum  requirement,  we  recom¬ 
mend  Zq  >10  Mfl.  As  previously  discussed,  however,  a  major  limitation  of  investi¬ 
gated  LCR  meters  is  poor  performance  in  four-electrode  conhguration  (for  resistive 
soils)  due  to  ehectively  high  contact  resistances  between  potential-current  electrode 
pairs.  On  the  basis  of  limited  product  research  and  related  inquiries,  it  is  our  impres¬ 
sion  that  the  vast  majority  of  multi-frequency  LCR  meters  and  impedance  analyzers 
are  based  on  the  same  or  similar  auto-balancing  bridge  designs  and,  consequently,  it 
is  presumed  that  similar  restrictions  apply.  Specihc  design  characteristics  and  related 
performance  specihcations  should  be  conhrmed  with  the  instrument  manufacturer. 

GVI  Instrumentation: 

Basic  instrumentation/apparatus  for  GVI  measurements  reported  here  is  displayed 
in  Figure  46.  The  selection  of  commercially  available  signal  generators  and  digital 
multimeters  (DMM)  is  yet  greater  than  for  impedance  analyzers  and  LCR  meters. 
However,  there  is  also  far  greater  latitude  in  related  performance  characteristics  and, 
consequently,  greater  care  is  required  in  selecting  appropriate  instrumentation. 

Signal  Generator:  In  practice,  given  relatively  modest  requirements  on  frequency 
range,  nominally  10  Hz  -  100  kHz,  and  noting  that  signal  stability  effects  are  elimi¬ 
nated  by  the  ratiometric  nature  of  GVI  measurements,  the  primary  requirement  on 
the  signal  generator  is  that  it  supply  a  sufficient  signal  level.  In  particular,  we  recom¬ 
mend  >10.0  Vpp  into  50  n  as  a  minimum  specihcation.  Actual  test  signal  level  and 
related  test  current  are  dependent  on  sample  load. 

Voltage/Current  Meters:  Compared  with  modest  requirements  on  signal  source  char¬ 
acteristics,  constraints  on  DMM  performance  for  voltage  and  current  measurement 


48 


are  substantially  more  restrictive  and  it  is  generally  recommended  that  a  matched 
pair  of  instruments  be  employed.  It  is  also  recommended  that  DMMs  be  battery 
powered  to  avoid  ground  loop  currents  and  optimum  isolation. 

Meters  should  read  true  RMS  voltage/current  over  a  minimum  frequency  range  of  10 
Hz  to  100  kHz  and  with  a  response  time  of  <1.0  s.  Note  that  the  foregoing  bandwidth 
requirement  is  particularly  demanding  for  current  and  only  higher  performing  meters 
will  meet  this  specihcation^*.  Recommended  measurement  ranges  are  Voltage  -  V: 
<0.1  mV  to  >10.0  V  and  Current  -  I:  <0.1  /rA  to  >10.0  mA,  with  minimum  4.5  digit 
(12-14  bit)  resolution. 

Again,  as  for  impedance  instrumentation,  sufficiently  high  input  impedance  is  key 
to  avoiding  loading  effects,  particularly  at  higher  frequencies.  Assuming  a  high 
impedance  isolation  amplifier  is  employed  at  front-end  of  volt  meter  (as  pre  previ¬ 
ous  description),  an  adequate  input  impedance  specification  is  >10  Mfl/<  100  pF. 
Alternatively,  in  absence  of  external  isolation,  an  minimum  input  impedance  of  >100 
MQ/<  100  pF  is  recommended. 

Isolation  Amplifier:  In  general,  the  isolation  amplifier  is  a  unity-gain  instrumentation 
amplifier  having  high  input  impedance  (>100  Mfl,  <20  pF),  broad  band-width  (DC- 
100  kHz),  linear  phase  characteristics  and  high  common-mode  rejection  (>100  dB, 
<1  kHz). 

4.4.3  Field  Methods  and  Procedures: 

Four-Electrode  Resistivity: 

In  principal,  methods  described  in  Section  4.4.1  can  be  also  be  employed  for  in  situ 
measurement  of  soil  electrical  conductivity /resistivity  using  a  quadrapole  electrode 
array  as  per  Figure  47.  In  fact,  the  scale  and  very  near-surface  focus  of  landmine 
detection  are  uniquely  conducive  to  the  objective. 

Again,  within  a  uniform,  continuous  and  isotropic  medium,  the  electric  held  E  and 
the  related  current  density  J  are  related  via  Ohm’s  law 

J  =  (tE,  (140) 

where  a,  denotes  the  electrical  conductivity.  Furthermore,  noting  that  the  electric 
held  is  the  gradient  of  a  related  scalar  potential  V, 

E  =  -VV,  (141) 

1 

Although  use  of  an  external,  current-sampling  shunt  resistor  is  generally  discouraged,  where 
employed,  it  is  imperative  that  the  resistor  have  a  sufficiently  low  temperature  coefficient  (<20 
ppm/degree  C)  and  should  handle  no  more  than  10%  of  rated  power.  It  should  also  be  appreciated 
that  related  current  bandwidth  is  inevitably  reduced  (due  to  stray  capacitance  effects)  in  comparison 
with  related  voltage  specification. 
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we  have 


J  =  -(xVy. 


(142) 


Assuming  conservation  of  charge  and  uniform  conductivity,  the  potential  held  arising 
for  a  point  (current  source)  electrode  delivering  steady-state  current  /  follows  as 


27rcj  r  ’ 

where  r  denotes  radial  distance  from  the  point  electrode.  As  illustrated  in  the  up¬ 
per  panel  of  Figure  48,  resulting  equipotential  surfaces  are  hemispheres  everywhere 
orthogonal  to  radially  directed  current  How.  Now,  assuming  a  second  point  (current 
sink)  electrode  of  equal  strength  and  opposite  polarity  at  some  Hnite  distance,  the 
combined  potential  at  an  arbitrary  point  Pi  is 


V 


I 

27rcj 


(144) 


where  ri  and  r2  denote  radial  distances  between  respective  current  electrodes,  Ii 
and  I2,  and  the  point  Pi.  Again,  the  resulting  current  How  pattern  and  orthogonal 
equipotentials  are  illustrated  in  Figure  48. 

Now,  consider  the  linear  four-pole  electrode  array  depicted  in  the  lower  panel 
of  Figure  48.  Note  that  the  conHguration  comprises  two  potential  electrodes  located 
within  the  interval  between  two  current  electrodes  driven  by  a  signal  source.  In  gen¬ 
eral,  the  separation  between  any  two  electrodes  is  arbitrary.  It  follows  from  foregoing 
standard  development  that  the  potential  difference  measured  between  the  interior 
electrodes  is 


Solving  for  conductivity,  we  obtain 
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or,  alternatively,  for  the  corresponding  electrical  resistivity 
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In  general,  the  result  in  Equation  (147)  is  referred  to  as  the  apparent  electrical  resis¬ 
tivity 

Pa  =  iR,  (148) 

where  R  =  AV / 1  is  the  measured  electrical  resistance  and 


7  =  27r 


(148) 
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is  a  geometrical  factor  accounting  for  the  influence  of  the  particular  electrode  ar¬ 
rangement.  The  most  common  and  routinely  employed  configuration  for  soils  char¬ 
acterization  is  the  Wenner  array,  having  a  constant  electrode  interval  o.  Substituting 
ri  =  r4  =  o  and  r2  =  rs  =  2a  in  Equation  (148)  yields  7  =  27ro  and  the  resulting 
apparent  resistivity  pa  =  27roi?.  It  emphasized  that  only  for  a  uniform,  isotropic 
ground  do  we  have  p  =  pa- 

As  regards  the  depth  of  investigation  for  in  situ  four-electrode  resistivity  measure¬ 
ments,  initial  insight  can  be  obtained  by  considering  the  distribution  of  subsurface 
current  density.  In  particular,  assuming  a  uniform  soil,  the  horizontally  directed 
current  density  at  arbitrary  point  P  at  depth  2  within  the  vertical  plane  of  current 
electrodes  Ci  and  C2  (Figure  48)  follows  from  Equations  (142)  and  (144)  as 


Jx 


dx 

I  d  /  1 
2tt  dx  \ri 


(150) 


where  ri  =  {x"^  and  r2  =  [{I  — xY  ,  I  denoting  the  separation  between 

current  injection  electrodes  and  x  represents  an  arbitrary  offset  as  depicted  in  Figure 
48.  Now,  on  substituting  ri  and  r2  and  for  point  P  located  at  depth  2  beneath  the 
dipole  midpoint  (x  =  £/2),  we  obtain 


I  £ 


(151) 


The  predicted  x-directed  current  density  for  a  fixed  Alength  dipole  is  displayed  in 
the  upper  panel  of  Figure  49  as  a  function  of  zjl  and  with  normalization  by  the 
corresponding  surface  current  density  Jq  =  Jx{0)  =  4//7r£^.  As  anticipated,  current 
density  decreases  monotonically  with  increasing  depth  and  is  reduced  to  30%  of  max¬ 
imum  surface  value  at  zjl  ~  0.56.  What  is  more  significant,  however,  is  the  variation 
of  horizontal  current  density  for  fixed  depth  ^  as  a  function  of  variable  dipole  length  1. 
This  result  is  also  displayed  in  Figure  49  with  normalization  by  the  maximum  value, 
occurring  in  general  for  zjl  =  l/v^  0.71.  In  effect,  this  result  reveals  the  consider¬ 
able  extent  to  which  the  sampling  of  a  specified  depth  interval  is  dependent  on  array 
aperture  and  the  degree  of  related  resolution.  Having  said  this,  however,  it  is  impor¬ 
tant  to  appreciate  that  the  foregoing  result  addresses  only  source  related  sensitivity 
( “illumination” )  and  ignores  the  influence  of  depth-selective  sampling  characteristics 
associated  with  a  particular  potential  electrode  arrangements  in  relation  to  the  source 
current  dipole. 

Corresponding  sensitivity  functions  (Frechet  derivatives),  reflecting  the  combined 
depth  sampling  characteristics  of  a  given  linear  four-pole  electrode  configuration  on 
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a  uniform  conductive  half-space  are  predicted  by  the  following  relation  reported  by 
Oldenburg  (1978) 


S{r,  z)  =  s{ri,z)  -  s{r2,  z)  -  s(r3,  z)  +  s(r4,  z),  (152) 


where 

,  ^  1  (2z/n)^ 

s(rk,z)  = - 

rk  [1  +  (2z/r,pf^ 

Computed  sensitivities  are  displayed  in  the  lower  panel  of  Figure  49  for  both  the 
Wenner  array  and  an  alternative  (Schlumberger)  conhguration,  having  r4  =  r2,  as  for 
the  Wenner  array,  but  with  potential  electrode  interval  £v  =  r2—ri  =  rs—r 4  «  £.  In 
particular,  Figure  49  depicts  the  case  for  £y !£  =  0.05  compared  with  £y !£  =  0.33  for 
the  Wenner  array.  It  is  notable  that  in  addition  to  a  substantially  reduced  sensitivity 
level  (reflecting  relatively  limited  voltage  difference  for  minimally  separated  potential 
electrodes),  peak  sensitivity  is  marginally  deeper  and  less  well  focused.  In  other 
words,  the  Schlumbeger  array  provides  marginally  lower  depth  resolution  than  the 
Wenner  conhguration.  However,  the  primary  advantage  of  the  Wenner  conhguration 
is  enhanced  sensitivity.  Note  that  maximum  sensitivity  for  the  Wenner  conhguration 
occurs  at  approximately  zl£  =  0.2,  consistent  with  a  commonly  accepted  investigation 
depth  of  roughly  o/2  (see  also  Roy  and  Apparao,  1971;  Barker,  1989). 

An  alternative  perspective  on  the  effective  depth  of  investigation  considers  the  cor¬ 
responding  integrated  sensitivities  given  by  T,{r,z)  =  fg  S(r,  x)  dx.^^  Corresponding 
results  are  also  displayed  in  the  lower  panel  of  Figure  49  as  functions  of  zl£.  It  is 
noted  that  normalized  (scaled  by  associated  maximum  value)  integrated  sensitivity 
for  the  Wenner  conhguration  only  marginally  exceeds  that  for  the  Schlumberger  array 
for  given  zl£.  More  signihcantly,  it  is  observed  that  peak  sensitivity  corresponds  with 
an  integrated  sensitivity  of  approximately  30%  of  maximum,  while  70%  of  maximum 
is  reached  only  for  2  approaching  0.5£  (z  ^I.5a).  Altogether,  foregoing  results  suggest 
that,  in  theory,  the  Wenner  conhguration  provides  relatively  enhanced  sensitivity  and 
an  effective  investigation  range  approximately  equivalent  to  electrode  interval  o. 

The  resulting  integrated  sensitivity  is  given  by 

E(r,  z)  =  S(ri,  z)  -  S(r2,  z)  -  S(r3,  z)  +  S(r4,  z), 


where 


S(r,2) 


{zjrk)"^  sinh  ^{^zjrk) 
Tk  {2/rkY 


z 

(AFTlF/W+ip 


52 


Figure  50  displays  corresponding  two-dimensional  sensitivity  distributions  for  the 
same  two  arrays,  albeit  with  a  substantially  larger  ly  1^  ~  0.175  ratio  for  the  Schlum- 
berger  array  (Dahlin  and  Zhou,  2004). 

In  practice,  soils  are  rarely  uniform,  particularly  surface  soils,  and  it  is  important 
to  appreciate  related  limitations  of  foregoing  development.  Where  signihcant  depth 
dependence  or  stratihcation  is  expected,  it  is  recommended  that  small-scale  vertical 
electrical  soundings  (VES)  be  acquired  and  that  appropriate  computer-based  inver¬ 
sion  methods  be  applied  to  yield  related  characterization  of  depth-dependent  electrical 
conductivity.  Qualitative  indication  of  non-uniformity  is  obtained  by  gauging  the  re¬ 
lated  extent  and  sense  of  measured  apparent  conductivity /resistivity  with  incremental 
increase  of  electrode  interval  o.  Where  signihcant  variation  is  observed,  an  alterna¬ 
tive  to  computer-based  inversion  is  to  acquire  a  sequence  of  small-scale  hxed-interval 
measurements  at  incremental  depth  within  one  or  several  test  pits. 

Although  only  limited  testing  has  been  carried  out,  initial  results  suggest  that  small- 
scale  quadrapole  measurements  employing  a  standard  function  generator  and  matched 
digital  multimeters  is  a  viable  approach  for  in  situ  characterization  of  frequency  de¬ 
pendent  electrical  conductivity /resistivity  of  near-surface  soils.  Related  research  and 
development,  however,  remains  a  work  in  progress. 

Preliminary  results  for  two  local  soils  are  displayed  in  Figures  51.  In  particular, 
apparent  electrical  resistivity  is  displayed  for  loam  and  sand  soils  under  generally 
moist  conditions  and  following  one  or  more  minor  precipitation  events.  Measured 
apparent  electrical  resistivity  for  Wenner  electrode  intervals  o=10,  20  and  40  cm 
indicates  that  depth-dependent  resistivity  is  signihcantly  different  for  the  two  soils. 
As  anticipated  the  sand  soil  is  generally  more  resistive  than  the  loam  soil.  Moreover, 
relatively  resistive  surface  conditions  (o=10-20  cm)  for  the  sand  are  consistent  with 
higher  porosity /permeability  and  lower  moisture  retention  capacity,  permitting  the 
relatively  rapid  downward  percolation  of  precipitation.  In  contrast,  the  loam  soil  is 
substantially  more  conductive  at  surface  and  resistive  at  depth  (o=40  cm). 

More  signihcantly,  and  consistent  with  general  expectation,  the  apparent  resistivity  of 
both  soils  decreases  signihcantly  with  increasing  frequency  (low-frequency  AC  resis¬ 
tivity  (10  Hz)  is  generally  in  close  agreement  with  commutated  DC  measurements). 
Moreover,  as  anticipated,  chargeability  (Equation  139)  for  the  loam  soil  (10  Hz-I 
kHz)  is  roughly  twice  that  for  the  sand.  As  for  initial  laboratory  measurements  using 
the  same  instrumentation,  however,  abrupt  and/or  accelerated  fall-oh  at  frequencies 
above  approximately  10  kHz  is  likely  attributable  to  loading  ehects  that  have  sub¬ 
sequently  been  addressed  by  providing  enhanced  isolation  at  potential  measurement 
inputs.  In  addition,  irregular  high  frequency  variation  (increased  apparent  resistivity 
in  the  range  of  5-20  kHz)  observed  for  the  sand  soil  with  o  =  0.4  m  may  rehect  similar 
instrumental  ehects. 
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Unfortunately,  due  to  the  preliminary  nature  of  these  experiments,  actual  soil  mois¬ 
ture  levels  and  precipitation  history  were  not  monitored.  In  general,  further  testing 
is  required  to:  1.  verify  the  effectiveness  of  enhanced  isolation  (as  per  laboratory 
measurements);  2.  investigate  potential  influence  of  limited  signal  levels  for  larger 
electrode  intervals;  and  3.  to  assess  performance  over  a  wider  range  of  soil  conditions, 
including  more  conductive  soils.  It  is  also  planned  to  investigate  the  application  of 
basic  ID  inversion  for  frequency /depth-dependent  resistivity /conductivity  structure. 


FDEM  Conductivity: 

Alternatively,  held  measurements  of  electrical  conductivity  can  be  acquired  via  a  range 
of  dual-coil  frequency  domain  electromagnetic  (FDEM)  devices  that  are  primarily 
designed  for  engineering  and  agricultural  investigations.  A  sample  of  representative 
instruments  is  pictured  in  Figure  52,  with  related  brief  description  in  Section  4.4.4. 

In  general,  FDEM  conductivity  meters  operate  in  the  so-called  low  induction  number 
(LIN)  regime.  In  particular,  with  reference  to  previous  related  development  in  Section 
3.1,  the  LIN  regime  is  dehned  as  77  1.0  and,  in  practice,  instruments  are  designed 

to  operate  within  the  range  77  <  0.5  as  outlined  in  Figure  53,  where  the  quadrature 
response  is  approximately  a  linear  function  of  conductivity.  More  generally,  as  em¬ 
phasized  in  Section  3.1,  the  response  characteristics  of  a  given  dual-coil  conhguration 
are  a  complicated  function  of  coil  orientation  and  frequency  (e.g.  Tabbagh,  1986). 
However,  commercially  available  systems  are  typically  based  on  a  hxed-offset  coplanar 
coil  system  that  is  operable  in  either  horizontal  coplanar  (HC  or  “vertical  dipole”)  or 
vertical  coplanar  (VC  or  “horizontal  dipole”)  via  90  degree  axial  rotation. 

As  per  Equation  (29)  in  Section  3.1,  the  induction  number  is  given  by 

rj  =  as  =  -^  =  (7rf/j,ay/^s,  (29) 


where  the  so-called  response  parameter  a  is  the  real  part  of  the  low-frequency  (quasi¬ 
static)  propagation  constant  7  ~  (727r///cT)^/^  (see  Chapter  2,  Equation  12)  and 
effectively  the  reciprocal  of  the  corresponding  skin  depth  dg  =  (tt/t/ct)"^/^.  s  denotes 
the  separation  between  transmitter  and  receiver  coils. 

Under  the  LIN  assumption  s  ^  dg  for  coplanar  coils,  the  quadrature  component  of 
the  secondary  to  primary  magnetic  held  ratio  Hg/Hp  (see  Footnote  8  -  Equation  37.5) 
at  the  receiver  coil  can  be  approximated  (McNeill,  1980,  Kaufman  and  Keller,  1983) 
as 

Hg  \  TTfiias^ 

Hp)q""  2 

with  result  that  an  appropriately  designed  FDEM  instrument  is  capable  of  reading 
electrical  conductivity  directly  via  appropriate  calibration  and  scaling. 
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In  complete  analogy  with  sensitivity  fnnctions  described  for  qnadrapole  electrode  ar¬ 
rays,  corresponding  relations  have  been  derived  (McNeill,  1980;  Kanfman  and  Keller, 
1983)  for  the  s-weighted  depth  sensitivity  of  horizontal  (HC)  and  vertical  (VC)  copla- 
nar  systems.  Resnlts  are  respectively 


Shc-q{Z) 


4.Z 

(4^2  +  1)3/2 


(155) 


and 

AZ 

Svc-q{Z)  =  2-  ^^^2  ^1)1/2  (156) 

where  Z  =  z/s.  Foregoing  resnlts,  depicted  graphically  in  Fignre  54,  illnstrate 
that  while  the  horizontal  coplanar  (“vertical  dipole”)  coil  conhgnration  yields  depth- 
focnsed  sampling  characteristics  similar  to  those  for  qnadrapole  electrical  arrays,  the 
vertical  coplanar  ( “horizontal  dipole” )  conhgnration  has  maximum  sensitivity  at  snr- 
face,  decreasing  monotonically  with  increasing  depth.  In  effect,  measnred  condnctivity 
represents  a  corresponding  depth- weighted  average  value.  In  the  lower  panel  of  Figure 
54  we  display  corresponding  integrated  sensitivities  S(Z)  =  S(x)dx  and  it  is 
notable  that  the  vertical  coplanar  coil  arrangement  attains  70%  of  integrated  response 
within  approximately  ^  =  0.75s  compared  with  roughly  twice  that  ^  =  1.5s  for  the 
horizontal  coplanar  conhgnration. 

Given  the  scale  of  smaller  FDEM  instruments  (s  ^1.0  m),  foregoing  results  suggest 
that  these  devices  are  not  ideally  suited  for  the  limited  depth  of  interest  in  connection 
with  metal  detector  performance  evaluation.  Although  depth  sensitivity  character¬ 
istics  can  be  selectively  focused  to  some  degree  by  appropriate  adjustment  of  sensor 
height  above  grade,  the  extent  of  potential  measurement  bias  due  to  deeper  soils 
should  be  given  due  consideration. 

As  described  above  for  four-electrode  electrical  conductivity  measurements,  vertical 
electromagnetic  soundings  can  be  acquired  by  recording  a  series  of  measurements  with 
the  instrument  raised  incrementally  above  grade.  Computer-based  inversion  methods 
are  applied  to  resulting  measurements  to  derive  an  optimized  estimate  of  actual  depth- 
dependent  electrical  conductivity.  Similarly,  a  rapid  semi- quantitative  indication  of 
depth-variability  is  obtained  by  noting  the  variation  in  instrument  readings  as  the 

Corresponding  integrated  sensitivities  follow  from  Equations  (155)  and  (156)  as 

Shc-q(Z)  =  1  -  1/(4Z^  +  1)^/2  Evc-q(Z)  =  1  +  2Z-  (4Z^  +  1)'/% 
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device  is  raised  from  grade  to  a  height  of  approximately  1.5s.  For  a  uniform  soil, 
displayed  conductivity  decrease  monotonically  with  increasing  height  in  accordance 
with  integrated  sensitivity  relations  in  the  lower  panel  of  Figure  54.  In  particular, 
s- weighted  height  dependent  response  above  a  uniform  ground  is  given  by  R{H)  = 

S{x)  dx  ,  where  H  =  h/s  and  h  is  height  above  grade.  Corresponding  response 
characteristics  are  depicted  in  Figure  54  (see  McNiell  and  Bosnar,  1999). 

A  similar  semi-quantitative  assessment  of  spatial  heterogeneity  is  obtained  by  com¬ 
paring  measurements  at  grade  with  the  instrument  in  two  orthogonal  orientations. 
For  a  uniform  or  vertically  stratihed  soil,  the  two  readings  are  identical  and,  in  gen¬ 
eral,  the  degree  of  departure  between  the  two  readings  is  directly  related  to  the  extent 
of  lateral  heterogeneity. 

Finally,  it  is  our  experience  that  FDEM  conductivity  meters  are  not  particularly  ac¬ 
curate  or  stable  in  an  absolute  sense,  with  related  error  due  principally  to  imperfect 
nulling  of  the  the  primary  held  and  related  thermal  drift.  Consequently,  where  in¬ 
struments  of  this  variety  are  to  be  employed,  the  device  should  be  operated  in  surface 
sensitive  conhguration  (i.e.  appropriate  coil  orientation  /  sensor  height)  and  manufac¬ 
turer  specihed  calibration  adjustments  should  be  carried  out  regularly  (for  optimum 
accuracy  in  advance  of  individual  measurements),  to  control  zero  drift.  In  general, 
while  FDEM  instruments  are  clearly  advantageous  for  mapping  spatial  variability  in 
soil  electrical  conductivity,  it  is  our  experience  that  four-electrode  measurements  are 
preferable  for  accurate  absolute  characterization,  including  for  provision  of  calibration 
and  constraints  on  FDEM  mapping. 

4.4.4  Field  Instrumentation: 

Four-Electrode  Resistivity: 

Instrumentation  requirements  for  in  situ  four-electrode  conductivity  (resistivity)  mea¬ 
surement  are  largely  the  same  as  described  in  Section  4.4.2  for  equivalent  laboratory 
measurements.  As  noted,  related  research  and  development  remains  ongoing.  How¬ 
ever,  it  is  clear  that  principal  restrictions  are  limited  source/transmitter  voltage  and 
the  input  impedance  of  potential  measurement  instrumentation.  As  a  practical  mat¬ 
ter,  even  reasonably  high-voltage  (>30  Vpp)AC  waveform  generators  are  a  rare  com¬ 
modity  and,  although  there  is  no  shortage  of  audio  frequency  amplihers,  amplihers 
having  an  upper  band-limit  beyond  20  kHz  are  also  relatively  scarce.  In  general, 
electrode  contact  resistance  is  the  chief  obstacle,  with  considerable  source  voltage 
required  to  establish  a  sufficient  test  current  level  (nominally  >0.1  mA).  Similarly, 
while  it  is  not  unusual  for  voltmeters  (and,  particularly,  “electrometers”)  to  provide 
a  DC  input  resistance  in  excess  of  100  Mfl  (and  as  high  as  100  GO),  it  is  rare  for  AC 

2 1 

Resulting  response  functions  follow  from  Equations  (155)  and  (156)  as 

Rhc-q{H)  =  I/(4H2  +  1)1/2  Ryc_Q{H)  =  (4H2  +  I)i/2  -  2H. 
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input  impedance  to  exceed  10  Mil,  typically  shunted  by  <100  pF,  due  to  limitations 
on  effective  bandwidth  (i.e.  time  constant).  For  this  reason,  an  isolation  ampliher 
(see  Section  4.4.2),  having  input  impedance  in  excess  of  100  Mfl,  is  a  necessity  to 
avoid  related  loading  effects. 

As  regards  the  electrode  array,  it  is  recommended  that  copper  or  stainless  electrodes 
be  used  with  a  minimum  diameter  >5.0  mm  and,  as  a  rule  of  thumb,  electrodes  should 
not  be  inserted  to  a  depth  exceeding  approximately  1/5  of  the  electrode  interval  o. 
Although  it  may  be  advantageous  to  employ  porous,  non-polarizable  electrodes  for 
potential  detection,  preliminary  comparison  between  stainless  and  Cu-CuS04  elec¬ 
trodes  has  revealed  no  indication  of  significant  electrode  polarization  in  tested  soils. 
Finally,  it  is  recommended  that  the  length  of  all  leads  be  kept  to  minimum  and  that 
coaxial  cabling  be  employed  to  as  great  an  extent  as  possible. 

Again,  as  for  laboratory  measurements,  it  also  possible  to  use  a  precision  impedance 
analyzer  or  LCR  meter  to  the  extent  that  it  is  sufficiently  portable.  However,  as 
discussed  in  Sections  4.4.1  and  4.4.2,  the  common  auto-balancing  bridge  design  is 
not  designed  for  four-electrode  measurements  and  related  measurement  error  can  be 
substantial,  particularly  at  high-frequency  and  for  resistive  soils.  In  general,  regardless 
of  instrumentation,  it  is  recommended  that  characterization  of  frequency  dependent 
soil  electrical  conductivity  be  constrained  by  coincident  DC  measurements  with  a 
integrated  and  purpose-built  resistivity  meter. 

FDEM  Conductivity  Meters: 

A  range  of  FDEM  conductivity  meters  are  commercially  available.  A  sample  of  com¬ 
mon  instruments  is  pictured  in  Figure  52,  together  with  related  instrument  design 
parameters.  As  noted,  the  vast  majority  of  commercially  produced  instruments  em¬ 
ploy  a  coplanar  coil  arrangement.  However,  we  are  aware  of  at  least  one  small-scale 
instrument  (DUALEM-IS:  /=9  kHz,  s=1.0  m)  that  combines  a  perpendicular  coil 
arrangement  with  vertical  coplanar  coils,  to  yield  complementary  surface  sensitive 
and  depth  focused  sampling  characteristics.  The  Geonics  EM38-DD  and  EM38-MK2 
provide  similar  dual  sensitivity  characteristics  by  conhguring  two  separate  coplanar 
instruments  in  orthogonal  orientation  about  axis  (the  latter  instrument  is  a  recent 
integrated  system  with  added  flexibility  of  dual  coil  separations  of  1.0  m  and  0.5  m). 
Note  (Figure  52)  that  in  addition  to  single  frequency  devices,  instruments  provid¬ 
ing  multiple  measurement  frequencies  are  also  available  (e.g.  Geophex  GEM-2  and 
GSSI  EMP-400),  allowing  assessment  of  related  frequency  dependence  and  associated 
viscous  influence  on  TDEM  sensors.  Fortunately,  under  LIN  conditions,  depth  sensi¬ 
tivities  given  by  Equations  (155-156)  are  not  significantly  dependent  on  frequency,  so 
that  variable  frequency  measurements  with  a  given  device  provides  direct  indication  of 
electrical  dispersion  (assuming  accurate  calibration  adjustments  and  linear  response 
characteristics). 
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4.5  Magnetic  Parameters  -  Measurement: 

Although  there  are  no  formalized  standards  for  measurement  of  soil  magnetic  prop¬ 
erties,  laboratory  measurement  of  magnetic  susceptibility  and  related  frequency  de¬ 
pendence  via  the  Bartington  MS2B  dual  frequency  system  is  a  de  facto  standard.  As 
described  in  Section  3.3,  magnetic  viscosity  and  related  frequency  dispersion  of  the 
magnetic  susceptibility  have  a  uniquely  significant  influence  on  the  performance  of 
pulse  induction  metal  detectors  and  consequently,  development  and  assessment  of  in¬ 
strumentation  for  characterization  of  frequency-dependent  magnetic  parameters  has 
been  a  principal  focus.  We  begin  by  describing  impedance  based  measurements  as 
previously  developed  in  Section  4.3  and  compare  results  with  corresponding  measure¬ 
ments  via  the  Bartington  MS2B  and  the  recently  available  SM-100  multi-frequency 
system  from  ZH  Instruments.  We  subsequently  provide  a  brief  description  of  portable 
systems  for  in  situ  susceptibility  measurement. 


4.5.1  Laboratory  Methods  and  Procedures: 

At  the  outset  of  work  reported  herein,  characterization  of  frequency-dependent  sus¬ 
ceptibility  relied  primarily  on  the  Bartington  MS2B  dual  frequency  system,  yielding 
magnetic  susceptibility  at  456  Hz  and  4.65  kHz.  Consequently,  in  connection  with 
measurement  of  frequency-dependent  electrical  conductivity  our  initial  aim  was  to 
investigate  the  practical  potential  of  multi-frequency  impedance-based  measurements 
of  complex  magnetic  susceptibility  as  developed  in  Section  4.3. 

Initially  and  in  analogy  with  the  soil-box  sample- fixture  described  in  Section  4.4.1,  a 
simple  solenoidal  sample  fixture  was  fabricated  to  permit  impedance  analysis  of  soil 
magnetic  parameters.  The  prototype  fixture  is  displayed  in  Figure  55  together  with 
pertinent  data  and  a  graph,  depicting  theory-predicted  axial  magnetic  flux  density  as 
a  function  of  position  along  the  coil’s  axis  (with  reference  to  coil  midpoint).  Axial 
flux  density  is  computed  via  the  following  approximate  relations 
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for  an  Alength  solenoid  about  the  x-axis,  centered  on  the  origin  and  comprising  N 
closely  spaced  turns  of  radius  R  (see  Purcell,  1963).  In  particular,  peak  flux  density 
at  the  center  of  the  test  fixture  is  approximately 
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yielding  B  ^  12  //T  and  corresponding  field  H  =  B/hq  ^  10  A/m  for  a  nominal  drive 
cnrrent  of  /  =  10  mA, 

A  rongh  estimate  of  corresponding  self  indnctance  follows  as 
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with  result  Lq  ^  110  /rH  for  the  air-cored  fixture. 
(Wheeler,  1982) 
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A  more  accurate  approximation 


(160), 


yields  Lq  =  98  //H,  in  reasonably  good  agreement  with  a  measured  inductance  of 
Lq  ~  102  /rH. 


Complete  frequency- dependent  complex  impedance  characteristics  of  the  fixture  (mea¬ 
sured  via  a  HP  4192A  low-frequency  impedance  analyzer)  are  displayed  in  Figure  56. 
By  design,  the  fixture’s  effective  resistance  and  inductance  are  roughly  frequency 
independent  over  the  desired  measurement  band  of  100  Hz  to  100  kHz. 


As  a  practical  matter  it  is  difficult  to  achieve  a  fiat  response  over  more  than  approxi¬ 
mately  two  or  three  decades.  Band  limits  involve  tradeoffs  between  resistance,  induc¬ 
tance  and  capacitance  of  coil  windings.  Although  the  useful  measurement  band  can 
extend  significantly  lower,  the  transition  frequency  fuL  =  Ro/‘^t^Lq  may  be  assumed 
as  a  conservative  lower  band-limit.  High-frequency  response  is  ultimately  restricted 
by  the  coil-specific  self  resonant  frequency  /g  =  l/27r(LoC'o)^/^  (approximately  3  MHz 
for  the  present  fixture). 

As  developed  in  Section  4.3,  soil  magnetic  susceptibility  is  determined  by  measuring 
the  deviation  in  the  inductance  AL  =  L  —  Lq  and/or  resistance  AR  =  R  —  Rq  oi  the 
test  fixture  on  inserting  the  sample  core.  As  pictured  in  Figure  55,  the  soil  sample 
is  contained  in  thin- walled  sample  tube,  having  outside  diameter  only  marginally 
less  than  the  inside  diameter  of  the  coil  former.  As  a  first  approximation,  in-phase 
and  quadrature  components  of  the  complex  soil  magnetic  susceptibility  are  given  by 
Equations  137  and  138  as  y'  =  AL/q^/zo  =  AL/Lq  and  y"  =  AR/2tt f^^HQ  = 
Ai?/27r/Lo,  respectively. 

More  precise  relations,  accounting  for  actual  core  radius  are 
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where  Acou  and  Acore  denote  respective  cross-sectional  areas  and  Kf  represents  an 
empirical  correction  factor  acconnting  for  shape-related  demagnetization  and  residnal 
hxtnre  effects  (Epstein,  1954).  In  practice,  given  ~  1-0  for  soils  and  assnming 
a  snfficiently  long  sample  core  (>~  3£),  demagnetization  effects  are  minimal.  The 
residnal  hxtnre  constant  Kf  is  determined  by  way  of  calibration  with  a  chemical 
componnd  having  known  magnetic  snsceptibility. 

Initial  calibration  of  the  prototype  hxtnre  coil  was  carried  nsing  reagent-grade  man¬ 
ganese  carbonate  (MnCOs),  an  antiferromagnetic  (parasitic  ferromagnetic)  componnd 
with  well  established  and  non-dispersive  mass  snsceptibility  Xp=l. 25x10“®  m^/kg  at 
20  °C  (Weast,  1971).  Resnlts  are  snmmarized  in  Fignre  57,  depicting  measnred  real¬ 
valued  mass  susceptibility  Xp  =  x'  / Pbi  where  pt  denotes  bulk  mass  density,  as  a 
function  of  test  frequency  between  100  Hz  and  100  kHz.  Inductance  measurements 
were  acquired  via  a  HP4274A  precision  LCR  meter.  In  general,  there  is  approximate 
agreement  with  the  tabulated  value  (Weast,  1971)  and  also  with  measurements  re¬ 
ported  by  Fukuma  (2000)  for  two  commonly  used  instruments  (Bartington  MS2B  and 
AGICO  Kappabridge  KLY-3).  Subsequent  independent  measurements  via  DRDC’s 
Bartington  MS2B  were  also  consistent  (xp=l. 212x10“®  m^/kg).  Accepting  the  tab¬ 
ulated  value,  a  mean  measured  mass  susceptibility  of  Xp=l. 1834x10“®  m^/kg  via 
measured  inductance  contrast  AL  implies  a  hxtnre  calibration  constant  R/=1.053. 

Subsequently,  in  connection  with  a  parallel  investigation  of  commercially  available 
instruments  for  measurement  of  soil  magnetic  susceptibility  (North,  2006),  an  addi¬ 
tional  calibration  check  was  undertaken  on  a  paramagnetic  Gd2  02  (gadolinium  oxide) 
standard  supplied  by  the  USAGE-ERDG  (North,  2005).  Again,  measured  suscepti¬ 
bility  was  on  average  in  approximate  agreement  with  results  of  an  inter-laboratory 
calibration  effort  reported  by  Sagnotti,  et  al.  (2003).  Results  are  displayed  in  Figure 
58. 

While  initial  calibration  results  establish  the  viability  of  direct  inductance  susceptibil¬ 
ity  measurements,  substantial  deviation  is  observed  about  mean  susceptibility  values. 
To  a  signihcant  extent,  observed  measurement  error  has  been  subsequently  reduced  by 
improvements  to  the  hxtnre  apparatus  and  through  development  of  computer-based 
measurement  control  software.  Details  are  provided  in  Section  4.5.2. 

Ultimately,  the  practical  limitation  of  direct  inductance  measurements  is  the  very 
considerable  dynamic  range  required  to  resolve  a  necessarily  small  deviation  in  a  rel¬ 
atively  large  quantity.  In  particular,  it  follows  directly  from  Equation  137  that  a 
basic  accuracy  of  0.01-0.1  %  for  highest  precision  impedance  measurement  devices 
imposes  an  effective  limit  on  volume  susceptibility  resolution  on  the  order  of  10“^  SI. 
In  practice,  an  effective  resolution  approaching  10“®  SI  can  be  achieved  for  real-valued 
susceptibility  via  measurement  averaging  and  with  real-time  computer-based  compen¬ 
sation  for  thermal  drift  associated  with  minor  temperature-dependent  huctuation  of 
hxtnre  geometry.  However,  because  the  quadrature  susceptibility  is  typically  an  order 
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of  magnitude  smaller  and  because  fixture  resistance  Rq  is  dramatically  more  sensi¬ 
tive  to  thermal  fluctuation,  corresponding  resolution  would  be  substantially  lower  and, 
consequently,  related  measurement  of  quadrature  susceptibility  has  not  been  pursued. 

As  we  shall  see  in  connection  with  subsequent  comparison  of  susceptibility  measure¬ 
ments  for  a  number  of  DRDC  reference  soils,  the  accuracy  of  special  purpose  suscep¬ 
tibility  meters,  employing  a  resonant  circuit  approach,  is  approximately  an  order  of 
magnitude  better. 

Resonant  frequency  susceptibility  meters,  including  the  Bartington  MS2B  and  ZH 
Instruments  SM-100  are  based  on  tuned  LC  oscillator  circuits  that  incorporate  the 
sensor  coil  as  the  inductive  circuit  element.  As  for  direct  inductance  measurements, 
magnetic  susceptibility  is  measured  via  a  two-step  measurement  cycle,  comprising  an 
initial  null  or  calibration  reading  for  the  air-cored  sensor  and  a  subsequent  reading 
with  the  sample  inserted  into  the  sensor  core.  Rather  than  measuring  the  resulting 
change  in  complex  impedance,  however,  these  instruments  detect  a  corresponding  shift 
in  the  resonant  frequency  A/*  =  /*  —  /o  ,  where  /g  =  l/27r(Lo^o)^^^  is  the  reference, 
open-coil  resonant  frequency  and  /*  =  l/27r(L(7o)^/^  is  the  corresponding  resonant 
frequency  on  insertion  of  the  sample  core.  Assuming  the  coil’s  effective  capacitance 
Co  (including  hxed  or  adjustable  tuning)  is  constant,  the  measured  shift  in  resonant 
frequency  yields  the  corresponding  inductance  contrast  AL  and  related  susceptibility 
estimate  y'  =  AL/j^hq,  where  7^  denotes  the  appropriate  coil  or  hxture  constant. 
The  advantage  is  that  frequency  (counting)  measurements  are  substantially  more 
accurate  than  complex  impedance  (voltage)  measurements,  particularly  with  limited 
dynamic  rage. 

Specihc  circuit  details  of  MS2B  and  SM-100  instruments  have  not  been  described 
by  the  manufacturers  and,  consequently,  it  is  uncertain  whether  the  tuned  circuit 
is  driven  by  a  reference  oscillator  or  whether  the  resonant  frequency  shift  is  mea¬ 
sured  directly,  as  a  deviation  from  balanced  mixed  output  or  via  detection  of  related 
phase-shift.  Precise  mode  of  frequency  selection/tuning  for  dual  or  multi-frequency 
measurement  capability  is  also  uncertain,  but  presumably  involves  manual  or  re¬ 
lay  controlled  switching  of  matched  LC  circuit  components.  Interestingly,  however, 
limited  measurements  of  frequency  dependent  susceptibility  presented  here  and  else¬ 
where,  suggest  that  dual  frequency  measurement  may,  ultimately,  provide  sufficient 
characterization  in  practice. 

Figures  59,  60  and  61  display  measured  frequency-dependent  magnetic  susceptibility 
for  three  DRDC  reference  soils  having  signihcant  but  substantially  different  magnetic 
characteristics.  In  general,  there  is  good  agreement  between  magnetic  susceptibility 
determined  via  complex  impedance  analysis  and  corresponding  measurements  with 
commercially  available  resonant  frequency  instruments.  Results  also  reveal  a  range 
of  interesting  and  potentially  signihcant  issues  regarding  the  inhuence  of  soil  texture 
and  heterogeneity  on  measured  susceptibility  and  the  extent  of  related  error. 
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To  facilitate  meaningful  comparison  of  results  for  a  giveu  soil  aud  betweeu  soils,  mass- 
specific  susceptibility  Xp  =  x/p  is  charted  to  accouut  for  variable  sample  compactiou 
aud  related  effect  of  mass  deusity  p.  lu  geueral,  it  is  recommeuded  that  mass-specihc 
susceptibility  be  the  primary  reported  quautity.  Alteruatively,  where  volume-specihc 
susceptibility  is  reported,  sample  deusity  should  also  be  specihed.  Ideally,  p  =  pb, 
where  pb  deuotes  dry  bulk  deusity.  However,  iu  geueral,  for  huite  moisture  couteut 
p  >  Pb  and  the  iuflueuce  of  soil  moisture  level  is  poteutially  siguihcaut. 

As  discussed  iu  Sectious  3.3.4  aud  3.3.6,  a  commou  measure  of  frequeucy  depeudeuce 
of  maguetic  susceptibility  is 


Xfd 


Xlf  -  Xhf 
Xlf 


X  100%. 


(163) 


However,  uulike  the  electrical  chargeability  dehued  iu  Sectious  3.3.7  aud  4.4.1,  xfd  is 
referred  to  a  specihc  measuremeut  baud,  with  xlf  aud  Xhf  deuotiug  maguetic  sus¬ 
ceptibility  at  measuremeut  frequeucies  /l  aud  fn  a  decade  apart  (i.e.  log^fn/fL)  = 
1).  Iu  particular,  Xfd  was  priucipally  established  iu  connection  with  the  Bartington 
MS2B  system  with  measurement  frequencies  =465  Hz  and  fn  =4.65  kHz.  As  dis¬ 
cussed  iu  Sectiou  3.3.6,  a  more  relevaut  iudicator  of  maguetic  dispersiou,  as  it  affects 
TDEM  seusors,  is  the  differeutial  susceptibility  Ax  =  Xlf  —  Xhf,  which  is  effec¬ 
tively  equivaleut  to  the  maguetic  viscosity  (see  Equatiou  103  aud  related  couditious 
iu  Sectiou  3.4) 

_  dx'if)  _  Xlf  -  Xhf  _  Ax  .  . 

^  dlog  f  logifn/fL)  logifn/fL) 
for  log(/H//L)  =  1-  ludeed,  it  is  the  maguetic  viscosity  that  is  respousible  for  the  sus- 
taiued  time-domaiu  respouse  that  effectively  masks  the  couductive  sigual  aud  thereby 
degrades  performauce  of  pulse  iuductiou  metal  detectors.  Cousequeutly,  appareut 
viscosity  as  iudicated  by  log  spectral  slope  u  =  dx{f)/dlog  f  is  the  ideal  practical 
predictor  of  soil  iuflueuce  ou  time-domaiu  seusors. 

Where  multi-frequeucy  characterizatiou  of  soil  maguetic  susceptibility  is  available,  aud 
where  log  spectral  slope  is  approximately  coustaut  (typically  the  case  as  discussed  iu 
Sectious  3.3.6  aud  3.4),  appareut  viscosity  is  determiued  directly  via  log-liuear  regres- 
siou  as  u  =  dx' if) /dlog  f .  Moreover,  assumiug  measuremeut  baudwidth  exceeds  a 
decade,  the  regressiou-determiued  viscosity  is  equivaleut  to  the  differeutial  susceptibil¬ 
ity  V  =  Ax  for  auy  two  measuremeut  frequeucies  separated  by  a  decade. Iu  geueral, 
for  dual  frequeucy  measuremeuts  at  arbitrary  frequeucies,  the  correspoudiug  relatiou 
is  u  =  Ax/log(/H//L)-  Note  that  for  mass-specihc  susceptibility  spectra  displayed 
iu  Figures  59-61,  related  maguetic  viscosity  is  u  =  where  Vp  =  dxp^fy /dlog  f 
represeuts  slope-predicted  viscosity  via  the  mass  susceptibility  spectrum. 

In  contrast,  it  is  noted  that  XFD  —  {^x/XLf)  X  100%  =  {f/xlf)  X  100%  depends  on 
reference  frequency  fj^  and  related  susceptibility  XLF- 
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Results  for  the  DRDC  Bosnia  soil  in  Figure  59  indicate  that  this  soil  possesses  a  con¬ 
siderable  magnetic  susceptibility  with  signihcant  frequency  dependence.  A  relatively 
non-magnetic  gravel  fraction  was  excluded  from  measurement  via  dry  separation  on 
a  ^3  mesh  (6.73  mm)  sieve.  Although  the  effective  sample  volume  of  the  solenoidal 
hxture  (approximately  100  cc)  for  direct  inductance  measurements  could  clearly  ac¬ 
commodate  larger  grains,  coarse-grained  soil  fractions  were  excluded  to  permit  con¬ 
sistent  comparison  with  measurement  systems  employing  a  relatively  limited  sample 
volume  (i.e.  10  cc)^^.  Measured  mass-specihc  susceptibility  at  1  kHz  is  in  the  range 
Xp  =4.0-4.75  xlO  ®  m^/kg  (mean  sample  mass  density  as  indicated). 

In  general,  there  is  good  agreement  between  measurements.  Five  repeat  susceptibility 
spectra  were  measured  via  the  direct  inductance  apparatus  for  a  single  sample.  Mea¬ 
surements  were  acquired  using  a  computer  controlled  ESI  2150  precision  LCR  meter 
at  ten  discrete  frequencies  between  100  Hz  and  100  kHz.  Measured  susceptibility  at 
each  frequency  represents  a  20-fold  average  of  50  ms  measurements  at  a  held  of  ap¬ 
proximately  50  A/m.  Log- linear  regression  yields  the  indicated  ht  and  observed  misht 
provides  an  indication  of  effective  measurement  precision.  As  anticipated,  measure¬ 
ment  error  is  evidently  most  signihcant  at  the  band- limits  of  100  Hz  and  100  kHz.  In 
general,  however,  results  are  in  good  agreement  with  measurements  acquired  via  the 
two  commercial  instruments  and  indicates  that  the  log-linear  decrease  in  magnetic 
susceptibility  with  increasing  frequency  extends  well  beyond  the  range  of  standard 
measurements . 

Although  susceptibility  measurements  via  MS2B  and  SM-100  instruments  are  mod¬ 
estly  higher,  the  deviation  is  likely  attributable  to  sample  heterogeneity  and  relatively 
smaller  sample  volumes.  In  addition  to  MS2B  measurements  at  465  Hz  and  4.65  kHz 
for  a  standard  10  cc  sample  volume,  corresponding  spectra  are  also  displayed  for  read¬ 
ings  at  500  Hz,  1  kHz,  2  kHz,  4kHz  and  8  kHz  as  measured  by  the  SM-100  instrument, 
using  both  an  instrument-standard  35  mm  hlm-canister  sample  holder  (approximately 
30  cc)  and  a  Bartington  standard  10  cc  pot  centered  within  a  purpose-built  35  mm 
adapter  (see  Section  4.5.2).  In  each  case,  10  susceptibility  spectra  were  acquired 
within  a  consistent  source  held  of  80  A/m.^^ 

Specification  of  a  ^3  mesh  (6.73  mm)  sieve  involves  a  subjective  trade-off  between  representa¬ 
tive  soil  characterization  and  the  need  for  standardized  and  comparable  measurements.  A  ^4  mesh 
(4.75  mm)  is  more  commonly  specified  and  defines  the  textural  transition  between  sand  and  gravel 
(United  Soil  Classification  -  USC  system).  However,  resulting  separation  is  substantially  more  re¬ 
strictive  and,  in  contrast,  a  3/8  in  (9.5  mm)  mesh  is  seemingly  too  large  in  relation  to  a  10  cc  sample 
volume.  Additional  investigations  and  consultation  are  required  to  establish  appropriate  standard 

practice. 

24 

Note  that  the  ZH  Instruments  SM-100  provides  for  variable  measurement  fields  between  10  and 
320  A/m.  Reported  spectra  were  measured  at  80  A/m  to  permit  consistent  comparison  between  the 
SM-100  and  the  single-field  (80  A/m)  Bartington  MS2B. 
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Initially,  for  SM-100  spectra  with  the  standard  35  mm  sample  holder,  the  approxi¬ 
mately  30  cc  sample  was  repeatedly  reconst itnted  from  a  single  100  ml  sonrce  sample 
to  assess  related  inflnence  of  sample  heterogeneity.  Corresponding  variation  in  mea¬ 
sured  susceptibility  (Figure  59)  is  approximately  ±0.25x10“®  m^/kg.  Following  these 
measurements,  the  hnal  sample  as  tested  was  set  aside  for  reference.  Subsequently, 
10  cc  measurements  were  similarly  acquired  on  samples  reconstituted  from  a  roughly 
50  ml  source  sample  with  expectation  of  proportionally  greater  measurement  fluc¬ 
tuation.  However,  despite  anticipation  that  sample  heterogeneity  would  have  more 
signihcant  influence,  it  is  observed  that  measurement  error  for  10  cc  sample  volumes 
is,  in  fact,  considerably  lower.  Unfortunately,  source  samples  were  not  consistent 
for  the  two  experiments,  suggesting  that  the  result  may  simply  reflect  lesser  hetero¬ 
geneity  in  the  10  cc  source  sample.  It  is  also  observed  that  the  SM-100  yields  a 
marginally  higher  susceptibility  than  the  MS2B.  To  ensure  consistent  comparison, 
subsequent  measurements  were  acquired  on  sub-samples  reconstituted  from  the  same 
100  ml  source  sample,  prior  to  sub-sample  archiving. 

As  expected,  results  for  the  DRDC  Austraila-A  soil  in  Figure  60  display  similar  and 
substantially  greater  susceptibility  variation  for  10  cc  measurements  with  both  SM- 
100  and  MS2B  instruments  compared  with  SM-100  measurements  for  the  instrument- 
standard  35  mm  sample  holder.  Here,  direct  inductance  measurements  are  also  in 
particularly  good  agreement  with  commercial  instruments,  yielding  a  measured  sus¬ 
ceptibility  of  approximately  Xp  =13.5x10“®  m^/kg,  roughly  consistent  with  the  mean 
of  MS2B  and  SM-100  measurements.  In  comparison  with  the  DRDC  Bosnia  soil,  the 
DRDC  Australia  soil  is  highly  susceptible  but  with  relatively  limited  frequency  de¬ 
pendence,  with  xfd  =1.0-1. 8%  versus  Xfd  =10.3-11.5%  for  the  Bosnia  soil.  It  is  no¬ 
table,  however,  that  respective  differential  susceptibility  (viscosity)  values  Ay;  =192- 
319x10“®  and  Ay;  =427-534x10“®  are  substantially  closer,  illustrating  the  combined 
influence  of  susceptibility  magnitude  and  the  related  level  of  dispersion.  With  respect 
to  CEN-WS7  (CEN,  2008)  guidelines  (Figure  62)  for  magnitude  of  susceptibility,  the 
predicted  inflnence  of  the  Australia-A  soil  (on  FDEM  detectors)  is  very  severe  while 
the  Bosnia  soil  is  only  moderate.  In  contrast,  according  to  related  guidelines  for 
differential  susceptibility  (influence  on  TDEM  detectors),  the  Australia-A  soil  has  a 
marginally  severe  to  very  severe  influence  while  the  Bosnia  soil  falls  clearly  within  the 
very  severe  category. 

Results  in  Figure  61  for  the  DRDC  Cambodia-2  reference  soil  reveal  yet  different  mag¬ 
netic  characteristics.  In  general,  it  is  observed  that  the  Cambodia-2  soil  possesses  a 
relatively  lower  mass-specihc  susceptibility  with  a  mean  value  less  than  Xp  =2.0x10“® 
m^/kg.  However,  a  substantially  higher  mass  density  p  ~1730  kg/m^,  compared 
with  p  ^985  kg/m^  for  the  DRDC  Bosnia  soil  means  that  corresponding  volume- 
specihc  susceptibility  x  ^346x10“®  is  considerably  higher  than  the  respective  value 
X~197  xlO  ®  for  the  Bosnia  soil,  although  both  soils  fall  within  the  moderate  CEN- 
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WS7  category  with  respect  to  influence  on  FDEM  metal  detectors.  Again,  it  is  in¬ 
structive  to  notice  that  while  the  associated  Xfd  =4.9-6.!%  is  intermediate  between 
DRDC  Bosnia  and  Australia-A  soils,  the  corresponding  viscosity  or  differential  suscep¬ 
tibility  Ax  =116-189x10“®  is  lowest.  The  explanation,  as  discussed  in  Section  3.6.6, 
is  that  despite  a  substantially  larger  underlying  frequency  dependence  than  for  the 
Australia-A  soil,  the  associated  susceptibility  is  nearly  an  order  of  magnitude  lower 
with  net  result  that  effective  apparent  viscosity  is  also  marginally  lower,  with  related 
effect  on  TDEM  sensors  categorized  as  severe  according  to  CEN-WS7  guidelines. 

Again,  as  expected,  variation  in  mass-speciflc  susceptibility  as  measured  by  MS2B 
and  SM-100  instruments  for  a  10  cc  sample  volume  ~  ±1.0x10“®  m^/kg  is  roughly 
similar  and  substantially  larger  than  the  corresponding  fluctuation  for  SM-100  mea¬ 
surements  with  the  instrument-standard  35  mm  (~30  cc)  sample  holder.  Measured 
susceptibility  via  the  direct  inductance  approach  is  again  in  good  agreement  with 
commercial  instruments.  In  general,  the  larger  sample  volume  of  the  solenoid  flxture 
limits  the  relative  influence  of  soil  heterogeneity  and  this  may  account  for  a  marginally 
lower  susceptibility  as  observed  for  the  DRDC  Bosnia.  However,  given  exclusion  of 
coarser-grained  soil  fractions,  there  also  remains  some  potential  that  the  assumed 
flxture  constant  A'/=1.053,  based  on  initial  MnCOs  calibration,  may  be  marginally 
low. 

Comparing  the  above  level  of  10  cc  measurement  error  with  a  corresponding  level 
of  ±10.0  xlO  ®  m^/kg  for  the  DRDC  Australia-A  soil,  it  is  tempting  to  conclude 
that  measurement  error  associated  with  sample  heterogeneity  is  directly  related  to  the 
magnitude  of  related  susceptibility.  Certainly,  there  is  a  direct  influence,  just  as  the 
influence  of  underlying  dispersion  on  net  magnetic  viscosity  is  modulated  by  the  level 
of  susceptibility.  However,  it  is  also  clear  that  the  scale  of  magnetic  heterogeneity, 
associated  texture,  sorting  and  related  compositional  characteristics  of  the  soil  can 
have  significant  influence. 

Related  effects  have  been  investigated  to  limited  extent  by  analyzing  the  susceptibil¬ 
ity  characteristics  of  specific  soil  textural  components.  In  particular,  a  direct  relation 
between  absolute  susceptibility  and  related  frequency  dependence  was  initially  noted 
in  connection  with  measurements  on  DRDC  Cambodia  reference  soils.  Although  the 
effect  is  notable  in  Figure  61  for  the  Cambodia-2  soil,  it  is  more  evident  for  the 
Cambodia-1  soil  in  Figure  63,  where  we  also  display  corresponding  results  for  five 
replicate  measurements  on  fine  and  coarse-grained  components  as  separated  via  a 
standard  ^10  (2.0  mm)  sieve.  As  anticipated,  it  is  observed  that  the  fine-grained 
soil  component  contributes  the  bulk  of  magnetic  viscosity,  having  both  a  higher  ab¬ 
solute  susceptibility  and  related  dispersion.  The  relationship  between  magnetic  char¬ 
acteristics  and  grain-size  is  characterized  in  the  lower  portion  of  Figure  63,  where 
density-normalized  viscosity  or  differential  susceptibility  Ax/p  and  frequency  depen¬ 
dence  parameter  xfd  are  charted  as  functions  of  mass-speciflc  susceptibility  (4.0  kHz 
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-  SMIOO  and  4.65  kHz  -  MS2B). 

As  anticipated  measnrements  for  coarse  and  fine-grained  separates  define  distinct 
fields,  with  the  hne-grained  fraction  yielding  a  more  tightly  dehned  data  held  than  the 
coarse-grained  fraction.  Of  conrse,  bnlk-sample  measurements  display  a  complicated 
but  systematic  mixing  relationship  with  respect  to  coarse  and  hne-grained  helds. 
Unfortunately,  however,  because  the  separates  were  not  derived  from  the  same  bulk 
sample,  the  range  of  bulk-sample  measurements  is  not  bracketed  by  end-member 
textural  helds  as  should  presumably  be  expected  for  consistent  samples. 

To  conhrm  this,  similar  experiments  were  undertaken  in  connection  with  characteri¬ 
zation  of  the  DRDC  Australia-B  soil,  with  all  measurements  acquired  on  subsamples 
from  a  consistent  bulk-soil  volume.  Results  are  displayed  in  Figure  64.  Here,  as  antic¬ 
ipated,  bulk-sample  measurements  are  roughly  bracketed  by  coarse  and  hne-grained 
end-member  helds.  Unexpectedly,  however,  it  is  the  coarse-grained  fraction  that  is 
the  primary  carrier  of  magnetic  viscosity  for  the  Australia-B  soil. 

To  further  investigate  the  ultimate  source  of  relatively  limited  magnetic  viscosity  for 
the  Australia-B  soil  and  noting  that  the  coarse-grained  fraction  included  soil  aggre¬ 
gates  and  soil-encrusted  rock  fragments,  the  coarse  grained  fraction  was  subsequently 
washed  via  distilled  water  to  remove  all  hne-grained  material.  Simultaneously,  the 
hne-grained  separate  was  further  separated  on  a  ^200-mesh  (75  //m)  sieve.  Result¬ 
ing  measurements  on  residual  fractions  and  adjusted  for  limited  sample  volume  are 
included  in  Figure  65  to  better  dehne  end- member  characteristics  associated  with 
bulk-sample  magnetic  properties.  In  general,  experimental  results  in  Figures  63-65 
illustrate  signihcant  variability  and  complexity  in  the  relationships  between  soil  textu¬ 
ral  composition,  related  magnetic  signatures  and  corresponding  inhuence  on  bulk  soil 
magnetic  properties.  Again,  it  is  noted  that  gravel  fractions  exceeding  approximately 
6.75  mm  (^3  mesh)  were  excluded  from  testing  and  that  further  consideration  must 
be  given  to  appropriate  instrumentation  (particularly  measurement  volume)  and  pro¬ 
cedures  to  ensure  that  laboratory  measurements  are  adequately  representative  of  in 
situ  soil  magnetic  characteristics. 

To  summarize  foregoing  assessment  of  DRDC  reference  soils  and  related  implications 
for  methodology  and  instrumentation,  we  highlight  the  following  hndings: 

1.  Tested  DRDC  reference  soils  appear  to  be  characterized  by  log-linear  frequency 
dependence  over  at  least  three  decades  of  frequency  between  100  Hz  and  100  kHz. 
Results  presented  here  and  by  others  (see  related  discussion  in  Section  3.4)  appear 
to  support  the  general  validity  of  basic  relations  set  out  in  Equation  (103)  and 
rehect  the  exponential  relationship  between  the  scale  of  magnetic  grains  at  the 
stable  single  domain,  superparamagnetic  transition  and  related  time-constants 
for  viscous  magnetization.  To  the  extent  that  emerging  consensus  on  this  issue 
continues  to  be  conhrmed,  dual  frequency  susceptibility  measurements  may  be 
perfectly  adequate  for  characterization  of  related  magnetic  viscosity. 
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2.  Soil  viscosity  v  =  Axj  iog{fH I fL)  for  Ih/Il  =10)  is  the  optimum 

performance  indicator  for  soil  influence  on  TDEM  sensors  and  reflects  both  the 
magnitude  and  frequency  dependence  of  soil  magnetic  susceptibility.  It  is  im¬ 
portant  to  appreciate,  however,  that  magnetic  soil  components  that  are  chiefly 
responsible  for  net  susceptibility  and  viscosity  are  not  necessarily  one  and  the 
same. 

3.  Soil  heterogeneity  can  yield  substantial  measurement  error  that  increases  with 

decreasing  sample  volume.  Consequently,  it  is  recommended  that  multiple  mea¬ 
surements  be  acquired  and  that  the  mean  values  be  reported  in  connection  with 
instruments  employing  smaller  measurement  volumes  (e.g.  10  cc).  It  is  also  rec¬ 
ommended  that  further  consideration  be  given  to  development  or  adaptation  of 
existing  instruments  for  accommodation  of  larger  volume  samples. 

4.  Influence  of  soil  heterogeneity  is  complicated  by  signihcant  complexity  in  the 
relationships  between  soil  textural  composition,  related  magnetic  signatures  and 
corresponding  influence  on  bulk  soil  magnetic  properties. 

5.  Despite  signihcant  exceptions,  frequency  dependence  of  susceptibility  is  less  prone 

to  volume-dependent  measurement  error  than  is  the  magnitude  of  susceptibility. 

4.5.2  Laboratory  Instrumentation: 

Following  sections  provide  additional  information  and  recommendations  regarding 
previously  described  instrumentation  for  multi-frequency  laboratory  characterization 
of  magnetic  susceptibility.  A  brief  description  of  alternative  instruments  is  also  pro¬ 
vided. 

Direct  Inductance  -  Sample  Fixture: 

Initial  design  of  a  prototype  solenoid  sample  hxture  for  direct  inductance  measure¬ 
ment  of  soil  magnetic  susceptibility  was  guided  by  two  primary  objectives:  a  sub¬ 
stantially  larger  sample  volume  than  provided  by  standard  instruments  and  as  broad 
as  possible  measurement  bandwidth.  As  previously  noted,  it  is  difhcult  to  achieve  a 
Hat  response  over  more  than  approximately  two  or  three  decades  and,  in  addition  to 
dimensional  characteristics,  band  limits  involve  tradeoffs  between  resistance,  induc¬ 
tance  and  effective  capacitance  of  coil  windings.  In  practice,  Hxture  design  is  guided 
by  use  of  appropriate  approximations  for  modelling  of  coil  impedance  parameters  and 
related  response  characteristics.  Fixture  design  should  also  take  advantage  of  unique 
performance  characteristics  of  the  speciHc  impedance  analyzer  or  LCR  meter  to  be 
used. 

In  addition,  it  should  be  appreciated  that  accurate  and  repeatable  complex  impedance 
measurements  require  frequent  compensation  adjustments  and  related  immobilization 
of  the  the  Hxture  and  connecting  leads  (particularly  beyond  the  extent  of  ground- 
referenced  shielding).  The  Hnal  rigidly  mounted  fixture  configuration  employed  here 
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is  illustrated  in  Figure  66,  including  expanded  view  of  a  special  purpose  connec¬ 
tor /switching  block  to  ensure  dependable  short  and  open  compensation  adjustments. 

Finally,  while  previous  results  illustrate  the  advantage  of  a  larger  measurement  volume 
in  limiting  measurement  uncertainty  associated  with  sample  heterogeneity,  a  smaller 
hxture  has  also  been  developed  for  use  where  relatively  limited  sample  is  available.  In 
particular,  the  small- volume  hxture  displayed  in  Figure  67  was  designed  to  accommo¬ 
date  the  same  35  mm  hlm-canister  sample  holder  employed  for  the  ZH  Instruments 
SM-100.  In  addition,  a  simple  adapter  has  been  fabricated  to  accommodate  the  stan¬ 
dard  10  cc  sample  vial  as  displayed  in  Figure  68  (also  for  use  with  ZHI  SM-100/105). 
Accurate  and  repeatable  open  and  short  compensation  adjustments  are  accomplished 
by  removing  the  socketed  coil  assembly  from  the  hxture  base  (open),  rotating  120 
degrees  and  reinserting  connectors  in  a  shorted  conhguration. 

Although  no  systematic  assessment  of  the  35  mm  hxture  has  been  possible  within  the 
scope  of  the  current  project,  an  initial  test  was  carried  out  on  a  volume-limited  sample 
of  Tiva  Canyon  Tuff  provided  by  USACE-ERDC  (North,  2006).  Results  are  depicted 
in  Figure  69.  First,  as  for  DRDC  reference  soils,  10  replicate  spectra  are  displayed  for 
SM-100  measurements  using  both  the  instrument- standard  35  mm  sample  hxture  and 
for  10  cc  samples  via  a  lOcc/35  mm  adapter.  Equivalent  dual-frequency  measurements 
are  displayed  for  the  MS2B  using  the  standard  10  cc  sample  vial  and  for  both  1.0  and 
0.1  ranges.  In  all  cases  measured  samples  were  reconstituted  from  the  same  75  ml 
source  sample.  Note  that  in  contrast  with  measurements  on  natural  soils,  replicate 
measurements  display  comparatively  little  variation  due  to  the  relatively  homogeneous 
nature  of  the  powdered  tuff. 

For  comparison,  hve  replicate  direct-inductance  susceptibility  spectra  were  acquired 
for  a  single  sample  of  the  Tiva  Canyon  tuff,  using  the  same  computer  controlled  ESl- 
2150  precision  LCR  meter,  in  this  case,  with  the  prototype  35  mm  sample  hxture. 
In  lieu  of  prior  calibration,  resulting  measured  spectra  were  subsequently  adjusted  to 
approximately  match  the  mean  susceptibility  measured  by  MS2B  and  SM-100  instru¬ 
ments  at  1  kHz.  In  effect,  an  approximate  calibration  was  obtained  for  the  35  mm 
hxture,  yielding  a  hxture  constant  of  R/  ~  1.5.  As  expected,  the  calibration  constant 
is  substantially  greater  than  unity  due  to  the  limited  aspect  (length/diameter)  ratio 
of  the  coil.  Assuming  effective  calibration  to  be  accurate,  however,  it  clear  that  the 
frequency  dependence  of  direct  inductance  spectra  are,  again,  in  very  good  agreement 
with  MS2B  and  SM-100  measurements. 

Specihcally,  it  is  observed  that  the  Tiva  Canyon  tuff  is  characterized  by  a  suscepti¬ 
bility  exceeding  Xp=1.0xl0“®  m^/kg  at  1  kHz  and  an  unusually  high  viscosity,  or 
differential  susceptibility,  approaching  Ax=2500xl0“®,  with  corresponding  xfd  ap¬ 
proaching  25  %.  Moreover,  there  is  some  indication  of  nonlinearity,  (i.e.  departure 
from  log-linear  relation)  particularly  at  higher  frequencies  that  was  not  observed  in 
connection  with  DRDC  reference  soils.  Viscosity  and  Xfd  indicated  in  Figure  69  are 
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determined  by  log-linear  regression  over  the  frequency  range  200  Hz  -  20  kHz.  Similar 
results  have  been  obtained  by  ZH  Instruments  (2006),  using  the  recently  available 
SM-105,  a  higher-frequency  equivalent  of  the  SM-100  with  measurement  frequencies 
extending  from  8  kHz  to  128  kHz. 

The  Tiva  Canyon  tuff,  also  referred  to  as  Yucca  Mountain  tuff,  is  a  well  characterized 
and  documented  volcanic  material,  incorporating  a  relatively  narrow  grain-size  dis¬ 
tribution  of  titanomagnetite  in  the  pseudo-single-domain/superparamagnetic  range 
(Worm  and  Jackson,  1999).  Both  single-domain  and  superparamagnetic  reference 
standards  have  been  prepared  via  selective  sampling  of  the  outcrop  (Jackson,  et  ah, 
2004;  Carter-Stiglitz,  et  ah,  2006).  It  is  our  understanding  (North,  2006)  that  tested 
material  is  a  powdered  and  homogenized  mixture  of  numerous  samples  and,  thus, 
presumably  has  a  broader  grain-size  distribution  than  previously  referenced  magnetic 
standards.  However,  in  general,  it  is  clear  that  the  frequency  extent  of  log-linear  fre¬ 
quency  dependence  decreases  as  the  magnetic  grain-size  distribution  narrows  and  Xfd 
increases.  It  is  notable  that  the  measured  value  of  Xfd  ~  25  for  the  Yucca  Mountain 
tuff  is  considerably  higher  than  the  commonly  accepted  upper  limit  of  Xfd  ~  15  for 
natural  soils  (see  related  discussion  in  Chapter  3).  Finally,  while  the  Tiva  Canyon  tuff 
is  certainly  a  useful  reference  material,  it  should  be  appreciated  (consistent  with  the¬ 
ory)  that  precise  magnetic  characteristics  are  signihcantly  dependent  on  temperature 
(e.g.  Worm  and  Jackson,  1999). 

Impedance  Analyzers  and  Multi-Frequency  LCR  Meters: 

As  noted,  direct  inductance  measurement  of  soil  magnetic  susceptibility  is  ultimately 
limited  by  the  dynamic  range  and  accuracy  of  inductance  measurements.  Conse¬ 
quently,  required  performance  specihcations  are  generally  more  restrictive  than  for 
electrical  conductivity  measurements.  In  particular,  we  recommend  that  the  instru¬ 
ment  provide  a  measurement  frequency  range  of  100  Hz  -  100  kHz,  with  minimum  of 
four  discrete  test  frequencies  spanning  the  range  (e.g.  100  Hz,  1  kHz,  10  kHz,  50-100 
kHz).  The  instrument  should  be  capable  of  providing  a  test  current  level  of  at  least 
10  mAjiMS  and  preferably  100  mAjiMS  for  specihed  open-coil  hxture  impedance. 
Although  it  is  desirable  that  the  test  current  level  be  regulated  with  sufficient  com¬ 
pliance  to  maintain  the  specihed  source  current,  it  is  more  often  the  case  that  only 
the  signal  voltage  level  (>1.0  Y rms)  is  controlled  with  resulting  test  current  depen¬ 
dent  on  load  impedance.  In  general,  it  is  recommended  that  the  source  signal  level 
(voltage/current)  be  adjusted  to  yield  a  moderate  low-held  intensity  of  approximately 
50-100  A/m. 

In  contrast  with  the  necessity  of  high  (>10.0  MQ)  input  impedance  for  measurement 
of  soil  electrical  conductivity,  the  instrument’s  output  impedance  characteristics  are 
more  signihcant  for  accurate  measurement  of  inductance  and  should  match  as  nearly 
as  possible  the  hxture’s  low-frequency  impedance.  Note  for  some  instruments,  output 
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impedance  varies  with  range  and  this  is  potentially  problematic  for  broadband  mea¬ 
surements  due  to  related  variation  in  accuracy  characteristics  of  the  instrument.  To 
assess  specihc  performance  potential  of  a  given  instrument,  however,  it  is  important 
to  consult  specihed  accuracy  formulae  or  tabulated  accuracy  data  as  a  function  of  hx- 
ture  inductance  and  measurement  frequency.  In  general,  the  instrument  should  have 
a  basic  accuracy  of  <  0.1%  and  a  display  resolution  of  at  least  5.5  digits.  In  practice, 
specihc  instrument  design  is  typically  tuned  to  yield  optimum  performance  for  either 
capacitance  measurements  or  inductance  measurements  rather  than  middle-ground 
characteristics  across  the  reactance  spectrum. 

Finally,  it  is  our  experience  that  computer  controlled  data  acquisition,  with  auto¬ 
mated  signal  averaging  and  drift  correction,  is  essential  for  successful  direct  induction 
measurement  of  soil  magnetic  susceptibility.  However,  it  should  be  cautioned  that 
the  implied  resolution  of  standardized  output  data  strings  can  exceed  instrument 
accuracy  dependent  on  range  and  frequency.  Again,  consult  related  accuracy  formu¬ 
las  or  tabulated  characteristics.  Ultimately,  as  previously  noted,  practical  limits  on 
measurement  accuracy  are  in  large  part  determined  by  the  precision  and  stability  of 
compensation  adjustments  and  related  manufacturer-specihed  procedures  should  be 
rigorously  implemented. 


Bartington  MS2B: 

The  Bartington  MS2B  susceptibility  meter  (Bartington,  2005)  is  depicted  in  Figure 
70.  The  standard  MS2B  system  comprises  a  MS2  measurement  console  and  a  MS2B 
dual  frequency  sensor,  with  measurement  frequencies  of  465  Hz  and  4.65  kHz.  Re¬ 
lated  held  intensity  is  80  A/m.  The  MS2B  is  a  stand-alone  instrument  with  integrated 
keyboard  and  display.  A  serial  interface  is  provided  for  computer  based  data  storage 
and  analysis.  As  previously  described,  the  MS2B  measures  magnetic  susceptibility  by 
gauging  an  associated  shift  in  the  selected  oscillator  frequency.  The  standard  mea¬ 
surement  sequence,  driven  by  front  panel  switches,  includes  an  initial  null  or  “ZERO” 
reading  with  the  sensor  coil  open,  followed  by  a  “MEASURE”  reading  with  the  sample 
positioned  within  the  sensor  coil.  For  improved  accuracy,  the  manufacturer  recom¬ 
mends  a  subsequent  open-coil  air  “MEASUREMENT”  reading  xa  following  the  sam¬ 
ple  “MEASUREMENT”  reading  xs  fo  control  instrument  drift  over  the  measurement 
interval.  The  corrected  volume  susceptibility  is  x  =  —  {xaI‘^)- 

Measurement  precision/resolution  is  selectable  from  xl  and  xO.l  ranges.  Corre¬ 
sponding  measurement  periods  are  1.2  s  and  12  s,  respectively  (presumably,  the  xO.l 
range  measurement  represents  a  10-fold  average  reading)  In  both  cases,  measured 
volume-specihc  susceptibility  is  displayed  in  units  of  xl0“^  SI  and  in  the  latter  case 
with  0.1x10“®  SI  accuracy  (readings  can  also  be  displayed  in  cgs  units).  Measure¬ 
ment  frequency  (/l=465  Hz  or  /h=4.65  kHz)  is  similarly  controlled  via  a  front  panel 
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switch.^® 


The  diameter  of  the  MS2B  sensor  core  is  approximately  36  mm  and,  as  per  previons 
discnssion,  the  standard  sample  holder  is  a  10  cc  cylindrical  vial.  However,  it  is  em¬ 
phasized  that  non-standard  containers  can  also  be  nsed  with  appropriate  adjustment 
of  resulting  measurements.  The  Bartington  operation  manual  provides  guidance  on 
the  use  of  non-standard  sample  holders  and  related  required  adjustments,  including 
centering  of  the  platen  on  which  the  sample  is  lowered  into  the  sensor  coil.  It  is  also 
noted  that  the  MS2  can  be  operated  with  a  larger  volume,  single-frequency  sensor. 
In  particular,  the  MS2B62  sensor  core  has  interior  diameter  of  approximately  62  mm 
and  accommodates  a  standard  250  cc  cylindrical  sample  holder.  The  MS2B62  sensor 
operates  at  a  nominal  frequency  of  1.64  kHz  with  a  held  intensity  of  approximately 
80  A/m,  consistent  with  the  MS2B. 

Figure  71  illustrates  the  results  of  a  simple  experiment  to  investigate  the  sensitivity 
characteristics  of  the  MS2B  sensor  and  related  signihcance  in  connection  with  use  of 
non-standard  sample  holders.  Initially,  as  per  standard  operating  procedure,  sensor 
platen  height  was  adjusted  to  yield  maximum  reading  (0.1  scale)  for  the  1%  Fe304 
calibration  standard  (it  is  noted  that  improved  adjustment  would  likely  be  obtained 
by  utilizing  a  more  highly  susceptible  standard).  Subsequently,  the  susceptibility  of  a 
standard  10  cc  sample  pot  was  measured  and  recorded  with  incremental  addition  of  a 
MnCOa  standard.  In  particular,  1.00  g  increments  were  added  and  leveled  with  due 
care  to  minimize  any  related  compaction  of  the  composite  sample.  Fortuitously,  1.00 
g  increments  of  the  MnCOs  standard  were  approximately  equivalent  to  1  cc  volume 
increments  (i.e.  pb  ~  1  g/cc).  The  resulting  linear  relation  between  volume  specihc 
susceptibility  (SI)  and  approximate  sample  volume  is  displayed  at  the  upper  right  of 
Figure  71. 

To  investigate  the  inhuence  of  a  non-standard  sample  vial,  the  experiment  was  sub¬ 
sequently  repeated  utilizing  the  marginally  smaller  diameter  and  taller  sample  pot 
pictured  at  the  upper  left  of  Figure  71.  Measured  physical  characteristics  of  the  two 
vials  are  displayed  for  reference.  In  this  case,  the  MnCOa  standard  was  added  in  0.50 
g  increments.  The  corresponding  increase  as  a  function  of  approximate  sample  volume 
is  displayed  at  the  lower  left  of  Figure  71,  with  the  standard  10  cc  pot  response  (solid 
line)  overlaid  for  reference.  It  is  noted  that  the  relation  is  again  linear  to  an  approx¬ 
imate  sample  volume  of  roughly  8  cc.  Note  that  for  the  non-standard  sample  vial,  a 
sample  level  of  2.3  cm,  the  nominal  height  of  the  standard  10  cc  pot,  is  equivalent  to 
a  sample  volume  of  approximately  8  cc.  Consequently,  it  is  apparent  that  the  height 

Note  that  a  sufHcient  stabilization  period  should  be  permitted  following  change  of  measurement 
frequency  (experimentation  suggests  approximately  1  min).  The  manufacturer  recommends  that  for 
dual  frequency  measurements  on  multiple  samples,  all  samples  should  be  tested  at  a  given  frequency 
and,  subsequently,  retested  at  the  second  frequency  following  adequate  settling  period. 
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of  the  Bartington  10  cc  sample  pot  was  designed  to  maintain  the  sample  within  an 
optimnm  range  of  the  sensor  midpoint  within  which  the  sensor’s  volnme  sensitivity 
is  ronghly  linear.  In  addition,  it  is  noted  that  the  slope  of  the  linear  portion  of  the 
sensitivity  cnrve  for  the  non-standard  pot  is  slightly  less  that  for  the  standard  10  cc 
pot,  snggesting  that  within  the  linear  range,  volnme  sensitivity  is  slightly  lower  for 
the  smaller  diameter  sample  vial. 

In  general,  experimental  results  indicate  that  for  a  standard  10  cc  sample  vial,  volume 
susceptibility  for  a  partially  hlled  sample  holder  can  be  corrected  by  the  factor  h/ho, 
where  h  denotes  an  arbitrary  height  and  ho  =  2.3  cm,  or  by  the  equivalent  volume 
factor  VjVQ.  In  general,  a  similar  approach  is  applicable  for  any  sample  holder  within 
the  linear  sensitivity  region  and  with  symmetry  about  the  sensor  axis.^®  However,  it 
is  emphasized  that  linear  volume  sensitivity  assumes  uniform  mass  density. 

Indeed,  within  the  linear  sensitivity  limit,  the  sensor  is  effectively  indifferent  to  volume 
or  equivalent  mass  sample  increments.  For  example,  with  reference  to  the  foregoing 
experiment;  if  additional  sample  increments  beyond  8  cc  were  compacted  into  the 
non-standard  vial,  such  that  the  composite  sample  level  remained  below  the  2.3  cm 
limit,  corresponding  susceptibility  readings  would  increase  an  amount  consistent  with 
effective  extension  of  the  linear  sensitivity  trend.  Consequently,  for  an  arbitrary 
sample  holder  within  the  linear  sensitivity  region  and  with  symmetry  about  the  sensor 
axis,  measured  susceptibility  is  linearly  related  to  sample  mass  as  well  as  sample 
volume.  In  fact,  in  view  of  the  volume  biasing  effects  of  variable  compaction,  it  is 
generally  recommended  that  raw  susceptibility  readings  be  converted  to  mass-specihc 
values  via  concurrent  measurement  of  sample  mass  (density)  and  that  mass-specihc 
susceptibility  be  the  primary  reported  parameter,  together  with  mass  density.  Having 
said  this,  however,  it  should  also  be  appreciated  that  sample  mass  can  be  skewed  by 
variable  soil  moisture  content  and,  consequently,  samples  should  be  thoroughly  air- 
dried  or  oven  dried  at  less  than  40°  C  prior  to  analysis. 

In  general,  a  non-standard  sample  vial  with  length  extending  beyond  the  linear  sen¬ 
sitivity  region  (assuming  symmetry  about  the  sensor  access)  can  be  utilized  with 
appropriate  calibration  via  a  known  reference  standard.  The  appropriate  correction 
factor  is  7  =  xixoi  where  y  is  the  measured  susceptibility  of  the  arbitrary  standard- 
hlled  sample  vial  and  xo  denotes  the  measured  susceptibility  of  the  standard-hlled  10 
cc  pot.  It  is  emphasized  that  calibration  measurements  for  arbitrary  and  standard 
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For  maximum  accuracy  with  partially  full  or  uou-staudard  sample  vials,  the  MS2B  plateu 
height  should  be  adjusted  for  maximum  readiug  prior  to  subsequeut  correctiou  by  fractioual  volume 

V/Vo- 
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Note  that  soil  maguetic  susceptibility  for  viscous  maguetic  soils  iucorporatiug  siguificaut  SD- 
SP  couteut,  is  temperature  depeudeut  aud  samples  should  be  allowed  to  stabilize  at  laboratory 
temperature  prior  to  aualysis  (also  to  avoid  thermal  affect  ou  measuremeut  iustrumeutatiou). 
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sample  pots  must  allow  for  readjustment  of  the  platen  level  for  appropriate  vertical 
centering.  Moreover,  the  maximum  height  of  an  arbitrary  sample  pot  is  restricted  by 
maximum  travel  of  the  platen  (±1.0  cm  travel  and  ho =2. 3  cm  yields  maximum  vial 
height  of  4.3  cm  or,  allowing  sufficient  play  for  vertical  centering,  approximately  4.0 
cm). 

To  further  establish  the  sensitivity  characteristics  of  the  MS2B  sensor,  both  axial  and 
transaxial  held  variation  were  investigated.  First,  axial  sensitivity  was  established 
by  incrementally  increasing  the  height  of  the  1%  Fe304  calibration  standard  above 
the  platen  (sensor  midpoint).  The  diagram  at  upper  right  of  Figure  72  displays  the 
variation  of  measured  volume  susceptibility  (SI)  with  increasing  sample  height  (base 
of  sample  pot  above  platen)  in  0.5  cm  intervals.  Note  that  when  the  base  of  the 
sample  pot  is  located  at  the  upper  housing  surface  level  (approximately  4.2  cm),  the 
measured  susceptibility  Xv  =  13.2  x  10“^  SI  is  approximately  4%  of  the  calibrated 
value  Xv  =  346.9  x  10“^  SI  at  the  sensor  midpoint  (platen  level). 

In  effect,  the  observed  reduction  in  measured  susceptibility  reflects  removal  of  an 
incremental  sample  volume  from  the  sensor’s  linear  sensitivity  region  as  previously 
discussed.  As  a  practical  matter,  it  is  noted  that  calibration  of  the  sensor’s  vertical 
sensitivity  permits  approximate  measurement  of  samples  having  volume  specific  sus¬ 
ceptibilities  beyond  the  standard  operating  range  of  the  instrument.  For  example, 
the  approximate  susceptibility  of  DRDC  magnetite  soils  was  measured  by  positioning 
a  standard  10  cc  sample  on  a  sheet  of  transparent  material  placed  over  the  sensor 
cavity  and  centered  on  the  sensor’s  axis.  Results  obtained  by  calibrated  scaling  of 
the  raw  measured  value  were  in  adequate  agreement  with  previous  direct-inductance 
measurements . 

Finally,  in  connection  with  Figures  71-72,  it  was  previously  noted  the  slope  of  the 
sensitivity  curve  for  the  non-standard  sample  pot  was  marginally  less  than  that  for 
the  standard  10  cc  pot.  It  was  anticipated  that  this  deviation  might  reflect  variation 
of  held  intensity  with  radial  offset  from  the  sensor  axis.  To  investigate  the  transaxial 
sensitivity,  a  small  yet  responsive  probe  was  fabricated  by  fixing  approximately  0.3  cc 
of  highly  susceptible  Fe203  in  the  cylindrical  reservoir  of  a  I.O  cc  syringe  (see  Figure 
73).  Related  susceptibility  measurements  were  subsequently  acquired  with  probe 
(oriented  vertically)  at  six  incremental  positions  across  the  diameter  of  the  platen 
as  illustrated  in  Figure  73.  The  process  was  subsequently  repeated  with  the  platen 
repositioned  0.5  cm  and  I.O  cm  above  the  sensor  midpoint.  Resulting  variation  of 
measured  susceptibility  is  displayed  at  the  lower  right  of  Figure  72.  It  is  observed  that 
in  all  cases,  apparent  sensitivity  is  minimum  at  the  axis  and  increases  marginally  with 
radial  offset  and  more  significantly  with  increasing  height  above  the  zero  referenced 
platen.  It  is  notable  that  measured  transaxial  sensitivity  variation  is  consistent  with 
the  marginally  lower  volume  sensitivity  observed  for  the  smaller  diameter  sample  pot 
in  Figure  71. 
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ZH  Instruments  SM-100/105: 

The  ZH  Instruments  SM-100  susceptibility  meter  (ZH  Instruments,  2007)  is  pictured 
in  Figure  74.  In  contrast  with  the  stand-alone  Bartington  MS2B  system,  control  and 
display  functions  of  the  SM-100  are  software-implemented  via  an  integrated  RS-232 
interface.  The  software  control/display  panel  is  displayed  in  Figure  75.  The  SM-100 
provides  substantially  more  flexibility  than  the  Bartington  system,  including  multiple 
test  frequencies  (500  Hz,  1  kHz,  2  kHz,  4  kHz,  8  kHz)  and  a  range  of  held  intensity 
levels  (10  A/m,  20  A/m,  40  A/m,  80  A/m,  160  A/m  and  320  A/m).  An  otherwise 
equivalent  system,  the  SM-105,  offers  an  extended  frequency  range  (8  kHz,  16  kHz, 
32  kHz,  64  kHz,  128  kHz)  and  a  more  restricted  range  of  held  intensity  levels  (10 
A/m,  20  A/m,  40  A/m,  80  A/m).  Test  frequency  and  held  intensity  parameters 
are  software  selectable  via  a  pull-down  menu.  In  addition,  SM-100/105  instruments 
include  an  integrated  scale  with  0.01  g  resolution,  allowing  real-time  measurement  of 
mass-specihc  susceptibility. 

Although  the  SM-100  ohers  substantially  greater  hexibility  in  test  held  characteristics, 
the  fundamental  design  and  measurement  sequence  is  similar  to  that  described  for  the 
MS2B.  In  particular,  like  the  MS2B,  the  SM-100  measures  magnetic  susceptibility  by 
gauging  an  associated  shift  in  the  selected  oscillator  frequency  that  results  on  inserting 
the  soil  sample  into  the  sensor  cavity.  Unlike  the  MS2B,  however,  the  measurement 
sequence  is  computer  automated  with  a  combination  of  display  and  audible  prompts 
cuing  the  operator  to  insert  and  remove  the  sample  and  with  automatic  drift  correction 
implemented  in  real-time.  The  corresponding  measurement  sequence  and  related  drift 
correction  is  diagrammatically  illustrated  in  Figure  75. 

As  previously  discussed,  the  SM-100/105  employs  a  standard  sample  holder  (35  mm 
him  canister)  having  approximately  3  times  the  volume  of  the  Bartington-standard  10 
cc  vial  and,  depending,  on  the  scale  and  extent  of  soil  heterogeneity,  this  is  potentially 
advantageous.  A  simple  adapter  to  accommodate  the  standard  10  cc  sample  holder 
is  displayed  in  Figure  68.  In  general,  although  project  scope  did  not  permit  a  parallel 
investigation  of  SM-100  sensor  characteristics,  it  is  reasonable  to  presume  that  the 
sensor  has  grossly  similar  sensitivity  characteristics  and  a  corresponding  linear  (i.e. 
uniform  held)  sensitivity  region.  Indeed,  measurements  acquired  for  10  cc  samples  via 
the  10  cc  to  35  mm  adapter  were  simply  scaled  by  the  relative  volume  factor  V35/V10, 
with  resulting  excellent  agreement  for  standards.  Again,  for  maximum  accuracy  it  is 
necessary  that  the  sample  volume  for  a  partially  full  or  non-standard  sample  holder 
be  approximately  centered  with  respect  to  the  sensor  axis  and  the  vertical  “plane 
of  symmetry”  as  illustrated  in  Figure  76.  As  indicated,  two  “weighting  plates”  are 
supplied  with  the  instrument,  one  that  centers  a  full  35  mm  sample  holder  and  one 
that  centers  a  ha, If- full  35  mm  sample  holder.  Note  that  the  10  cc  to  35  mm  adapter 
was  designed  to  center  the  10  cc  sample  pot  within  the  standard  35  mm  holder  and 
is,  consequently,  used  with  the  standard  (full  holder)  “weighting  plate”. 
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Finally,  no  qnantitative  assessment  has  been  made  of  accnracy  for  the  integrated 
strain-gange  mass  scale  and  all  reported  susceptibility  measurements  are  based  on 
independent  measurements  of  sample  mass  via  a  Mettler  PE-3600  (0.01  g)  bench  scale. 
However,  general  comparison  with  corresponding  SM-100  measured  mass  suggests 
agreement  to  0.05-0.1  g. 

Alternative  Laboratory  Instrumentation: 

Mention  should  also  be  made  of  another  recently  developed  and  commercially  avail¬ 
able  system,  the  Magnasat  developed  by  Ultradynamics  in  Australia  (Ultradynamics, 
2008,  North,  2006)  and  a  prototype  system  developed  by  the  Geophysics  Laboratory 
at  the  University  of  Toronto  (West  and  Bailey,  2005).  Both  instruments  employ  a 
null  or  balanced  coil  system  and  yield  frequency-dependent  complex  susceptibility  via 
phase- sensitive  detection  of  the  out-of-balance  signal  produced  on  insertion  of  the  soil 
sample.  Similar  instruments  have  been  previously  described  by  Strangway  (1967), 
Mullins  and  Tite  (1973)  and  Dabas,  et  al.  (1992).  In  contrast  with  the  direct  in¬ 
ductance  approach,  a  null  or  balanced  coil  apparatus  offers  enhanced  resolution  by 
devoting  the  full  extent  of  system  dynamic  range  to  the  out-of-balance  signal,  rather 
than  detecting  a  corresponding  deviation  of  similar  scale  in  a  primary  signal  several 
orders  of  magnitude  larger.  A  simple  proof  of  concept  experiment  has  demonstrated 
that  similar  balanced-coil  measurements  can  be  acquired  using  an  appropriate  dual¬ 
coil  sample  hxture  and  an  off-the-shelf  gain-phase  analyzer. 

There  is  also  a  range  of  single  frequency  magnetic  susceptibility  bridges  available  that 
could  be  used  in  combination  to  assess  frequency  dependence  of  soil  magnetic  suscep¬ 
tibility.  However,  where  such  an  approach  is  taken,  instruments  should  accommodate 
the  same  sample  holder  and  due  consideration  should  be  given  to  appropriate  and 
precise  inter-calibration.  In  addition,  measurement  held  intensities  should  be  roughly 
consistent. 

Finally,  recent  development  of  special  purpose  soil  magnetic  viscosity  meters  should 
be  noted.  In  contrast  with  foregoing  frequency- domain  instrumentation,  viscosity  me¬ 
ters  are  time-domain  devices  and  measure  viscosity  directly  by  exposing  the  sample 
to  a  sustained  primary  magnetic  held,  subsequently  terminating  the  primary  held 
and  gauging  the  rate  of  decay  of  the  transient  secondary  held  associated  with  elec¬ 
tromagnetic  polarization  of  the  sample.  A  prototype  system  has  been  described  by 
Siegenfeld  (2002)  and  a  commercially  developed  instrument,  including  a  held  sensor, 
is  more  recently  available  (Pulsepower  Developments,  2008).  Measurement  windows 
are  in  the  range  1-100  /rs. 

Although  direct  time-domain  measurement  of  electromagnetic  viscosity  is  clearly  of 
signihcant  interest  in  connection  with  assessment  of  TDEM  sensor  performance,  it 
is  noted  that  foregoing  devices  provide  no  direct  indication  of  absolute  magnetic 
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susceptibility.  No  testing  of  these  instruments  has  been  undertaken  as  part  of  the 
current  study. 

4.5.3  Field  Methods,  Procedures  and  Instrumentation: 

Hand-Held  Susceptibility  Meter: 

Magnetic  susceptibility  and  its  frequency  variation  can  be  characterized  in  situ  using 
of  a  range  of  special  purpose,  commercially  manufactured  instruments.  In  general, 
held  susceptibility  meters  are  based  on  the  same  tuned  LC  oscillator  design  and  as 
employed  by  MS2B  and  SM-100  laboratory  systems.  In  fact,  the  Bartington  MS2 
console  can  be  conhgured  with  a  handle-mounted  MS2D  search  coil  (Figure  77  -  left) 
for  single  frequency,  in  situ  measurement  of  soil  magnetic  susceptibility.  A  range  of 
hand-held  instruments,  including  the  Kappameter  KT-6  (Figure  77  -  right)  are  also 
available. 

In  place  of  the  cylindrical  cavity  sensor  coil  employed  for  laboratory  systems,  a  pla¬ 
nar  or  low-prohle  sensor  coil  comprises  the  inductive  element  in  the  meter’s  oscillator 
circuit.  The  instrument  is  initially  nulled  in  air  and  subsequently  brought  into  planar 
contact  with  the  soil  under  test  for  relative  measurement  of  associated  susceptibility. 
Although  we  are  not  aware  of  dual  or  multi- frequency  devices  for  in  situ  measurement, 
instruments  are  available  having  a  range  of  operating  frequencies  so  that  frequency 
variation  of  susceptibility  can  be  established  by  separate  measurements  with  instru¬ 
ments  operating  at  two  or  more  frequencies. 

Unfortunately,  due  to  the  variable  geometry  and  limited  scale  of  typical  sensor  coils 
(<  10-20  cm),  effective  measurement  volume  can  differ  signihcantly  and  is  generally 
restricted  to  <  5-10  cm.  For  example.  Figure  78  displays  the  measured  depth  sensi¬ 
tivity  for  MS2D  and  KT-5  (KT-6  equivalent)  as  reported  by  Lecoanet,  et  al.  (1999). 
As  a  consequence  of  restricted  measurement  volume,  the  effective  precision  of  these 
instruments  is  highly  sensitive  to  both  soil  heterogeneity  and  surface  irregularity.  To 
minimize  related  influence,  the  soil  surface  should  be  prepared  to  ensure  uniform  and 
continuous  sensor  contact  and  measurements  should  avoid  (as  possible)  anomalous 
soil  heterogeneity  to  yield  readings  reflective  of  mean  soil  characteristics. 

In  general,  as  for  in  situ  electrical  conductivity  measurements,  it  is  recommended  that 
multiple  measurements  and  appropriate  spatial  sampling  be  pursued  to  characterize 
mean  susceptibility  and  related  spatial  variation.  Specihcally,  a  series  of  repeat  mea¬ 
surements  at  each  sampling  location  with  the  sensor  position  re-located  randomly 
by  a  distance  comparable  to  sensor  diameter  yields  useful  indication  of  a  represen¬ 
tative  value  and  associated  deviation.  To  establish  the  depth-dependent  variation 
of  magnetic  susceptibility  a  sequence  of  measurements  should  be  acquired  at  incre¬ 
mental  depth  within  one  or  several  test  pits.  Finally,  Although  bulk  soil  magnetic 
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susceptibility  is  typically  controlled  by  magnetic  content  of  the  fine-grained  soil  frac¬ 
tion,  the  relative  susceptibility  of  larger  gravel  and  cobble-sized  grains  should  also  be 
investigated  to  assess  potential  for  related  magnetic  clutter. 

In  contrast  with  electrical  conductivity  measurements,  variable  soil  moisture  has  only 
limited  diamagnetic  influence  on  soil  magnetic  properties  and,  consequently,  assess¬ 
ment  and  monitoring  of  short-term  temporal  variation  is  not  ordinarily  required.^®  In 
general,  temperature  dependence  is  also  limited,  although  potentially  signihcant  for 
highly  viscous  soils  under  strong  seasonal  temperature  variation. 


FDEM  Susceptibility  Meters: 

Magnetic  susceptibility  can  also  be  measured  via  larger  dual-coil  FDEM  devices  as 
previously  described  for  measurement  of  soil  electrical  conductivity  (see  Sections  4.4.4 
and  Figure  52).  Instead  of  gauging  the  quadrature  response,  however,  magnetic  sus¬ 
ceptibility  is  derived  from  the  in-phase  signal  component  under  the  same  low  induction 
number  (LIN)  assumption.  In  particular,  as  discussed  in  Section  3.1  and  with  refer¬ 
ence  to  Equation  (29)  and  Figure  53,  it  is  observed  that  dependent  on  coil  separation 
s  and  for  a  sufficiently  low  frequency- conductivity  fa  product,  the  induction  number 
rj  falls  within  the  so-called  resistive  limit  and  the  in-phase  response  is  almost  entirely 
dependent  on  magnetic  susceptibility  and  independent  of  other  parameters.  In  gen¬ 
eral,  for  soils  having  ordinary  electrical  conductivity  {a  <100  mS/m)  and  elevated 
magnetic  susceptibility  (x  >  1-0  ^  10”^);  small-scale  FDEM  devices  are  typically 
operating  within  the  resistive  limit,  particularly  at  lower  operating  frequencies,  and 
related  magnetic  susceptibility  is  derived  directly  as 


X  =  li 


(165) 


where  qj  denotes  the  appropriate  instrument-specihc  scaling  factor. 

Again,  it  should  be  appreciated  that  the  depth  sensitivity  of  related  measurements  is 
strongly  dependent  on  coil  orientation  and  separation  and  that  the  effective  measure¬ 
ment  volume  varies  accordingly.  In  analogy  with  Equations  (155)  and  (156)  for  the 
quadrature  sensitivity  of  coplanar  coil  systems,  corresponding  expressions  for  in-phase 

QO 

Note  also,  experimental  investigations  by  Maier,  et  al.  (2006)  appear  to  confirm  that  eddy- 
current  effects  associated  with  moisture-dependent  soil  electrical  conductivity  have  little  or  no  prac¬ 
tical  influence  on  the  operation  of  common  hand-held  susceptibility  meters  operating  at  frequencies 
below  10  kHz. 

OQ 

For  conductive  soils  having  induction  numbers  significantly  beyond  the  resistive  limit  for  given 
instrument  parameters,  the  in-phase  response  [Hg / Hp') j  (8/15)  is  relatively 

limited  and  effectively  frequency-conductivity  fa  dominated  (see  McNeill  and  Bosnar,  1999). 
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sensitivity  are  (McNeill  and  Bosnar,  1999) 


^  12Z  (8Z2  _  3) 

ShC-i{  )  ^^^2  +  1)7/2 

(166) 

SvC-l{Z)  -  ^  ^^5/2 

(167) 

for  horizontal  coplanar  (HC)  and  vertical  coplanar  (VC)  systems,  respectively.  As¬ 
sociated  depth  sensitivities  are  depicted  in  the  npper  panel  of  Fignre  79,  again,  as 
fnnctions  of  s-normalized  depth  [Z  =  z/s).  Notably,  and  consistent  with  Fignre  53, 
the  horizontal  coplanar  sensitivity  is  negative-valned.  More  signihcantly,  however, 
in-phase  depth  sensitivities  for  both  horizontal  and  vertical  coplanar  conhgnrations 
are  well  focnsed,  providing  considerably  better  resolntion  of  near-surface  conditions 
than  corresponding  quadrature  conductivity  sensitivities  in  Figure  54.  Correspond- 
ing  integrated  sensitivities  S(Z)  =  S(x)dx  in  Figure  79  indicate  that  70%  of 
cumulative  sensitivity  for  horizontal  coplanar  coils  is  attained  for  .2  <  0.25s  and  for 
vertical  coplanar  coils  only  marginally  in  excess  of  ^  =  0.5s.  Consequently,  compared 
with  sensitivity  characteristics  for  quadrature  conductivity  measurements  (Figure  54), 
depth  of  in-phase  investigation  for  standard  FDEM  instruments  is  better  suited  for 
present  purposes. 

As  described  in  Section  4.4.3  (in  connection  with  FDEM  conductivity  measurements 
and  in  analogy  with  quadrapole  electrical  soundings),  vertical  electromagnetic  sound¬ 
ings  can  be  acquired  by  recording  a  series  of  measurements  with  the  instrument  raised 
incrementally  above  grade.  Computer-based  inversion  methods  are  applied  to  result¬ 
ing  measurements  to  derive  an  optimized  estimate  of  actual  depth-dependent  magnetic 
susceptibility.  Similarly,  semi-quantitative  assessment  of  depth- variability  is  obtained 
by  recording  the  variation  in  measured  in-phase  response  as  the  device  is  raised  from 
grade  to  a  height  of  approximately  1.5s.  For  a  uniform  soil,  in-phase  response  should 
decrease  with  increasing  height  in  accordance  with  integrated  sensitivity  relations  in 
the  lower  panel  of  Figure  79.  In  particular,  s-weighted  height  dependent  response 
above  a  uniform  ground  is  given  by  R{H)  =  S{x)  dx,  where  H  =  h/s  and  h  is 
height  above  grade  (corresponding  response  characteristics  are  also  depicted  in  Figure 
79). 

Finally,  although  the  in-phase  FDEM  response  (Figure  53)  due  to  induction  is  effec¬ 
tively  frequency  independent  for  soils  within  the  resistive  limit,  the  intrinsic  frequency 

on 

Integrated  sensitivities  follow  from  Equations  (166)  and  (167)  as 

1:HC-i{Z)  =  [(1  -  8Z2)/(4Z2  +  1)5/2]  _  1  i:yc-i{Z)  =  -1/(4Z2  +  1)3/2. 

O  1 

Resulting  response  functions  follow  from  Equations  (166)  and  (167)  as 

Rhc-i{H)  =  {8H^  -  +  1)^/2  Rvc-i{H)  =  1/{AH^  +  1)3/2. 
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dependence  of  soil  electromagnetic  properties,  particnlarly  the  magnetic  snsceptibility, 
can  have  a  snbstantial  inflnence.  Conseqnently,  as  for  FDEM  condnctivity  meters, 
mnlti- frequency  FDEM  susceptibility  measurements  can  potentially  provide  assess¬ 
ment  of  magnetic  dispersion,  associated  viscosity  and  related  influence  on  time-domain 
(TDEM)  sensors  as  described  in  Section  3.3. 
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5.0  SOIL  NOISE  -  SUMMARY  AND  CONCLUSIONS: 


As  described  in  Section  3.0  and  with  reference  to  Fignre  4,  the  net  response  of  both 
freqnency- domain  (FDEM-CW)  and  time-domain  (TDEM-PI)  metal  detectors  re¬ 
flects  both  backgronnd  and  anomalons  components,  with  the  anomalons  component 
comprising  both  signal  associated  with  targets  of  interest  (mines)  and  noise  or  clnt- 
ter  dne  to  localized  soil  heterogeneity.  Moreover,  to  the  extent  that  compensation 
circuitry  is  not  provided  or  is  rendered  ineffective  due  to  exceptionally  strong  back¬ 
ground  response  or  large-scale  gradients,  the  background  component  is  also  effectively 
noise  and  potentially  imposes  a  limiting  restriction  on  system  performance. 

In  general,  reduced  signal  to  noise  ratio  implies  degraded  performance.  However,  it  is 
important  to  appreciate  that  signal  to  noise  ratio  is  reduced  both  by  decreasing  signal 
and/or  increasing  noise.  In  practice,  signal  levels  associated  with  landmines  depend 
principally  on  the  amount  (surface  area)  of  incorporated  metal  and  depth  of  burial. 
Composition,  shape  and  orientation  of  metal  components  also  have  influence.  Clearly, 
deeply  buried,  low-metal  mines  are  the  ultimate  test  of  metal  detector  performance. 
However,  accepting  practical  limitations  on  signal,  the  present  study  has  focused  on 
the  origin,  nature  and  extent  of  soil  noise,  methods  for  characterization  of  associated 
soil  electromagnetic  properties  and  related  prediction  of  metal  detector  performance. 

Background  Soil  Noise: 

Again,  there  are  two  basic  components  of  soil  noise,  background  and  anomalous.  By 
background  noise,  we  refer  to  the  broad-scale  response  associated  with  mean  soil  elec¬ 
tromagnetic  properties  and,  specihcally,  to  background  levels  that  are  not  effectively 
cancelled  or  controlled  by  compensation  circuitry  or  related  adjustments  (manual  or 
automated).  In  general,  theoretical  analysis  and  related  modelling  presented  in  Chap¬ 
ter  3.0,  has  conhrmed  the  primary  role  of  soil  magnetic  characteristics  in  limiting  the 
performance  of  both  FDEM  and  TDEM  metal  detectors.  In  particular,  it  is  observed 
that  for  typical  coil  geometries  and  FDEM  operating  frequencies,  the  related  induc¬ 
tion  number  (Equation  29)  satishes  r]  «  1.0  for  the  vast  majority  of  natural  soils 
(Figure  9).  Moreover,  for  soils  having  substantial  magnetic  susceptibility  exceeding 
X  =  1.0  X  10“^  —  1.0  X  10“^  SI,  the  induction  number  commonly  falls  within  the 
resistive  limit  (Figure  10),  with  related  FDEM  response  dominated  by  the  in-phase 
component  and  proportional  to  magnetic  susceptibility.  Clearly,  where  soil  magnetic 
susceptibility  is  sufficiently  high,  compensation  is  ineffective  and  the  background  soil 
noise  level  exceeds  normal  signal  ranges.  Ultimately,  with  increasing  susceptibility,  the 
background  response  effectively  reduces  and  ultimately  saturates  available  dynamic 
range,  rendering  the  detector  practically  inoperable. 
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As  regards  practical  prediction  of  related  soil  influence  on  FDEM  metal  detectors, 
findings  reported  here  generally  support  guidelines  for  soil  magnetic  susceptibility 
adopted  by  CEN-WS7  in  Table  1  of  Figure  62  (CEN,  2006,2008;  Billings,  et  ah, 
2003).  Although  the  relative  influence  of  electrical  conductivity  is  clearly  less  signih- 
cant,  there  is  also  a  basis  in  theory  for  expecting  that  anomalously  high  background 
conductivity  levels  can  too  have  signihcant  influence,  particularly  at  higher  FDEM 
operating  frequencies.  Unfortunately,  there  remains  insufficient  empirical  basis  for 
establishing  corresponding  gradational  guidelines.  However,  theoretical  results  in 
Section  3.1  provide  some  guidance  on  signihcant  electrical  conductivity  levels.  For 
example,  on  the  basis  of  modelled  FDEM  response  characteristics  in  Figure  10  (Won 
and  Huang,  2004),  the  upper  diagram  in  Figure  80  displays  the  induction  number 
for  which  the  quadrature-phase  response  is  equivalent  to  the  in-phase  (resistive  limit) 
response  as  a  function  of  soil  magnetic  susceptibility.  The  corresponding  relation  be¬ 
tween  soil  magnetic  susceptibility  and  electrical  conductivity  is  depicted  in  the  lower 
chart  with  operating  frequency  as  a  parameter. 

It  is  notable  that  for  the  moderate  susceptibility  range  (x=50xl0“^  -  500x10“^  SI) 
set  out  in  CEN-WS7  (Table  1  of  Figure  62),  required  soil  electrical  conductivity  for 
equivalent  quadrature-phase  response  ranges  from  approximately  (t=100  mS/m  to 
roughly  ct=10  S/m,  depending  on  magnetic  susceptibility  and  operating  frequency. 
More  signihcantly,  results  in  Figure  80  conhrm  that  only  rare  soil  environments  (e.g. 
salt-affected  soils)  with  electrical  conductivity  in  excess  of  cj=1.0  S/m  can  have  compa¬ 
rable  severe  inhuence  and  that  even  the  most  conductive  soils  are  virtually  incapable 
of  producing  the  very  severe  effects  of  highly  magnetic  soils  (x  >2000x10“®  SI)  at 
metal  detector  frequencies. 

Similarly,  as  regards  background  soil  influence  on  TDEM  sensors,  analysis  presented 
in  Section  3.3  demonstrates  that  soil  magnetism  is  again  the  primary  cause.  How¬ 
ever,  unlike  the  case  for  FDEM  systems,  it  is  the  frequency  dispersion  of  the  magnetic 
susceptibility  and  equivalent  time-domain  viscosity  that  are  the  key  parameters.  In 
particular,  as  discussed  in  Sections  3.3  and  3.4,  only  hnite  intrinsic  dispersion  (as 
indicated  by  xfd  and  commonly  associated  with  the  colloidal  magnetic  soil  fraction) 
is  required  to  trigger  a  sustained  viscous  response.  The  precise  level  of  viscosity 
and  related  background  response  (noise),  however,  depends  on  the  volume  concentra¬ 
tion  and  mineral  composition  of  dispersive  magnetic  material  comprising  the  soil  as 
reflected  in  the  net  absolute  susceptibility  (in  general,  includes  contributions  due  to 
non-dispersive  as  well  as  dispersive  magnetic  soil  fractions).  Again,  hndings  reported 
here  generally  support  related  guidelines  established  in  CEN-WS7.  However,  interest¬ 
ingly,  it  is  noted  that  specihc  viscosity  {u  =  Ax/  log{fH / fr))  thresholds  specihed  in 

Note  that  results  in  Figure  80  assume  hjv  =  h/o,  =  0.6  and  O  =  0.15  m  for  equiva¬ 
lent  conductivity  calculations.  It  is  emphasized  that  required  electrical  conductivity  for  equivalent 
quadrature-phase  response  increases  with  decreasing  h n  and  with  decreasing  O. 


81 


Table  2  of  Figure  62  (established  with  beneht  of  related  empirical  data  and  together 
with  practical  constraints  on  intrinsic  magnetic  dispersion)  imply  related  susceptibil¬ 
ities  on  the  order  of  y  =  1.0  x  10“^  SI  or  higher  for  even  moderate  influence.  Clearly, 
yet  higher  susceptibilities  would  be  required  for  soils  with  limited  intrinsic  dispersion. 

In  addition,  while  flndings  conflrm  the  uniquely  signiflcant  influence  of  magnetic  vis¬ 
cosity  on  TDEM  sensor  performance,  theoretical  analysis  in  Section  3.3  also  estab¬ 
lishes  the  potential  signiflcance  of  dispersive  soil  electrical  conductivity.  Modelling 
results  summarized  in  Figures  26  and  27  demonstrate  that  appreciable  dispersion 
of  soil  electrical  conductivity  produces  an  anomalous  and  sustained  response  within 
the  typical  measurement  range  of  hand-held  TDEM  systems.  In  particular,  signifl¬ 
cant  electrical  chargeability  yields  a  characteristic  polarity  reversal  and  subsequent 
reduced  decay  rate  with  dependence.  However,  while  the  resulting  signa¬ 

ture  can  potentially  dominate  the  TDEM  response  for  weakly  magnetic  and  strongly 
conductive  soils,  flndings  suggest  that  the  natural  range  of  soil  electrical  conductiv¬ 
ity /chargeability  is  unlikely  to  yield  the  moderate  background  response  level  implied 
by  the  related  magnetic  viscosity  threshold  {u  >  5.0x10-^  SI)  specifled  in  CEN-WS7 
(Table  2,  Figure  62).^^  Indeed,  it  is  anticipated  that  only  abnormally  elevated  soil 
conductivities  (ex  >1.0  S/m)  associated  with  salt-affected  soils  could  yield  a  similar 
response  and  related  polarization  effects  would  likely  occur  at  delay  times  exceeding 

100  ITS. 

Consequently,  the  current  study  generally  conflrms  that  background  noise  levels  con¬ 
templated  by  guidelines  specifled  in  CEN-WS7,  for  both  FDEM  and  TDEM  metal 
detectors,  are  primarily  associated  with  soil  magnetic  properties.  Results  do  suggest, 
however,  that  soil  electrical  conductivity  can  have  comparable  influence  at  moderate 
and  marginally  lower  background  levels  and  could  potentially  constitute  a  signiflcant 
source  of  background  soil  noise  where  soil  magnetic  influence  is  limited. 

Anomalous  Soil  Noise: 

By  anomalous  soil  noise,  we  refer  to  localized  anomalous  response  signatures  or  “clut¬ 
ter”  due  to  small-scale  soil  heterogeneity  or  localized  contrasts  in  larger-scale  soil 

33 

It  is  notable  (and  consistent  with  previous  discussion)  that  viscosity  thresholds  specified  by 
CEN-WS7  guidelines  for  TDEM  sensors  (Table  2,  Figure  62)  imply  very  considerable  absolute 
susceptibility  levels.  For  example,  with  reference  to  Figure  21,  a  strongly  dispersive  magnetic 
soil  (ln(T2/Ti)  =  20,  xfd  10%)  must  also  possess  a  substantial  absolute  susceptibility 
iXdc  =0.001  SI)  to  yield  the  threshold  viscosity  of  P  =  Xdc/ =  5.0  X  10  ^  SI  for 
moderate  effect.  Of  course,  a  soil  having  lesser  dispersion  (e.g.  ln(T2/Ti)  =  100,  XFD  ~  2%) 
requires  still  higher  absolute  susceptibility  (e.g.  Xdc  =0.005  SI)  to  satisfy  the  same  criterion.  Signif¬ 
icantly,  the  same  soil  would  have  severe  influence  on  an  FDEM  metal  detector  according  to  related 
CEN-WS7  guidelines  (Table  1,  Figure  62). 
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electromagnetic  properties  that  are  not  effectively  controlled  by  compensation  cir¬ 
cuitry  or  related  adjustments.  In  contrast  with  background  soil  noise,  anomalous 
noise  does  not  generally  restrict  dynamic  range  or  otherwise  degrade  metal  detector’s 
function.  Rather,  anomalous  noise  is  signal-like  in  character  (at  levels  comparable 
to  weaker  signals)  and  reduces  the  operator’s  target  recognition  and  discrimination 
capacity.  In  general,  while  background  noise  restricts  the  effective  detection  limit 
(probability  of  detection),  anomalous  noise  increases  the  false  alarm  rate.  It  is  also 
possible  for  anomalous  soil  noise  to  mask  lower  level  signals  from  underlying  targets 
on  a  localized  spatial  scale. 

As  a  practical  matter,  the  influence  of  anomalous  soil  noise  is  primarily  significant 
for  low-metal  or  deeply  buried  mines  for  which  related  signal  levels  are  comparable  to 
soil  generated  response  variation.  Unfortunately,  expected  signal  levels  and  related 
detection  limits  are  also  substantially  dependent  on  the  operating  characteristics  of  a 
given  metal  detector  and,  additionally,  on  background  soil  properties.  Consequently, 
it  is  difficult  to  establish  broadly  applicable  guidelines  (e.g.  CEN-WS7  provides  no 
related  guidance  comparable  to  Tables  1  and  2,  Figure  62)  for  influence  of  anomalous 
soil  noise  as  functions  of  specific  soil  electromagnetic  parameters.  Nevertheless,  it 
would  appear  reasonable  to  anticipate  that  a  substantially  wider  range  of  soil  electro¬ 
magnetic  properties  (magnetic  and  electrical)  could  be  signihcant. 

In  general,  additional  studies,  both  laboratory  and  held  are  required  to  confirm  ex¬ 
isting  guidelines  for  performance  prediction  and  to  further  investigate  and  establish 
the  practical  influence  of  soil  electromagnetic  parameters.  As  possible  and  through 
whatever  agencies  or  organizations  might  facilitate  or  sponsor  it,  a  sustained  cam¬ 
paign  should  be  pursued  to  encourage  the  measurement,  sampling  and  reporting  of 
soil  electromagnetic  properties  in  connection  with  humanitarian  demining  operations, 
particularly  where  soil  conditions  are  perceived  to  be  problematic.  Compilation  and 
analysis  of  resulting  data,  together  with  results  of  detector-based  performance  eval¬ 
uation  (as  per  CEN-WS7  recommendations)  would  further  validate  and  advance  our 
understanding  of  soil  influence  on  metal  detector  performance  and  could  potentially 
lead  to  significant  technological  innovation.  We  trust  that  the  foregoing  review  of  rel¬ 
evant  theory  and  measurement  methodology  provides  a  useful  foundation  for  related 
future  research  and  development. 
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NOTE:  Coincident  horizontal  Tx/Rx  (a  =  b  =  0.1  m)  at  soil  surface. 

v1  -  non-magnetic,  non-dispersive  conductive  soil 
v2  -  non-magnetic,  non-dispersive  conductive  soil  (late  time) 
v3  -  non-dispersive  magnetic,  non-dispersive  conductive  soil 
v4  -  non-conductive,  dispersive  magnetic  soil 
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NOTE:  Coincident  horizontal  Tx/Rx  (a  =  b  =  0.1  m)  at  soil  surface. 

v1  -  non-magnetic,  non-dispersive  conductive  soil 
v2  -  non-magnetic,  non-dispersive  conductive  soil  (late  time) 
v3  -  non-dispersive  magnetic,  non-dispersive  conductive  soil 
v4  -  non-conductive,  dispersive  magnetic  soil 
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NOTE:  Coincident  horizontal  Tx/Rx  (a  =  b  =  0.1  m)  at  soil  surface. 

v1  -  non-magnetic,  non-dispersive  conductive  soil 

v2  -  non-magnetic,  non-dispersive  conductive  soil  (late  time) 

v3  -  non-dispersive  magnetic,  non-dispersive  conductive  soil 

v4  -  dispersive  magnetic,  non-dispersive  conductive  soil 

v5  -  dispersive  magnetic,  non-dispersive  conductive  soil  (approximation) 
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NOTE:  Coincident  horizontal  Tx/Rx  (a  =  b  =  0.1  m)  at  soil  surface. 

v1  -  non-magnetic,  non-dispersive  conductive  soil 
v2  -  non-magnetic,  non-dispersive  conductive  soil  (late  time) 
v3  -  non-dispersive  magnetic,  non-dispersive  conductive  soil 
v4  -  dispersive  magnetic,  non-dispersive  conductive  soil 

v6-p  -  non-magnetic,  dispersive  conductive  soil  (positive-polarity  -  terminated  source) 
v6-n  -  non-magnetic,  dispersive  conductive  soil  (negative-polarity  -  terminated  source) 
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NOTE:  Coincident  horizontal  Tx/Rx  (a  =  b  =  0.1  m)  at  soil  surface. 

v1  -  non-magnetic,  non-dispersive  conductive  soil 
v2  -  non-magnetic,  non-dispersive  conductive  soil  (late  time) 
v3  -  non-dispersive  magnetic,  non-dispersive  conductive  soil 
v4  -  non-conductive,  dispersive  magnetic  soil 

v6-p  -  non-magnetic,  dispersive  conductive  soil  (positive-polarity  -  terminated  source) 
v6-n  -  non-maanetic.  disoersive  conductive  soil  (neaative-oolaritv  -  terminated  source) 
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NOTE:  Coincident  horizontal  Tx/Rx  (a  =  b  =  0.1  m)  at  soil  surface. 

Effective  Decay  Rate 

v6-n  -  non-magnetic,  dispersive  conductive  soil  (negative-polarity  -  terminated  source) 
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Series  Impedance  (Z=V/I) 


Complex  Impedance:  Z 
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Z=R+i  [27tfL-(27tf  Cr  V  0  =  tan-^X/R) 

6  =  tan"'(R/X) 
Complex  Admittance:  Y=Z“^=(  R+iX)“^ 
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SOIL  BOX-  PHYSICAL  DATA* 

SMALL  LARGE 


Length  (L-cm) 
Width  (W-cm) 
Height  (H-cm) 
XI  (cm) 

X2  (cm) 

Mass  (g) 
Volume  (ml) 


Body 

Pin  Electrodes 
Plate  Electrodes 


22.0 
3.8 
3.0 
12.4 
4.6 
240 
250 

MATERIALS 

Lexan 

Brass 

Stainless  Steel 


11.0 

3.0 

2.4 

7.1 

1.8 

130 

75 


TOP  VIEW 


END  VIEW 


*  NOTE:  Physical  data  are  nominal  values  -  actual  data  deviate  within  approximately  +/-  0.1  cm 


FIGURE: 

^^KTerrascan 

Geophysics 

MC  MILLER  SOIL  RESISTIVITY  BOX 

34 

PROJECT: 

DRDC  Soils 

DRAWN  BY: 

GMC 

DATE: 

14  December,  2003 

— >  Primary  Field 

Polarization  Field 


Xy 


Terrascan 

Geophysics 


FOUR  ELECTRODE  MEASUREMENT 

FIGURE: 

ELECTRODE  POLARIZATION 

35 

PROJECT: 

DRAWN  BY: 

DATE: 

CCM AT  Soils 

GMC 

15  May,  2007 

□□□□ 


Lc  Lp  Hp  He 


FOUR  TERMINAL 
TWO  ELECTRODE 


Lc  Lp  Hp  He 


FOUR  TERMINAL 
FOUR  ELECTRODE 


SOIL  BOX  RESISTIVITY  MEASUREMENT  CONFIGURATIONS 
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SOIL  BOX  RESISTIVITY  MEASUREMENT  CONFIGURATIONS 
GENERATOR  VOLTAGE  CURRENT  (GVI )  WITH  ISOLATION 
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Resistivity  (Ohm-m) 


O  2-electrode  GVI  •  2-electrode  GVI  -  iso 
□  4-electrode  GVI  ■4-electrode  GVI  -  iso 
A2-electrode  HP4192A  A4-electrode  HP4192A 
O  2-electrode  HP4274A 


Frequency  (kHz) 
75  ml  Soil  Box 


O 2-electrode  GVI  •2-electrode  GVI  -  iso 

□  4-electrode  GVI  ■4-electrode  GVI  -  iso 

A  2-electrode  HP4192A  A  4-electrode  HP4192A 
O  2-electrode  HP4274A 


Frequency  (kHz) 
250  ml  Soil  Box 
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□  4-electrode  GVI  ■4-electrode  GVI  -  iso 

A  2-electrode  HP4192A  A  4-electrode  HP4192A 
O  2-electrode  HP4274A 


o 

o 

3 

Q. 

C 

o 

< 

'■< 

3 

(f) 

3 


Frequency  (kHz) 


Frequency  (kHz) 


250  ml  Soil  Box  Vancouver  Tap  Water 


250  ml  Soil  Box  Vancouver  Loam  Soil 


PARAMETERS:  w  -  gravimetric  noisture  content  (refered  to  "lab-dry”  condition. 

m2E  ■  electrical  chargeability  (2-electrode)  HP4192A  - 10  Hz,  40  kHz. 
m4E  -  electrical  chargeability  (4-electrode)  GVI  - 10  Hz,  40  kHz. 
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Frequency  (kHz) 


250  ml  Soil  Box 


PARAMETERS:  w  -  gravimetric  noisture  content  (refered  to  "lab-dry”  condition. 

m2E  -  electrical  chargeability  (2-electrode)  HP4192A  - 10  Hz,  40  kHz. 
m4E  -  electrical  chargeability  (4-electrode)  GVI  - 10  Hz,  40  kHz. 
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Frequency  (kHz)  Frequency  (kHz) 


250  ml  Soil  Box  250  ml  Soil  Box 


PARAMETERS:  w  -  gravimetric  noisture  content  (refered  to  "lab-dry”  condition. 

m2E  ■  electrical  chargeability  (2-electrode)  HP4192A  - 10  Hz,  40  kHz. 
m4E  -  electrical  chargeability  (4-electrode)  GVI  - 10  Hz,  40  kHz. 
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Frequency  (kHz) 
250  ml  Soil  Box 


Frequency  (kHz) 
250  ml  Soil  Box 


PARAMETERS:  w  -  gravimetric  noisture  content  (refered  to  "lab-dry”  condition. 

m2E  -  electrical  chargeability  (2-electrode)  HP4192A  - 10  Hz,  40  kHz. 
m4E  -  electrical  chargeability  (4-electrode)  GVI  - 10  Hz,  40  kHz. 
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GENERATOR  VOLTAGE  CURRENT  (GVI ) 
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IN  SITU  ELECTRICAL  RESISTIVITY  MEASUREMENT 
SUBSURFACE  POTENTIAL  /  CURRENT 
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SOURCE:  Dahlin  and  Zhou,  2004 
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Source:  Won  and  Huang  (2004) 
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Normalized  Response 


PROTOTYPE  SOLENOID  FIXTURE 


Coil  Parameters: 
N=110 
R=0.016m 
1=0. 11m 

Model  Parameters: 
ls=8.5  mA 

l.ir=1.0 
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Mass  Susceptibility 


CRC  Reference 


Prototype  Single  Coil  Fixture 
HP4274A  Direct  Inductance  Susceptibility  Analysis 
MnCOs  Standard 

1.500 

1.400 

1.300 

S  1.200  ' 

1.100 
^  1.000 
°  0.900 

i  0.800 

0.700 

0.600 

0.500 

1.E+02  1.E+03  1.E+04  1.E+05 


Mass  Specific  Susceptibility 
1.2463  X 10'^  m^/kg 

Fukuma  -  Bartinaton  MS2 

Mass  Specific  Susceptibility 
1.1418  xia^  m^/kg 

Fukuma  -  Kappabridue 

Mass  Specific  Susceptibility 
1.2472  X  10'^  m^/kg 


Frequency (Hz) 


Terrascan  -  HP4274A 


NOTES: 

1.  Dashed  line  indicates  reference  value  as  tabulated  CRC  Handbook  of 
Chemistry  and  Physics,  51st  Edition  (1971). 

2.  Bartington  MS2  and  Kappabridge  data  reported  by  K.  Fukuma, 
Department  of  Environmental  Sciences,  Kumamoto  University,  Japan. 

3.  Mean  value  measured  via  HP4274A  yields  fixture  calibration  constant 
Kd=1.053 


Mass  Specific  Susceptibility 
Mean  (100  Hz -100  kHz) 
1.1834  X 10'^  m^/kg 
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Mass  Susceptibility 


CRC  Reference 


Prototype  Single  Coil  Fixture 
HP4274A  Direct  Inductance  Susceptibility  Analysis 
Gd203  Standard 


Mass  Specific  Susceptibiiity 
1.8442  X  10'^  m^/kg 

DRDCS  -  Bartinaton  MS2 

Mass  Specific  Susceptibiiity 
1.5787  X  10'^  mVkg 

Sagnotti,  et  al.  -  Bartington 


Mass  Specific  Susceptibiiity 
Mean  (Table  2) 

1.6651  xia^  m^/kg 


Sagnotti  et  al.  -  Kappabridge 


Frequency (Hz) 

NOTES: 

1.  Dashed  line  indicates  reference  value  as  tabulated  CRC  Handbook  of 
Chemistry  and  Physics,  51st  Edition  (1971). 

2.  Bartington  MS2  and  Kappabridge  data  reported  by  L.  Sagnotti  et  al., 

2003,  Physics  of  the  Earth  and  Planetary  Interiors,  138,  25-38. 

3.  Mean  value  measured  via  HP4274A  (o)  adjusted  (m)  via  fixture  calibration 
constant  Kd=1.053  determined  for  MnC03  standard. 


Mass  Specific  Susceptibiiity 
Mean  (Table  2) 

1.7437  xlQ-^  m^/kg 

Terrascan  -  HP4274A 

Mass  Specific  Susceptibiiity 
Mean  (100  Hz -100  kHz) 

O  1.5648  X  10'^  m^/kg 
•  1.6477  X  10'^  /kg 


-®-  ZHI  SM-100(35mm)  |  p  *  0.99  g/cm^  |  Ax  =  479x10'®  |  xfd=  10.5  % 

-O-  ZHI  SM-IOO(IOcc)  I  p  *  1 .02  g/cm®  |  Ax  =  491x10'®  |  xfd=  10.3  % 

Bartington  MS2B(10cc)  |  p*1.02g/cm®  |  Ax  =  534x10'®  |  xfd=  1 1 -5  % 

-0-  ESI  2150  LCR  I  p  *  0.91  g/cm®  |  Ax  =  427x10'®  |  xfd=11.1% 

Solenoid  fixture 


100 


1000  10000 
Frequency  (Hz) 


100000 


DRDC  BOSNIA  REFERENCE  SOIL 

FIGURE: 

fiv  Terrascan 
Geophysics 

MASS-SPECIFIC  MAGNETIC  SUSCEPTIBILITY 

59 

PROJECT: 

COM  AT  Soils 

DRAWN  BY: 

GMC 

DATE: 

15  February,  2007 
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-©-  ZHI  SM-100(35mm)  |  p«1.46g/cm^  |  Ajc  =  210x10'®  |  xfd=  1-1  % 

-O-  ZHI  SM-100(10cc)  I  p*1.53g/cm^  |  Ax  =  232x10'®  |  xfd=  13  % 

— Bartington  MS2B  (10  cc)  |  p  a  1.53  g/cm®  |  Ax  =  319x10'®|  xfd=1-8% 
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Table  1  —  Indicative  values  of  magnetic  susceptibility  that  can  be  expected  for  soil  effects 
for  single  frequency  continuous  wave  metal  detectors,  source:  cwa  14747-1,  A3 


Soil  effect  class 
(Defined  in  Error!  Reference 
source  not  found.) 

Indicative  values  of  magnetic  susceptibility 

10  ®  SI 

Neutral 

Below  50 

Moderate 

50  to  500 

Severe 

500  to  2  000 

Very  severe 

Above  2  000 

NOTE  These  values  are  indicative  oniy  because  the  effect  of  a  soil  on  a  metal  detector  depends  on  the  metal 
detector,  and  the  measured  values  of  susceptibility  depend  on  the  measuring  instrument. 

Table2  —  Indicative  values  of  frequency  variation  magnetic  susceptibility  that  can  be  expected  for  soil 
effects  for  pulsed  induction  and  most  continuous  wave  metal  detectors.  SOURCE:  Billings,  etai.  (2003) 


Soil  effect  class 
(Defined  in  Error! 
Reference  source  not 
found.) 

Indicative  values  of  frequency  variation  of  magnetic  susceptibility 
(458  Hz  and  4  580  Hz) 

10  ®  SI 

Neutral 

Below  5 

Moderate 

5  to  15 

Severe 

15  to  25 

Very  severe 

Above  25 

NOTE  These  values  are  indicative  only  because  the  effect  of  a  soil  on  a  metal  detector  depends  on  the  metal 
detector,  and  the  measured  values  of  susceptibility  depend  on  the  measuring  instrument. 

Table  A.1  —  Definition  of  classes  for  soil  effects  on  metal  detectors 


Neutral 

A  soil  has  a  neutral  effect  on  a  metal  detector  response  if  it  has  a  no  influence 
on  the  performance  of  the  metal  detector  even  without  ground  compensation. 

For  the  detector,  such  a  soil  is  equivalent  to  air. 

Moderate 

A  soil  has  a  moderate  effect  on  a  metal  detector  response  if  its  influence  on  the 
detector  performance  is  noticeable  but  the  metal  detector  can  be  used  without 
ground  compensation. 

Severe 

A  soil  has  a  severe  effect  on  a  metal  detector  response  if  it  makes  the  use  of 
ground  compensation  necessary. 

Very  severe 

A  soil  has  a  very  severe  effect  on  a  metal  detector  response  if  the  metal 
detector  cannot  be  used  even  with  ground  compensation. 
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NOTES: 

1.  Bottom  12  mm  spacer  comprised  of  dense  non-magnetic  foam. 

2.  Standard  10  cc  sample  vial  lid  centralizes  vial  in  35  mm  holder. 

3.  Adapter  specific  fixture  constant  required  for  direct-inductance 
susceptibility  measurements  with  10  cc  -  35  mm  adapter. 
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ZHISM-100(35  mm)  |  p  0.91  g/cm"  |  Ax  =  2465x10'”  |  xfd  =  23.9  % 

ZHI  SM-1 00(10  cc)  I  p*0.84g/cm^  |  Ax  =  2270x10'®  |  xfd- 23.9% 

BartingtonMS2B(10cc-1.0)  |  p  *  0.91  g/cm®  |  Ax  =  2587x10'®  |  xfd  =  25.5  % 

Bartington  l\/IS2B  (10  cc-0.1)  |  p  *  0.92  g/cm®  |  Ax  =  2468x10'®  |  xfd  =  24.3  % 

ESI  2600  LCR  |  p^0.92g/cm®  |  Ax  =  2378x10'®  |  xfd  =  22.9% 

35  mm  fixture 
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y  =  1 2. 333x  -2.0933 


STANDARD 
10  cc  POT 


NON-STANDARD: 
Height  =  4.9  cm 
Diameter  =  2.1  cm 
Volume  =  16.97  cc 

NOTE: 

Volume  at  2.3  cm 
Sample  level  =  7.97  cc 
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Kappaneter  KT-6  Magnetic  Susceptibility  Meter 
Operating  Frequency:  10  kHz 

Sensitivity:  IxlO'"’  SI 


Bartington  MS2D  Field  Magnetic  Susceptibility  System 
Operating  Frequency:  0.958  kHz 
Sensitivity:  2x10'^  SI 
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